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An Investigation of the Effect of a Varying 
Tip-Weight Distribution on the Flutter 
Characteristics of a Straight Wing 


* 


WILLIAM H. GAYMANt 
Northrop Aircraft, Inc. 


ABSTRACT 


Experimental and theoretical investigations into the antisym- 
metric flutter characteristics of a straight wing having a number 
of tip-pod mass, size, and shape configurations are summarized. 
The experimental phase, conducted with a low-speed flutter 
model of the U.S.A.F.’s Northrop Scorpion F-89 wing, is discus- 
sed from the standpoint of a design development study. Model 
design features, testing techniques, and test results are described, 
with emphasis on an observed need for considering small varia- 
tions of mass distribution within a pod having expendable con 
tent. 

Methods and results of theoretical flutter analyses of several 
tested pod configurations are discussed. Computations from 
parallel analyses using different assumptions of pod aerodynamic 
eflectiveness are compared with test data. 

It is concluded from these investigations that the aerodynami- 
tally effective tip pod may be small in comparison with the pod 
geometric plan form but, further, that an accurate representation 
of this effective tip pod may be essential to a successful flutter 


Received September 5, 1951. 

* The investigation described herein is the product of a team 
Peffort. Although it is not possible to acknowledge all individual 
contributions, particular recognition is due personnel of the Air 
craft Laboratory, Wright Air Development Center. The basic 
model was designed and instrumented by L. A. Tolve and R. P 
Peloubet and was constructed by F. F. Hense. These and other 
Aircraft Laboratory personnel, notably W. J. Mykytow, are 
thanked for their generous cooperation with Northrop personnel 
during the transition of test operations. Within the Dynamics 
Section of Northrop Aircraft, Inc., the major part of the theoreti 
| cal investigation has been carried out by J. K. Latzer, who also 
has conceived the inverse solution of the flutter determinant as a 
possible means of establishing the aerodynamic effectiveness of a 
tip pod from measured flutter data. D. R. Branchflower has 
made this approach practical from a computational standpoint 
by a problem formulation amenable to solution by IBM comput- 
ing facilities. D. T. Egbert has directed the course of much of the 
theoretical work Finally, W. E. Cox, in his supervision of all 
Phases of the Northrop-conducted investigations, has given in- 
Valuable aid through his many suggestions and his encourage- 
ment. 

t Dynamics Engineer. 


analysis. It is also indicated that analyses of wing-and-pod con- 
figurations having the pod center of gravity on or near the wing 
elastic axis are more adversely affected by errors in pod aerody- 
namic representation than are those analyses applying to cases of 


large pod overbalance. 


NOTATION 


wing-element weight, Ibs. 
wing-element unbalance about wing elastic axis, in. 
Ibs. 
wing-element moment of inertia about wing elastic 
axis, lb. in.? 
tip-pod weight, lbs. 
tip-pod unbalance about wing elastic axis, in.lbs 
tip-pod moment of inertia about wing elastic axis, 
Ib. in.? 
wing flexural rigidity, lb. in.? 
wing torsional rigidity, lb. in.? 
frequency of first antisymmetric uncoupled bending 
mode, cycles per min. 
frequency of first antisymmetric uncoupled torsion 
mode, cycles per min 
frequencies of first two antisymmetric coupled modes 
normalized wing deflection in the rolling degree of 
freedom 
normalized wing deflection in the first antisymmetric 
uncoupled bending mode 
a reference h 
normalized wing deflection in the first antisymmetric 
uncoupled torsion mode 
normalized wing deflection in the second antisym- 
metric uncoupled torsion mode 
flutter-mode damping coefficient 
flutter frequency, cycles per min 
critical flutter speed, m.p.h 
wing semichord, in. 
tip-pod semichord, in. 
tip-pod span, in. 
a} 
location of the wing 
the lifting-surface quarter-chord, as a fraction of 
when aft of the 


elastic axis measured from 


the local semichord; positive 


quarter-chord 
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INTRODUCTION 


K ibenigenrgin FLUTTER ANALYSES of wings carrying 
tip pods have not been consistently successful in 
giving reliable predictions of flutter speed. A few cases 
are known wherein, contrary to expectations based on 
analyses, aircraft with tip fuel tanks have experienced 
wing flutter. More numerous instances are known of 
disagreement between results of analysis and of flutter- 
model tests, in which the incidence of flutter has been 
observed much more often than in flight. 

The inadequacy of analytical treatment is believed 
to lie primarily in the absence of a rational accounting 
for the aerodynamic forces on a tip pod oscillating in 
translation and in pitch. Furthermore, in the case of 
tip tanks, the effective pitching moment of inertia of 
the fuel at the flutter condition is not predictable to 
the degree of accuracy that may be demanded in a given 
situation. Indeed, it is in recognition of these and 
other related limitations in the flutter analysis of wings 
with tip pods that the Aircraft Research and Testing 
Committee of the Aircraft Industries Association is 
favoring a research program aimed at meeting the spe- 
cific needs of the problem. 

Some of the spade work has already been done by 
various investigators. The N.A.C.A., for example, has 
conducted both theoretical and experimental investiga- 
tions into the flutter characteristics of straight uniform 
cantilevered wings having variable placement of a con- 
centrated mass. However, it is not attempted here to 
review the progress in this particular field, nor is it 
attempted to assess the relation between other investi- 
gations and that of the subject discussion. 

About a year ago the Dynamics Branch of the Wright 
Air Development Center's Aircraft Laboratory initiated 
a series of wind-tunnei tests of a flutter model of the F- 
89 wing, which has an integral tip fuel pod. Several 
months later Northrop Aircraft acquired the model for 
use in a tip-pod design development study. To date 
about 300 pod mass configurations have been tested 
among six tip pods differing in size and shape. The pur- 
pose of this paper is to summarize the salient features 
and results of these test programs and to report on the 
progress of a complementary analytical approach to 
this flutter problem. 


SCOPE OF THE INVESTIGATION 


The flutter-model test programs described herein 
were begun by the W.A.D.C. with a specific objective 
of determining the flutter safety of an existing wing- 
aileron design. This program was extended by the 
W.A.D.C. to include preliminary series of tests of tip- 
pod installations. These tests, and other tip-pod in- 
vestigations subsequently conducted by Northrop Air 
craft, have had a similar motivation, with the result 
that the entire experimental phase has assumed the as 
pect of a design development study rather than a re- 


search program. 
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This philosophy of approach is emphasized because 
of its dominant influence on flutter-model construction, 
on the establishment of test configurations, and, for 
reasons to be discussed presently, on the quest for ac- 
ceptable methods of theoretical analysis. 

At the outset it was considered that preliminary coy 
firmation of the flutter safety of specific wing-and-pod 
design configurations should rest entirely on model 
tests, with final confirmation of the safety of any chosen 
configuration depending on flight flutter testing judi. 
ciously based on model experience. Therefore, the basic 
requirements to be met in a flutter-model test program 
so directed are as follows: (1) attainment of dynamic 
similarity between the model and its prototype and 
(2) coverage of critical weight configurations of wing 
and tip pod. 

In general, the critical weight régimes cannot be 
foreordained. Therefore the test objectives, shifting as 
the program progresses, may pursue such a chronologi 
cal course as follows: 

(1) The initial objective is the definition of flutter 
characteristics as functions of specific wing-and-pod 
weight configurations that follow the expendable-load 
sequencing to be preferred from aircraft design and 
operational considerations. 

(2) In the event that an unsafe flutter condition is 
indicated by the foregoing, avoidance of major design 
changes might be sought by any of the following expedi- 
ents: (a) repetition of (1) with alternate usage se- 
quences that are feasible from a design standpoint and, 
at the same time, that do not penalize the operational 
utility of the aircraft; or (b) investigation of the ef- 
fects of incorporating ballast weights in the pod; or (c 
imposition of tentative aircraft operational restric- 
tions, to be evaluated by subsequent flight flutter 
testing, provided that such limitations, again, constitute 
no serious sacrifice to the operational utility. 

(3) If none of the alternatives of (2) is acceptable 
from a design or a flutter standpoint, further testing 
may be directed toward examination of more drasti 
remedial measures, such as a change in the pod-wing 
geometric configuration, addition of pod fins, relocation 
of pod fuel bulkheads, added or revised fuel compart 
mentization within the pod, or elastic mounting of the 
pod on the wing tip. 

The experimental investigations to be described have 
pursued this chain of objectives to the extent required 
with each of several proposed design configurations. 
However, certain limitations in the scope of the tests 
have been imposed by the practical considerations 0! 
model design and of available test facilities. These 
considerations have dictated the choice of a low-speed 
flutter model, leaving the evaluation of high Mach 
Number effects on flutter characteristics to a delicate 
blend of theory, experience, and judgment. Moreover, 
use of a single model in constant-density wind tunnels 
has precluded investigation of altitude as a flutter 








pa 


to 


nal 
are 
the 
apy 
the 
at 

bet 
mo 
fre 
tio1 


for 


mo 


to | 
suc 
ma 
tiot 
act 
con 
tra 
ant 
ade 
wit 
frec 


fab 
ent 
all 


ribs 





because 
ruction, 
ind, fc ir 


for ac- 


iry con- 
ind-pod 

mode] 
Chosen 
ig judi- 
he basic 
rogram 
ynamic 
pe and 
1 wing 


not be 
ting as 
nologi- 


flutter 
nd-pod 
le-load 
n and 


tion is 
design 
xpedi- 
ge se- 
t and, 
tional 
he ef- 
or (Cc 
estric- 
utter 
titute 


table 
esting 
rastic 
-wing 
‘ation 
ipart 
f the 


have 
uired 
tions. 

tests 
ns of 
[hese 
speed 
\Iach 
licate 
over, 
nnels 
utter 





FLUTTER CHARACTERISTICS OF 


parameter. Therefore the test results apply directly 


to only a subsonic sea-level environment. 


THE EXPERIMENTAL INVESTIGATION 


Description of the Flutter Model 


Wing.—The design criteria for attainment of dy- 
namic similarity between a low-speed flutter model and 
its equivalent prototype are rather basic and therefore 
are not restated here. It is noted, however, that in 
the design of the F-89 wing flutter model, scale factors 
applying to geometry and to air speed were selected by 
the W.A.D.C. to permit use of the 5-Ft. Wind Tunnel 
at Wright Field. These choices fixed the relationship 
between the frequencies of corresponding oscillation 
modes of model and prototype. Attainment of the 
frequency scale factor was sought by appropriate reduc- 
tions of airplane-wing mass and rigidity distributions. 


The F-89 wing flutter model is a semiwing, designed 
for vertical support from a faired structure that may be 
mounted against the ceiling of a wind-tunnel working 
section, as shown in Fig. 1. The wing root is attached 
toa hinged plate provided with a locking arrangement 
such that, optionally, the wing may be cantilevered or 
may be freed in roll. This provision permits investiga- 
tion of both symmetric and antisymmetric flutter char- 
acteristics, although it excludes from the former those 
contributions by the rigid-body symmetric modes of 
translation and pitch that exist in free flight. For the 
antisymmetric case, however, free-flight simulation is 
adequate, inasmuch as the rolling frequency of the 
wing, as a pendulous rigid body, is much lower than the 
frequency of the first elastic mode. 

The basic model semiwing consists of nine balsa-and- 
fabric-covered skeletal wooden sections, each independ- 
ently attached at its bounding ribs, to an aluminum- 
alloy spar. The covering between the closely spaced 
ribs of adjacent sections has been cut around the entire 
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Installation of flutter model in Convair Wind Tunnel 


Fic. 1. 


airfoil, and the resulting gaps are sealed by thin latex 
strips. This type of construction has been used to 
minimize contributions of the wooden sections to the 
flexural and torsional rigidity distributions of the as- 
sembled wing. Thereby, the structural design of the 
wing was simplified considerably; the centerline of the 
spar was located to approximate the elastic axis location 
of the prototype wing, and the spar was designed to pro- 


The W.A.D.C. 


conducted static tests on the assembled wing in order to 


vide nearly all of the required rigidity. 


TABLE | 
Model-Wing Data 


a) Geometric: 
Semispan (less tip pod ) 
Aspect ratio 
Taper ratio 


0.51 
NACA 0009-64 





Airfoil : 
Sweepback of 30 per cent chord: 0 
b) Inertial and Structural: 
Wing Section 
Midspan M Sa la EI X 10 GJ X 10 
Number station (Lbs. ) (In. Lbs (Lb. In.? (Lb. In.? Lb. In.? 
l 2:0 0.679 0.872 10.3 0.98 0.98 
2 6.0 1.482 2.038 27.9 0.98 0.98 
3 10.0 1.511 l 332 18.4 0.98 0.98 
4 14.0 1.642 2.198 14.9 0.90 0.89 
5 18.0 1.484 1.922  - 0.76 0.77 
6 22.0 1.217 1.773 14.9 0.61 0.59 
7 26.0 1.268 1.930 14.5 0.40 0.40 
8 30.0 1.022 1.192 9.1 0.23 0.22 
g a.7 0.377 0.433 3.7 0 1] 0 O8 
10.682 13.690 134.8 
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Profiles of tested pod configurations relative to the wing 
tip chord. 


Fic. 2. 


indicate any final adjustments needed to realize ade- 
quate representation of design // and GJ distributions. 

tach of the nine wing sections was ballasted with lead 
weights in a manner to simulate a lumped distribution 
of prototype-wing weight, static unbalance about the 
elastic axis, and moment of inertia about the elastic 
axis. The mass parameters of each section were de- 
termined experimentally prior to wing assembly in 
order to assure reasonable agreement with design 
values. 

Geometric, inertial, and structural data applying to 
the model wing are given in Table 1. The mass dis- 
tribution shown corresponds to a condition with full 
wing fuel tanks; other distributions, not listed, were 
also used. 

The model was fitted with an aileron, which was 
locked in its neutral position during the subject series 
of tests. 

Tip Pods. the 
tested tip pods are shown in relation to the tip wing sec- 
tion in Fig. 2. Pods V and VI were also tested in 
slightly forward positions, the displacements being 
about 18 per cent tip chord ahead of the positions 


The geometric configurations of 


shown. 

Pod construction details were dictated by a consider- 
ation of test loading schedules simulating fuel usage. 
Several problems were presented. First, in a partially 
filled tank of an aircraft in flight, the instantaneous 
disposition of fuel is determined by a time history of 
airplane attitude and acceleration. Even if the dy- 
namic effects of fuel sloshing are to be excluded from a 
flutter investigation, there still remain an infinite num- 
ber of instantaneous configurations, several régimes of 
which may be of concern from a flutter standpoint. 
Therefore, it was sought to simulate a practical number 
of discrete configurations representing a liberal sam- 
pling of those possible. 

The next problem was concerned with the means of 
such simulation. General use of a liquid was precluded 
by intolerable complications in model construction, 
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inasmuch as the wind-tunnel working section could no 
be rotated in the gravity field. 
variable quantity of liquid in the nose or in the tail of g 


For example, use of a 


tip tank would require a rather elaborate system of bulk. 
heads that were changeable during the course of the 
tests. Thus, as a general expedient, approximations 
of desired fuel conditions were effected by use of suit. 
ably disposed lead slugs that, in all but one of the pods 
so loaded, had centroidal pitching moments of inert 


representing solid or ‘‘frozen’’ fuel. 


Use of frozen fuel in flutter tests introduces errors 
How- 
ever, on the assumption that tip-pod-fuel sloshing fre. 


that are not amenable to theoretical analysis. 


quencies will generally be far below any flutter fre. 
quency, the existence of a significant time-varying fuel 
inertia can be dismissed, together with the added de 
gree of freedom of fuel motion. Therefore, the frozep. 
fuel errors may be assessed qualitatively as a conserva 
tive omission of energy dissipation within the pod 
and, of potentially greater importance, as an improper 
representation of pod moment of inertia. No general 
criteria have been evolved for determining whether 
such representation is conservative or unconservative 
from a flutter standpoint. However, in the correlation 
of data obtained in the subject tests, it has generally 
been possible to indicate whether the use of liquid fuel 
would raise or lower the flutter speed observed with 
frozen fuel. Such indications are predicated on a neg- 
lect of added damping and on a reduction of pod pitch- 
ing moment of inertia. The order of magnitude of the 
reduction is bracketed by pod inertia values calculated 
on the extreme bases of solid fuel and of zero centroidal 
moment of inertia of the fuel within any one compart 
ment. 

Pods I and II represent airplane tip pods having 
single-compartment fuel tanks. For these tank shapes, 
errors accruing from the frozen fuel assumption may be 
expected to be large. In order to permit a direct evalu- 
ation of these errors as applying to a particular type ol 
fuel disposition, the W.A.D.C. constructed a thin fiber 
glas model of Pod I having internal bulkheads corre- 
sponding with those of the prototype fuel tank. A 
liquid having the proper specific gravity was used for 
the loading medium. Pod II was made by attaching a 
paper tail cone to Pod I and ballasting as required to 
meet the “‘empty-pod”’ condition. 


Model pods of each configuration, except that oi 
Pod II, were made of balsa hollowed out and fitted with 
provisions for attaching the lead weights intended to 
simulate expendable load. For each pod, design mass 
parameters—weight, c.g. location, and centroidal pitch 
ing moment of inertia—were calculated for all chosen 
test configurations, and the pod ballasting arrangement 
was effected accordingly. Actual mass parameters of a 
representative number of pod configurations were de- 
termined experimentally to ascertain reasonable agree- 


ment with design values. Table 2 lists the model-pod 
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TABLE 2 
Tip-Pod Mass Parameters for ‘‘Empty”’ and ‘‘Full’’ Conditions 


—‘‘Empty”’ Condition 


Sa, 


Pod M7 T 

Designation (Ibs. ) (in.Ibs. ) 

I 0.261 —(0.211 
ll 0.457 0.765 
III 1.189 —2.154 
IV 1.825 —2.353 
V (aft position ) 1.835 —4.,244 
(forward position ) 1.835 —6.630 

VI (aft position ) 2.495 —6.509 
(forward position ) 2.495 —10.25 


“Full” Condition 


Ta, M, Sa, la T 
(Ib. in.?) (Ibs. ) (in.Ibs. ) (Ib. in.?) 
7.35 3.364 —2.677 69.77 
13.67 3.560 —1.701 — 
37.62 5.285 —2.166 144.5 
68.83 6.384 0.9838 219.7 
71.60 6.412 —2.113 230.7 
85.74 6.412 — 10.45 247.0 
128.1 7.072 —6.028 316.3 
153.2 7.072 — 16.67 350.5 


** The moment of inertia of Pod II, filled with a liquid, was not measured. 


mass parameters corresponding to ‘‘empty”’ and “‘full’’ 
conditions. 

Pod attachment to the wing was made by three ma- 
chine screws, one entering a tapped hole in the end of 
the spar and the other two engaging plate nuts at for- 
ward and aft locations on the tip rib. Pillow blocks 
placed between the pod and the rib served to stiffen the 
mounting in roll. Quantitative representation of the 
rigidity of the prototype pod attachment was not at- 
tempted; it was merely noted that such mounting 
would be sufficiently rigid to obviate the need of the 
added degree of freedom of pod motion relative to the 
wing tip. 

Pods III through VI represent pod designs having 
several fuel compartments, and, therefore, frozen fuel 
errors could be shown to be less important than those 
applying to Pods I and II; _ the importance is governed 
to some extent by fuel sequencing. 

The simulation of fuel in one of the tested pods was 
based on measured moments of inertia of a model cell 
containing water. This cell was a right cylinder with 
length equal to radius and with the rest position of its 
axis horizontal. Moments of inertia of the cell in 
several conditions between empty and full were meas- 
ured about a horizontal centroidal diameter of the cell 
by the method of timed oscillations. The test results, 
as applying to the water only, are shown nondimen- 
sionally in Fig. 3 as the ratio of water /, to the calcu- 
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Fic. 3. Measured moment of inertia of a liquid in a cylindrical 


cell as a fraction of the moment of inertia of a solid having the 
same density and volume disposition. 


lated /, of a solid having the same density and the same 
volume disposition as the water. 

These model tests are mentioned not for emphasis 
of a particular engineering result, but rather, because 
it is believed that such testing may be a generally useful 
means of determining effective moments of inertia of 
oscillating fuel tanks as functions of shape, frequency, 
and, perhaps, amplitude. Some question logically may 
be raised about scale effects, which are certainly subject 
to investigation. However, from the limited experi- 
ences obtained in the subject programs, it is anticipated 
that useful results may be had from model tests con- 
ducted with oscillation amplitudes sufficiently small to 


prevent sloshing. 


Instrumentation and Test Techniques 


The dynamic testing of the model was conducted in 
two phases: (1) vibration tests, to determine fre- 
quencies and node lines of the lowest natural modes 
of vibration in still air, and (2) wind-tunnel flutter tests, 
to determine flutter frequency, flutter speed, and mode 
damping as a function of air speed. 

For the majority of the tests the wing instrumentation 
was simple. The 
bridges were fed through Brush amplifiers into a two- 
The strain gages were lo- 


outputs of two d.c. strain-gage 


channel Brush recorder. 
cated on the wing spar near the midspan of Wing Sec- 
tion 1. One two-gage bridge had gage placement selec- 
tively sensitive to spar bending; the other bridge, using 
four gages, was made responsive to spar torsion. 

For the vibration tests excitation was provided by a 
motor-driven variable-speed mechanical oscillator at- 
tached to the model tip pod by a small elastic cord. 
Frequencies of coupled modes were obtained from the 
Brush recordings. Node lines were located by the 
sense of touch or by visual observation of relative am- 
plitudes. 

The wind-tunnel test phase involved relatively simple 
In testing each model configuration, the 
During 


procedures. 
tunnel speed was increased in increments. 
speed changes and at stabilized air speeds, small wires 
attached to the tip rib and extending through the tun- 
nel wall were jerked manually to excite a transient re- 
sponse. The character of the damping indicated on the 
Brush tape served as a guide in estimating the proximity 
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TABLE 3 
Geometric Shape Effects on Flutter Characteristics (Wing Free to Roll) 


Measured Pod Parameters peace “ lutter 

: , equency 

Pod Mr Sap lay vy os 

Designation (Ibs. ) (in. Ibs.) (Ib. in.?) (M.P.H.) (Cycles per Min ) 

II 3.441 —4.447 125.4 160 Ben 
V (aft position) 3.477 —4.369 124.5 156 594 
V (forward position) 3.469 — 4.352 124.3 154 591 
III 3.134 —6.255 114.5 143 624 
V (aft position) 3.157 — 6.256 11] : 139 630 
V (forward position ) 3.151 —6.280 112.5 141 628 


Note: A wing weight distribution somewhat different from that shown in Table 1 was used in these tests. 


to flutter and, hence, the magnitude of the succeeding 
speed increment. 

The establishment of 
inclined to be somewhat arbitrary, 
resolution varying with different pod loading conditions. 
This is, of course, to be expected, since flutter is qualita- 
tively described in categories between “mild” and “‘ex- 
plosive,’’ these adjectives relating to the slope of the 
damping-air-speed curve in the region of marginal 
stability. With some of the tested configurations the 
tunnel speed could be varied over a 10-m.p.h. range with 
no appreciable change in flutter amplitude. With a 
few other configurations, a 1-m.p.h. air-speed increment 
produced a definite transition from stable to divergent 
oscillation. Whenever instability occurred, the control 
wires were tightened in an attempt to limit flutter am- 
plitude while the tunnel speed was being reduced. 
vibration tests of only a few of the pre- 


the critical flutter speed was 
the sharpness of 


Generally, 
selected configurations were run prior to tunnel testing, 
inasmuch as the test program was continually subject 
to change. However, after each tunnel test program, 
natural mode data were measured with all model ar- 
rangements for which flutter data had been obtained. 
Thus the repetition of certain of these vibration tests 
served to indicate whether the model had suffered 
structural damage in any of the more violent cases of 
flutter. 

Initial tests of Pods I and III were conducted in the 
5-Ft. Wind Tunnel at Wright Field. All subsequent 
flutter tests were conducted in the Convair Wind Tun- 
nel at San Diego. For the majority of the tests, the 
limiting air speed of interest was arbitrarily placed at 
165 m.p.h. 

Observations made during the earlier tests indicated 
that the lowest flutter speeds were to be expected with 
antisymmetric flutter modes. For this reason most of 
the subsequent testing was conducted with the wing 
free to roll. Limited testing with the wing cantilevered 
was programmed on the general observance that sym- 
metric flutter had occurred with only certain of those 
tested pod configurations that had centers of gravity 
aft of the wing elastic axis. In these few cases flutter 
speeds were only slightly below the preselected maxi 
mum test speed. Consequently, the test results con 
sidered to be of greatest interest are those pertaining to 


antisymmetric flutter modes. 


Test Results 


No attempt is made herein to present a systematic 
review and interpretation of all test data. Rather, 
chosen results illustrating types of flutter behavior are 
offered primarily to emphasize a general need for thor 
ough consideration of the quantity and the disposition 
of the expendable content within a tip pod. 

Several tests were conducted to explore the efficacy 
of pod ballast weights in improving a marginal flutter 
situation. 
tential sensitivity of flutter speed to minute changes in 
pod loading. Arbitrarily, 
tion of Pod IV—one with which the observed flutter 
speed was 138 m.p.h. 
figuration. A weight equal to | per cent of this pod 
weight was located successively at different positions 


These tests give direct evidence of the po- 
particular loading condi 


was taken as the reference con- 


along the axis of the pod, and the flutter speeds corre- 
sponding thereto were noted to range between 134 
m.p.h. and a speed above 165 m.p.h., as shown in Fig. 
t. 

Adherence to the objectives of the test programs pre- 
systematic isolation of pod aerodynamic ef 
fects from the effects of mass distribution. However, 
four special tests were made as follows: Pod V, in 
each of its two positions on the wing, was loaded to re 


cluded 


produce the mass parameters of two arbitrarily chosen 
test configurations of Pod III, and flutter data were ob 
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LOCATION OF “ONE PERCENT WEIGHT”, INCHES AFT WING ELASTIC AXIS 


Fic. 4. Variation of flutter speed with location of a small mov 
able weight in Tip Pod IV. 
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al TABLE 4 
Tip Pod IV Mass Parameters for Conditions Represented in Fig. 6 


Nominal 


Fuel 
Content, 

Test (Per Cent My 
Number Capacity) (Lbs. ) 
9 0 1. 825 
62 11 2.167 
61 22 2.509 
60 33 2.850 
59 4 3.192 
57 56 3.533 
56 67 3.874 


Note: The wing mass parameters are given in Table 1. 


tained for each condition. The test results, given in 
Table 3, show that the changes in pod geometry relative 
to the wing had small effects on flutter characteristics, 
the largest differences shown being only slightly greater 
than limits of test accuracy. While no general opinions 
can be based on these limited data, the results support 
the tentative conclusion of Sewall and Woolston! that 
changes in pod aerodynamic shape have little effect on 
flutter characteristics. 

An unwanted by-product of the design-development 
type of approach pursued in the subject tests has proved 
to be a complication in correlation and interpretation 
of test data. Typically, a change in pod loading was 
accompanied by changes in all three pod mass param- 
eters, masking a clear-cut relationship between cause 
and effect. For example, in interpreting the results of 
a series of tests involving changes only in pod loading, 
it was sought to correlate observed flutter speeds with 
the three independent variables, 1/7, Se;, and Jay. 
Four-dimensional plots proved to be too bewildering to 
assimilate, even with the aid of notes and color coding. 
A three-dimensional plot using a contour representation 
of flutter speed in a planar coordinate system having 
Tap Say and Sa,/My axes has not inspired confidence 
in either its applicability or its utility because of its 
resemblance to a topographical map of extremely rugged 
terrain. Had a research approach been used in these 
investigations, the trends in flutter behavior might have 
been made more apparent by the systematic control of 
one variable ata time. Even so, from the unproductive 
attempts at cross-correlation of data applying to each of 
several pods, it is contended that useful empirical flut- 
ter-speed relationships cannot be expressed in terms of 
pod mass parameters only and that insight into the 
theoretical aspects of the problem is needed to evolve 
suitable correlation bases. 

This contention is supported by several cases of ap- 
parent discontinuity in flutter behavior. A case of 
this type was observed in the testing of Pod IV. For 
the purpose at hand it is sufficient to represent this pod 
as one containing two fuel compartments having a com 
mon bulkhead just aft of the wing elastic axis. 
the fuel sequencing programs of interest involved an 


One of 


empty forward compartment and various amounts of 


Sa, Ta, 
(in.Lbs.) (Lb. In.?) 
—2.353 68.83 
—} 845 69.93 
—1.148 71.40 
—0.189 74.16 

1.054 78.74 

2.570 85.53 

4.352 94.89 


fuel disposed against the forward bulkhead of the aft 
compartment, such as would exist in an unaccelerated 
The pod mass parameters applying to 
Fig. 5 shows the 


vertical dive. 
this sequencing are given in Table 4. 
variation of observed flutter speed with aft-compart- 
ment fuel content. An abrupt reduction from a flutter 
speed above 165 to 99 m.p.h. accompanied a decrease 
from 67 to 56 per cent of tank capacity. Violent flutter 
occurred in the testing of the 56 per cent fuel condition 
because of an inadvertent overshooting of the flutter 


speed by 11 m.p.h. 
170 
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compartment empty. 
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TABLE 5 
Data Applying to Illustrative Test Conditions of Tip Pod IV 


Test Test 
No. 40 No. 62 
Mr (\bs.) 5.019 2.167 
Say (in.Ibs.) —7.764 —1.845 
Tay (1b. in.?) 162.4 69.93 
w,; (cycles per min.) 510 675 
we (cycles per min.) 670 830 
wa (cycles per min.) 615 717 
wa (cycles per min.) 536 765 
vy (m.p.h.) 158 163 
wy (cycles per min.) 568 700 


Note: The wing mass parameters are given in Table 1. 


The damping records of three of the tests in this series 
have been analyzed. The results, plotted against air 
speed, are shown in Fig. 6, wherein the damping coef- 
ficient, g, is the measured logarithmic decrement di- 
vided by x. While some liberties have been taken in 
fairing curves through the test points, there nevertheless 
appears to be a family of curves having similar charac- 
teristics. In particular, the curves applying to the 56 
and 67 per cent configurations show the same type of 
slope reversal in the region of small damping, and, in the 
sense of tendency to flutter, the differences in behavior 
can hardly be called discontinuous. However, since 
flutter stability is not regarded in a qualitative sense, a 
true point of discontinuity exists at that pod weight 
configuration having its damping curve tangent to the 
air-speed axis. 


THE THEORETICAL INVESTIGATION 


A theoretical investigation, complementary to the 
experimental investigation, has been directed toward 
evolving a method of flutter analysis applicable to wing- 
and-pod configurations similar to the F-89 flutter model. 
The criteria of acceptability are considered to be the 
practicability of examining a relatively large number of 
pod weight conditions and, of course, the reliability of 
flutter prediction. 

This investigation is not completed as of this writing. 
Moreover, the results obtained to date are not suf- 
ficiently conclusive to establish the validity of the 
chosen approach. For this reason a formal mathemat- 
ical presentation is deleted from this discussion, and, 
in its place, a qualitative summary of the studies is 
given with emphasis on correlations between results 
and assumptions. 

The chosen method of analysis is conventional, con- 
forming in all major respects with the three-dimensional 
treatment set forth by Smilg and Wasserman.” Only 
the antisymmetric case is examined. The analyses are 
based on uncoupled modes of vibration in bending and 
in torsion. It has been assumed that three properly 
chosen degrees of freedom are sufficient to determine 
the two lowest modes of free vibration, as well as any 
flutter mode of interest. 
taken as follows: 


These degrees of freedom are 
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6: a zero-frequency, rigid-body rolling mode 

h: the first uncoupled antisymmetric bending mode 
a: the first uncoupled antisymmetric torsion mode 


A few exploratory analyses have also been made with 
the added degree of freedom of the second torsiona| 
mode. 

The discussion is presented in two parts. (1) the 
free vibration analyses and (2) the flutter analyses. J. 
lustrations are focused on the analyses of two tested 
configurations of the model wing with Pod IV attached. 
These conditions are selected to represent considerably 
different régimes of pod mass parameters; data per- 
taining thereto are given in Table 5. 


The Free-Vibration Analyses 


Free-vibration analyses are made prior to flutter 
analyses to assure compatibility between uncoupled. 
mode representation and measured coupled-mode data, 
The uncoupled modes are defined in terms of frequen- 
cies and normalized mode shapes, which are established 
as follows: 

Uncoupled-Mode Shapes. 
by a wing rolling displacement giving unit deflection 


The 6 mode is represented 


at the tip-pod centerline. 

The a mode shapes applying to various wing-and-pod 
configurations may be obtained by independent Holzer 
analyses of each. However, in all analyses dealing with 
Pod IV, a reduction in labor has been effected by the 
general adoption of the particular shape calculated for 
the case of a mean moment of inertia of this pod. The 
justification of this approximation is supported by 
Holzer analyses of the cases with extreme pod inertias; 
the normalized a mode shapes do not change appreci- 
ably as Jay is varied from a mean value to either ex 
treme. 

The 4 and a modes are orthogonal by definition. The 
h mode is defined so as to be orthogonal with the # 
mode. This definition is effected by approximations 
that permit relatively simple accounting for the change 
in tip-pod mass. Specifically, for a particular wing 
mass distribution containing a preselected reference tip 
mass, the normalized mode shape of the first antisym- 
metric bending mode is obtained as a reference mode 
shape, h,, by any of the conventional iteration proce- 
dures using mass and rigidity distributions. Thereafter, 
with any other tip mass, the normalized mode shape 1s 
approximated by a linear combination of the /, and 
@ modes in a manner to satisfy the orthogonality condi- 
tion. 
gives good agreement between approximate and exact 


As evidenced by special studies, this technique 


mode shapes, particularly if the reference tip mass 1s 
about midway between the extremes to be considered 
with a given tip pod. 
Fig. 7 shows the uncoupled-mode shapes used in 
analyses of Pod IV Test Conditions 40 and 62. 
Uncoupled-Mode Frequencies.—Uncoupled-mode fre- 
quencies are obtained from a free-vibration analysis 
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using measured coupled frequencies and model mass 
parameters. This analysis gives a frequency equation 
containing w,”, @a’, and terms defining the inertia 
coupling between the h and a modes and the a and @ 
modes. This equation (a quadratic in w*) is solved for, 
say, wa? in terms of these same inertia coupling param- 
eters and the roots, w;,and we, measured in vibration 
tests. The solution yields two roots for wa*. The 
ambiguity of which of these roots applies is eliminated 
statistically in the following manner: For each of a 
number of the tested configurations, both roots are 
used in calculating the effective torsional stiffness of 
the wing in the a mode. Examination of all of these 
stiffness values shows that half, represented by one 
from each test condition, are closely grouped, whereas 
the other half differ widely. The effective torsional 
stiffness is also obtainable from wing GJ and a distri- 
butions; therefore, if the use of a single a representa- 
tion is accepted, this stiffness may be considered as a 
physical constant of the wing. Thus it is apparent 
that the small spread of the better grouped set of stiff- 
ness values obtained from the coupled-mode analyses is 
attributable to errors in measured coupled frequencies 
and mass parameters, to the adoption of an average a 
mode shape in place of the exact shape, and to neglected 
coupling with higher vibration modes. The cumula- 
tive errors from these sources are considered to be 
smaller than those incident to static measurements of 
torsional rigidity, and, hence, the coupled-mode ap- 
proach to the determination of uncoupled frequencies is 
preferred. This preference follows from the belief 
that, in the subject application, at least, errors in fre- 
quency are more likely to have adverse effects on the 
flutter analyses than are errors in mode shape when both 
frequency and mode shape are determined by, say, a 
Holzer analysis using a measured GJ distribution. 

The technique of stiffness calculation by coupled- 
mode analysis not only provides a basis for determining 
which of the roots of wa’ applies to any one test configur- 
ation but also affords a means of spotting erroneous or 
inconsistent data applying to any configuration. In 
practice, an average torsional stiffness has been used 
in the determination of the wa value for each test condi- 
tion. The uncoupled bending frequency, w,, is calcu- 
lable as soon as wa is established. 

This procedure has been followed with 30 of the tested 
conditions of Pod IV. The computing labor entailed 
has been kept relatively small by suitable tabular setup. 
Average stiffnesses in bending, as well as in torsion, 
have been used in the calculation of uncoupled fre- 
quencies. 

Prior to flutter analysis, the adequacy of the treat- 
ment thus far is judged by the correspondence of calcu- 
lated and measured node line locations in each of the 
coupled modes. Fig. 8 and Fig. 9 show the node-line 
plots for Test Conditions 40 and 62, respectively. In 
both cases the agreement is considered sufficiently close 
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POD IV- TEST NO. 40 
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Measured and calculated node-line locations, first 
coupled vibration modes of Test Condition 40. 
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Measured and calculated node-line locations, first 
coupled vibration modes of Test Condition 62. 
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to warrant confidence in the assumptions and the tech- 
niques of the free-vibration analyses. 


The Flutter Analyses 


The formulation of modal aerodynamic energies ac- 
counts for the effects of wing taper with the method 
recommended by Smilg and Wasserman. That is, a 
reference semichord is chosen, and the aerodynamic 
coefficients applying to tapered airfoils are used in the 
appropriate aerodynamic integrals. The contributions 
of wing and tip pod are evaluated independently 
and are subsequently combined in the flutter determi- 
nant. 

The aerodynamic effects of an oscillating tip pod are 
not accounted for in available flutter theory. There- 
fore, as an expedient to permit the application of exist- 
ing theory, it is sought to replace the tip pod with its 
flat-plate ‘‘aerodynamic equivalent.” The term “‘aero- 
dynamic equivalent” is defined herein as a rectangular 
wing extension having a semichord, a span extent, and 
a quarter-chord location relative to the wing elastic 
axis which give the same aerodynamic effect under 
analysis by the method of reference (2) as does the ac- 
tual pod under wind-tunnel test conditions. It must 
be granted that pod effects may not be confined within 
the pod span but may extend over the outboard portion 
of the wing; a significant end-plate effect may exist 
in the oscillating case, as well as in the stationary case. 
At the same time, it is considered that the presence of 
the tip pod may strengthen the validity of the treatment 
accorded the wing, inasmuch as the method of reference 
(2) makes no rational provision for finite aspect-ratio 
effects. 

The aerodynamic coefficients contained in reference 
(2) are given as functions of v/bw. Since all parts of 
the wing are operating at the same air speed and flut- 
ter-mode frequency, the coefficients applying to the 
aerodynamically equivalent tip pod are made compat- 
ible with those applying to the wing by use of v/bw 
values consistent with the relation between pod semi- 
chord and wing reference semichord. It follows that 
the relationships between the moduli and the phase 
angles of the aerodynamic vectors associated with the 
pod and with the wing are sensitive to the choice of pod 
semichord. 

One method of determining the parameters of the 
tip-pod aerodynamic equivalent is the method of direct 
substitution. That is, the semichord, the span, and 
the quarter-chord location are arbitrarily chosen, and 
flutter analyses are made in the conventional manner. 
A failure of this tip representation to give acceptable 
determination of flutter characteristics demands the 
choice of a new set of parameters and a repetition of the 
analyses. In principle, this process is repeated until a 
unique pod configuration is found to give reasonable 
correspondence between calculated and observed flutter 
behavior in each of the examined test cases. 
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TABLE 6 
Assumed Flat-Plate Aerodynamic Equivalents of Tip Pod ]\ 
(Pod Diameter at Maximum Section: 4.73 In.) 


Quarter-Chord 
Location 
br{( 1/2) + a} 


Semichord Span (In. Forward 
Case br (In.) sr (In.) Elastic Axis 
1 5.5 4.73 0.88 
2 16.5 r ie if 6.33 
3 16.5 3.60 0) 
} 5.5 4.73 6.33 
5 8.7 3.00 6.27 
6 5.5 4.73 3.63 


This approach was adopted in preliminary flutter 
analysis directed at establishing the sensitivity of cal- 
culated flutter characteristics to the choice of widely 
different pod plan-form arrangements. Successively. 
six different aerodynamic equivalents were assumed to 
apply to Pod IV. 


given in Table 6. 


The geometries of these tip cases are 
Case | is merely an extension of the 
wing tip chord over the span of the pod. Case 2 js 
synthesized from estimates of the force and moment 
coefficients applying in the nonoscillatory condition 
Case 3 is a modification of Case 2 in the sense that a 
fictitious horizontal fin is added to move the center of 
pressure aft to the wing elastic axis; however, the ef 
fective chord is unchanged, and the added lift is ac 
counted for by an increase of pod span. Case 4 utilizes 
the plan form of Case | and the quarter-chord location 
in Case 2. Case 5 is essentially a modification of Case 
!, having the same plan-form area and nominally the 
same quarter-chord location but an increased semi 
chord. Case 6 reverts to the plan form of Cases | 
and + but adopts an intermediate quarter-chord posi- 
tion. 

None of the six arbitrarily chosen tip-pod equivalents 
listed in Table 6 has given consistent agreement be 
tween calculated and observed flutter characteristics 
This failure has prompted a re-evaluation of the choice 
of degrees of freedom to the end that a few supplemen 
tary analyses have been made with the addition of the 
second uncoupled torsion mode. 

Some results of the flutter analyses are illustrated for 
Test Conditions 40 and 62. Only those calculations 
pertaining to the flutter mode of interest are examined 
Fig. 10, applying to Test Condition 40, is a conven 
tional plot of damping coefficient against air speed, 
where, in contrast to the connotation used in Fig. 6, 
the damping is that required to sustain simple harmonic 
motion in the flutter mode. The critical flutter speed 
is that at which the required damping is equal to the 
structural damping, taken, typically, as 0.02 on the 
Case 1, with either the 


three or the four degree of freedom analysis, gives ac- 


basis of model measurements. 


ceptable agreement with the observed flutter speed. 


It may be conceded that the better agreement pro- 


vided by the three degree of freedom analysis is virtu- 


ally coincidental in view of the indicated sensitivity 
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FLUTTER 


of critical flutter speed to the assumption of pod effec- 
tiveness. Of equally important note, both analyses 
using Case 1 give flutter frequencies within 2 per cent 
of the observed flutter frequency, whereas the Cases 2 
and 4 show a progressive reduction in the flutter-mode 
frequency with increasing air speed, pointing toward a 
substantially low flutter frequency at an excessively 
high flutter speed. 

Fig. 11 shows damping-air speed plots calculated for 
each of five tip cases applied to Test Condition 62. 
Here, Case | fails to give agreement with observation, 
whereas Cases 2 and 5 give excellent correspondence 
with observed flutter speed. However, the flutter fre- 
quencies in these latter two cases are calculated as 320 
and 400 cycles per min., respectively, as compared with 
a measured flutter frequency of 700 cycles per min. 
The flutter speed computed with Case 6 is only about 
10 per cent high, but the corresponding flutter frequency 
is 450 cycles per min., reflecting only a slight improve- 
ment over Case 5. Thus none of the five aerodynamic 
tip cases tried with Test Condition 62 gives acceptable 
agreement between calculated and observed flutter 
characteristics. 

In general, the results of the analyses of these and a 
few other test configurations indicate both that the tip- 
pod representation is an extremely important factor in 
the flutter analyses and that the number of degrees of 
freedom chosen is of considerably lesser importance. 
Further, it appears that (1) the aerodynamic equivalent 
of a pod is small as compared with the pod plan form 
and (2) the importance of accurate pod aerodynamic 
representation is greatest for those configurations hav- 
ing relatively small pod mass unbalance about the wing 
elastic axis. 

In view of these indications, effort is currently being 
directed toward evolving a more positive means of de- 
termining the aerodynamically equivalent tip pod. 
The approach is as follows. 

A “sensitive”’ test configuration— that is, one having 
a small pod-unbalance—is subjected to an analysis 
using the wing reference v bw agreeing with measured 
flutter speed and flutter frequency. The flutter de- 
terminant for the three degrees of freedom, 6, h, and a, 
is set up to retain the pod aerodynamic terms as de- 
pendent variables. The measured flutter frequency 
and the model structural-damping coefficient are used 
in forming a known characteristic complex root of this 
determinant. This root is then substituted in the de- 
terminant to leave the pod aerodynamic terms as the 
only unknowns. Arbitrary assignment of a trial value 
of equivalent-pod semichord permits evaluation of these 
aerodynamic terms as complex functions of span, s7, 
and dimensionless quarter-chord location, [(1/2) + a]. 
Subsequent manipulation leads to a third-order deter- 
minant of rank two having characteristic roots, 1/sr. 
Mathematically, these roots are determined by the par- 
ticular values of 1/sr and [(1/2) + a] that equate the 


determinant to zero. Physically, the roots are deter- 
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to Test Condition 60 with several of the tip cases of Table 6. 


mined by the condition that they be positive real num- 
bers of reasonable magnitude. 

With a preselected pod semichord, 67, the solution for 
A trial 
A series of physically 


an applicable root, 1/s7, proceeds as follows. 
value of [(1/2) + a] is assigned. 
reasonable | /s7 values is used successively in evaluating 
the determinant as a complex number. Moduli of each 
determinantal value are plotted against 1/s7, and the 

This process is repeated 
+ a], and the minimum 


minimum modulus is sought. 
with other values of [(1/2) 
moduli so obtained are plotted against [(1/2) + a]. 
If a real solution exists, the minimum modulus of this 
plot will be zero or a value differing from zero only be- 
cause of significant-figure limitations in the computa- 
tion stages. If a real solution is not obtained with the 
preselected pod semichord, then a new value of 67 is as- 
signed and the process is repeated. The values of sr 
and [(1/2) + a] which produce a zero modulus of the 
determinant complete the delineation of the aerodynam- 
ically equivalent tip pod applying with the chosen 
semichord. 

It is to be expected that a range of br values, each as- 
sociated with particular values of sr and [(1/2) + a], 
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can be found by this method to apply to any one test 
configuration. The particular set of these parameters 
that represents the tip pod is, in theory, determined by 
its satisfying the flutter-analysis requirements with 
other tested configurations. Further, if the assump- 
tions of these flutter analyses are valid, the character 
of a damping-air-speed curve obtained by use of other 
v/bw values should be consistent with the complemen- 
tary damping relationship obtained by test. 

This procedure is admittedly a brute-force method of 
ascertaining the existence of a tip-pod aerodynamic 
equivalent conforming with the definitions and hypoth- 
eses of the theoretical treatment. However, the 
labor entailed in the evaluation of the modified deter- 
minants has been made practical by the adoption of 
IBM computing facilities. 

The present status of the investigation does not per- 
mit a prediction of the overall reliability and practica- 
bility of the method. It is to be hoped, however, that 
the calculations now in progress will lead to a positive 
evaluation of the validity of the chosen approach. 

Some preliminary results of interest have been ob- 
tained by analyzing Test Condition 62 for its equivalent 
tip-pod parameters, sy and !(1/2) + a], as functions of 
measured flutter data and a selected pod semichord of 
5.5 in. This value of br was chosen merely because it 
was a Case | parameter contributing to a prior success 
in calculating flutter frequency and critical flutter speed 
of Test Condition 40. The analysis shows that the ob- 
served flutter conditions are satisfied with a pod semi- 
chord of 5.5 in., a pod-span of 1.5 in., and a quarter- 
chord location 5.0 in. forward of the wing elastic axis. 

This aerodynamically equivalent plan form, called 
Case 7, has been subsequently applied in a reanalysis of 
Test Condition 40. The results show a calculated 
flutter speed of 132 m.p.h. as compared with an ob- 
served flutter speed of 158 m.p.h., both values applying 
with a structural-damping coefficient of 0.02. The 
calculated flutter frequency is within | per cent of the 
measured flutter frequency. 

Results of a few other flutter analyses using the Case 
7 pod plan form are as follows: 

The three test conditions (60, 57, and 56), represented 
in Table 4 and Fig. 6, have been analyzed. For Test 
Condition 60, the calculated flutter speed of 109 m.p.h. 
compares favorably with the observed speed, 118 
m.p.h.; corresponding flutter 
and 700 cycles per min. However, the analyses of 
Test Condition 56 and 57 have failed to show agreement 
with the observed flutter behavior. Yet, for the latter 
of these conditions, the shape of the calculated rela- 
t.onship between damping coefficient and air speed is 
consistent with the measured relationship of Fig. 6 in 


frequencies are 699 


the respect that a point of minimum stability is shown 
at 108 m.p.h.; the damping coefficient is computed to 
be —0.001 at this speed, or incrementally, about 0.021 
more negative than that indicated by test. 
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Test Condition 37 has pod mass parameters, i, = 
3.995 Ibs., Sap = —11.48 in.Ibs., Jay = 148.3 Ib, in? 
The calculated flutter speed obtained with the Case 7 
plan form is 11 m.p.h. below the observed flutter 
speed of 157 m.p.h.;_ the calculated flutter frequency js 


3 per cent higher than the measured frequency, 639 


cycles per min. 

Two of the test conditions represented in Fig. 4 also 
have been subjected to flutter analysis. The condition 
having the incremental weight located 8.3 in. aft of the 
wing elastic axis has a calculated flutter speed 15 m.p.h, 
below the observed speed of 138 m.p.h.; agreement be- 
tween computed and measured flutter frequencies js 
good. Calculations for the condition having the incre- 
mental weight 17.1 in. aft of the elastic axis show a flut- 
ter speed in excess of 165 m.p.h. but a point of minimum 
stability (g = +0.016) at 130 m.p.h. 
the damping records obtained during the wind-tunnel 


Significantly, 


test of this condition show behavior of a similar charac- 
ter, with minimum damping occurring somewhere in the 
range between 135 and 145 m.p.h. 

The Case 7 tip-pod equivalent appears to be far supe- 
rior to any of the wholly arbitrary Cases 1 through 6 be 
cause of its leading to reasonably good agreement be- 
tween calculated and observed flutter characteristics 
in six of the eight analyses made to date. Moreover, it 
is considered unlikely that this plan form is the best 
representation that may be made, since the semichord 
was selected more or less arbitrarily. Therefore, it may 
be expected that an extension of the treatment could 
lead to a somewhat different set of parameters, 57, 57, 
and [(1/2) + a] giving an overall improvement in 
analytical results. 


The limited test results presented in Table 3 are not 
inconsistent with the results of the flutter analyses de- 
scribed herein. For these test cases of moderately 
large pod overbalance, the inertia coupling among 
modes apparently predominates over the changes in 
aerodynamic coupling to reflect a relative insensitivity 
of flutter characteristics to the imposed changes in pod 
shape or pod location on the wing. Had this type of be- 
havior been sensed at the time of these tests, additional 
tests of the same nature would have been conducted 
with pod configurations having centers of gravity close 
to the wing elastic axis. It is surmised that such tests 
might reveal a greater sensitivity of flutter speed to the 
Sewall and Woolston have 
made limited comparison tests of this type on a straight 
Their data, covering two cases of 


changes in pod geometry. 


cantilever wing. 
large and small unbalance of the concentrated tip 
weight, indicate that the flutter speeds and frequencies 
are relatively unaffected by radical changes in the aero- 
However, they do 
Indeed, 


dynamic shape of the tip weight. 
not conclude that such is generally the case. 
it is contended herein that entirely different degrees of 
importance of pod aerodynamic effects may be antici- 
pated in different régimes of pod weight distribution 
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and, for that matter, in symmetric and antisymmetric 


flutter cases. 
CONCLUSION 


The potential sensitivity of wing-and-pod flutter 
characteristics to relatively small changes in pod mass 
parameters has been demonstrated by test. 


The theoretical studies indicate that, while the aero- 
dynamically effective tip pod may be small in compari- 
son with pod geometric plan form, an accurate represen- 
tation of this effective tip pod may be essential to a suc- 
cessful flutter analysis. Further, it is indicated that 
analyses of wing-and-pod configurations having the pod 
center of gravity on or near the wing elastic axis are 
more adversely affected by errors in pod aerodynamic 
representation than are those analyses applying to cases 


of large pod overbalance. 


Conclusions on the general reliability and the practi- 
cability of the investigated methods of flutter analysis 
are deemed premature in consideration of the fluid sta- 


tus of the development. From a positive point of view, 
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there is some basis for expecting an improvement in 
analytical predictions of flutter speeds, although the 
embryonic method of analysis requires the use of appli- 
cable test data. 

Accordingly, the view is currently held that flutter- 
model testing offers the more certain and the more di- 
rect method of evaluating the flutter safety of wing-and- 
tip pod configurations prior to flight test of the aircraft. 
Adoption of this view forces attention on the problems 
incident to attainment of dynamic similarity between 
model and prototype not only with respect to simula- 
tion of elastic behavior but also with respect to simula- 
tion or conservative representation of fuel inertia prop- 
erties. 
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The Effect of Nonlinear Aerodynamic 
Characteristics on the Dynamic Response 
to a Sudden Change in Angle of Attack 


TELFORD W. OSWALD* 
Douglas Aircraft Company, Inc. 


ABSTRACT 


The effect of nonlinear pitching moment and lift variations 
with angle of attack on the dynamic response to a sudden change 
in angle of attack is considered. An approximate solution to 
the nonlinear equation of motion is developed. Several numeri 
cal examples are considered, and the results of the approximate 
solution are compared with the accurate results of numerical 
integration, as well as the classical linearized solution. The 
effect of a nonlinear moment curve on the determination of sta 
bility derivatives from flight-test data is discussed in the light 


of these examples. 


(I) INrRODUCTION 


bes SOLUTIONS TO THE EQUATIONS OF MOTION of an 
airplane when subjected to small disturbances 
from steady flight have been well known to aeronauti 
cal engineers for the past several years. The classical 
treatment of the problem is due largely to G. H. Bryan, 
who applied the theory of small oscillations previously 
developed by Routh. Bryan assumed that the forces 
and moments due to a slight disturbance from a state 
of equilibrium depend linearly on the disturbance. The 
solution is then shown to depend on a number of con- 
stants called stability Until recently, 
there has been little or no change in the theory, and it 
is in use today essentially in the form developed by 


derivatives. 


Bryan. 

It should be noted that the classical theory applies 
strictly only in the case of infinitesimal disturbances. 
However, the results have been found to apply with 
reasonable accuracy to finite disturbances of the mag- 
nitude experienced in the flight of conventional sub- 
sonic airplanes. In other words, the second-order 
terms neglected in the theory do not have any appre- 
ciable effect on the motion. 

However, in the case of high-speed airplanes and 
missiles, serious deviations from the conditions of the 
classical theory occur when one attempts to extend the 
results to the case of finite disturbances. In most 
nissile configurations, small aspect ratio wings are 
mounted in tandem on proportionately large diameter 
bodies. This exaggerates the aerodynamic interfer- 
ence effects, and such missiles exhibit markedly non- 
linear dependence of the aerodynamic forces and mo- 
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ments on the angle of attack. Consequently, the am. 
plitude of disturbances for which one would expect th, 
classical theory to be valid is considerably reduced 
At the same time, in the normal operation of such 
missiles, large sudden changes in the angle of attack ar 
frequently required. Hence, it would appear that, j: 
the case of supersonic missiles, the extension of th, 
theory of small oscillations to finite disturbances js 
questionable. 

In the present paper, the effect of nonlinear variation 
of lift and pitching moment with angle of attack on the 
dynamic response to a step function input in angle of 
attack is investigated. The missile is assumed to be 
flying at constant speed and altitude. 
reduced to an ordinary second-order nonlinear differ- 


The problem is 


ential equation that cannot be solved exactly by meth 
ods known today. However, approximate solutions 
may be obtained, and, by the method of numerical 
integration, a solution to any desired degree of accuracy 
can be obtained. 

The stability of the nonlinear system is determined 
from considerations of the solution in the neighborhood 
of the critical points of the differential equation, and the 
investigation of the dynamic response is confined to 
systems that have been previously determined to be 
stable. 

In its application to a guided missile the problem of 
dynamic response properly involves the combination of 
the missile as an aerodynamic body and the control 
system. But since ultimately any motion of the mis 
sile must be brought about by aerodynamic forces, it 1s 
desirable to know the dynamic behavior of the missile 
alone. In particular, the ‘‘overshoot’’ characteristic 
of the response to a sudden change in trim position 1s 
important in the design of the control system, since 1l 
indicates the amount of “‘feedback’’ necessary to pro 
vide satisfactory steering of the missile. 

In the present analysis, an attempt is made to ap- 
proximate this overshoot without resorting to the 
tedious method of numerical integration. A compari 
son of this approximation with the numerical solution 
and the linearized solution is presented. 

An important problem confronting the aerodynamt- 
cist working in the field of guided missiles is the reduc- 
tion of flight-test data in order to evaluate the sta- 
bility derivatives. One of the procedures followed is to 
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analyze the transient response of the missile to a step- 
function input in the control-surface deflection. This 
is equivalent to a step-function input in the angle of 
attack. In Appendix (II) of this paper the effect 
of a nonlinear moment curve on this procedure is 


discussed. 


I] COORDINATES, NOTATION, AND SYMBOLS 


The equations of motion for a missile can be referred to a set 
of body axes that are fixed in the missile and move with it. An 
orthogonal set of principal body axes with origin in the missile 
center of gravity (the x-axis the longitudinal centerline, the y-axis 
in the normally horizontal plane, and the z-axis in the normally 
vertical plane) is selected. This coordinate system is shown in 
Fig. 1. The 
ments, and angular displacements are all defined in accordance 


velocities, moments, angular velocities, displace 


with the right-hand rule 

Throughout the text, the usual convention of placing a dot 
over a quantity to denote differentiation with respect to time is 
used. Whenever a fractional power of a quantity is involved, 
the positive real branch of the multivalued function is taken 
For the inverse trigonometric functions, the principal branch is 
taken. In the following list the principal quantities used in this 


paper are defined: 


1, B, C, D, E, F = coefficients of nonlinear equation of mo 
tion —defined in text 

a = speed of sound, ft. per sec. 

( = reference chord for moment coefficients, 
ft 

Cr = 1/(1/2)pV%s = lift coefficient 

Cy = M/(1/2)pV%sc = pitching moment co- 
efficient 

CL. = 0C//da, evaluated at trim position 

Cu. = O0Cys/Oda, evaluated at trim position 

Cu, = 0Cy/O a 27 ), evaluated at a = 0 

Cu, = 0Cy/0(qc/27), evaluated at gq = 0 

iy, d = incomplete elliptic integral of first kind 

g = acceleration of gravity, ft./sec.? 

I. = S (9? + s?)dm = moment of inertia 


about x-axis (roll), Ibs. sec.?/ft. 
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moment of inertia 


S (x2 + s*)dm = 


about y-axis (pitch), Ibs. sec.?/ft 


S (x? + 


about z-axis (yaw), lbs. sec.?/ft. 


y?)dm = moment of inertia 


constants used to define nonlinear lift 
curve 
constants used to define nonlinear mo 


ment curve 
complete elliptic integral of first kind 
lift, Ibs. ; 


mass, lbs. sec.?/ft 


rolling moment, Ib. ft. 
pitching moment, lb. ft.; Mach Numbei 
accelerometer reading 


yawing moment, lb. ft 


¢ = O¢/Of = perturbation in roll rate 
(angular velocity about x-axis), rad 
per sec 

&@ = 0/df = roll rate, rad per sec 

steady-state value of P 

6 = 06/df = perturbation in pitch rate 
(angular velocity about y-axis), rad. 
per sec 


0 = 00/Of = pitch rate, rad. per sec. 
steady-state value of Q 


(1/2)p 72 = dynamic pressure, Ibs. per 
sq.ft. 

y = Oy/df perturbation in yaw rate 
(angular velocity about axis), rad. 
per sec 

W = OW/df = yawrate, rad. per sec 

steady-state value of R, rad. per sec 

radius of curvature of flight path, ft 
reference area for aerodynamic coeffi- 


cients, sq.ft. 

time, sec 

velocity perturbation in x-direction, ft 
per sec 

velocity component in x-direction, ft 
per sec, 

steady-state value of lL’, ft. per sec 

velocity perturbation in y-direction, ft 
per sec 

velocity component in y-direction, ft. per 
sec 

steady-state value of V, ft. per sec 

missile velocity, ft. per sec 

velocity perturbation in z-direction, ft 
per sec 

velocity component in z-direction, ft. per 
sec 

steady-state value of W’ 

coordinate’ axis, 

body axes (roll, pitch, yaw axis, re- 


right-handed principal 
spectively ) 

force in x-direction, Ibs 

force in y-direction, Ibs 

force in s-direction, lbs 

OX/Oun = 
change of force in x-direction with ve- 


resistance derivative, rate of 


locity in x-direction, lbs. sec. /ft. 


)p, v, 0, r = resistance derivatives, Ibs 
sec. /ft 
N \p,v,v,r = Totary derivatives, Ibs. sec 
(M lw, w,q ly 


angle of attack, rad 
damping constant, 1 /time 


w/U, = incremental angle of attack, rad. 


perturbation angular displacement in 


pitch, rad. 
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0 = angular displacement in pitch measured 
from horizontal, rad. 

0; = steady-state value of 0 

p = mass density of air, lbs. sec.?/ft. 

T = time, sec. 

w = circular frequency, rad. per sec. 

Wn = natural undamped circular frequency, 
rad. per sec. 

Q = phase angle, rad. 


(III) Eguations or MOTION 


The equations of motion are referred to a set of prin- 
cipal body axes; hence, all product of inertia terms 
vanish. It is assumed that the missile is initially in 
a state of steady flight, that the linear and angu- 
lar velocities and displacements can be represented 


as 
U=Uritu P=Pitp, 0 =%+80) 
V=aNt+e, O=-At+a ®=&4+¢e) (1) 
W=W4+4, R=k,+,, v=v,4+y) 


and that P; = R,; = 2 = VY = 0. It is also assumed 
that 6 is sufficiently small that sin @ may be taken equal 
to @ and cos @ may be taken equal to | and that u, 
v, w, p, g, r are sufficiently small that their squares 
and products may be neglected. The equations of 


motion are, then, 


m(i— Vir+ Wi(Qit+gq)]) = X — mg sin 0, — 
mg cos 0,6 
mo — Wip + Ur) Y + mg cos Oi¢ (2) 
m(w — U\(Qi + g) + Vip] = Z + mg cos 0; 
LP= L; 1q=M; lr =N 


If the motion is confined to infinitesimals, each of the 
air reactions can be expressed in the form 


tn hs w+ wes a+ 
—— ow ‘ av‘ 
ar aX aX 
10 iP IR 
a ae * on 


where the second-order terms are neglected. The par- 
tial derivatives are then written in the form 0X /OU’ = 
X, and are known as “‘stability derivatives.’ To 
investigate the effect of the nonlinear variation of 
pitching moment and lift with angle of attack on the 
dynamic behavior of a missile, the second- and third- 
order terms in dW are included in the expressions for 
M and Z. Accordingly, instead of the usual notation 
M,,w, Z,w, these terms are written M(w), Z(w), respec- 
tively. 

If the missile has a longitudinal plane of symmetry 
(which is usually the case) and if it is flying at zero 
yaw, many of the stability derivatives are zero, and the 
set of six simultaneous equations [Eqs. (2)] reduce to 
two sets of three equations each, 


m(u + Wig + g cos 0,0) = X,u + X,w + 
Xq + X1 — mg sin 0; — mW1Q; 


(3.1) 
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m(w— Uig+gsin00) = Z,u + Z(w) 4 
Z.9 + Zyw + Z; + mg cos 0; + mUiOQ, (32 


Lg = Mu-+ Mw) + Mg++ Myw+ M, (3: 


mo + Uy — WiP — gcos Qg)= Vv + Y,p + 
Y,r am Yj = Y; (4, 


I,p = Lw = Lyp + Ly + Ly (4.9 
Li = Nwvu+NepP+Nr+NG@+tM, (43 


It is to be emphasized that these equations are re 
ferred to principal body axes instead of the more com 
monly used ‘“‘wind’’ axes or so-called ‘‘stability” axes 
and that the appropriate expressions for the stability 
derivatives must be used. 

Considering the longitudinal Eqs. (3), it is observed 
that, if Z, = M, = 0, Eqs. (8.2) and (3.3) become 
independent of Eq. (3.1), and the pitch motion can be 
described by a two degree of freedom system. For a 
missile that is symmetrical with respect to the yaw 
plane and which is trimmed at zero angle of attack, it 
is apparent that the lift and moment on the missile are 
zero and that their derivatives with respect to u will be 
identically equal to zero. Hence, it is to be expected 
that at small angles of attack the coupling of u with 
w and @ is small and that the system can be adequately 
treated as a two degree of freedom system. This means 
that the influence of the long-period motion on the 
short-period motion is negligible. Calculations for 
representative missiles have been carried out, and this 
is indeed found to be the case. The results of these 
calculations are presented in Table |. 


TABLE | 
Effects of Long-Period Oscillation and Gravity on the Charac- 
teristics of the Short-Period Oscillation for a Representative 


Missile 

Mach Number 2.45 1.42 
Altitude (ft.) 40,000 70,000 
Angle of climb 83° 15° 
Three degrees of freedom including gravity 

Damping constant 0.76 0.12 
Circular frequency 4.78 1.38 
Two degrees of freedom including gravity 

Damping constant 0.76 0.12 
Circular frequency Wes 1.38 
Two degrees of freedom neglecting gravity 

Damping constant 0.75 0.12 
Circular frequency 4.78 1.38 


It is noted that the long period for a supersonic 
missile is of the order of 100 sec., whereas the short 
period is of the order of 1 sec. Therefore it is to be ex- 
pected that the long period should have no appreciable 
effect on the short period. Furthermore, it is obvious 
that the short-period motion is the one that is impor- 
tant for steering and control. The long-period motion 
would appear to be of little practical importance, since 
the assumptions made will certainly be violated during 
the course of the long period. 
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For equilibrium: 
F=-Z,cos a + X,sina-mg cos(@-a) 
X,cos a=mg sin (@-a) 


For small values of a , neglecting squares 


F =-Z,-mg cos®, 


Fic. 2. Forces acting on missile turning in vertical plane. 


If it is assumed that the missile is stabilized in roll, 
the lateral Eqs. (4) reduce to exactly the same form as 
the longitudinal equations except for the gravity terms. 

Returning to the longitudinal equations, dividing 
Eq. (3.2) by m and Eq. (3.3) by J, and designating by 
primes the quantities so divided, the equations become 


(1 — Z,’)w — Z'(w) — (Ui + Z,')q = Z1' — 
4 sin 0,0 + g COs 6, + UiQi 


M’'(w) — M,’qa+q@ = M)’ 


(5.1) 


—M,/w — (5.2) 


There remains a bothersome term in these equations 

namely, g sin 0,0. If the missile is flying at a small 
angle of climb, sin 0, = 0 and this term may be neg- 
lected. However, as is often the case, if the missile is 
flying in a nearly vertical direction, it is not apparent 
that this term may properly be neglected. Further- 
more, its presence gives rise to an additional root to the 
characteristic equation, and this additional root is 
divergent. The divergence is slow, however, and it 
can be shown that it represents the curvature of the 
flight path due to gravity. This curvature will have a 
negligible effect on the oscillations about the flight path 
during the time interval for which the motion is being 
Computations for representative cases have 
These 


studied. 
been made, and the results are shown in Table 1. 
results indicate that there is no appreciable error in- 
volved in neglecting gravity. 

Noting that Z,’ < U, and Z,’ < 1, hence Z,’ and 
Z,' can be neglected and Eqs. (5) are written 


w— Z'(w) — Uig = Z)' + UiQi + g cos Oy 
—M,'’w — M'(w) — M,'q+q¢= M,’ 


(6.1) 
(6.2) 


Solving Eq. (6.1) for g and substituting into Eq. (6.2), 
the system is reduced to the one equation of motion, 
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w Z'(w) : w  Z'(w) 

=) _ '(w) — M’q + _- — ) _ 
l 1 l 1 l 1 l 1 


/ 


M 
M,,'w = M,’ = 7, (Z,’ + UiQ) + g COs 01) ( 
1 


“J 


) 


If it is desired to compute the response to a step 
function input in control-surface deflection, the forcing 
functions, 4,’ and Z,’ can be written as step functions 
and the response computed directly from Eq. (7). The 
incremental velocity, w, is measured from the original 
trim position, and Q, and ©, correspond to the steady- 
state values at the original trim position. Alterna- 
tively, a linear transformation in w can. be made such 
that w = 0 corresponds to the trim position at the new 
control-surface deflection, whence Q; and 0, are referred 
to the new trim position. The solution of the trans- 
formed equation is found subject to the initial condi- 
tions, w(0) = wy) and w(0) = 0, where w» corresponds 
to the old trim position measured in the new coordinate 
system. If this latter point of view is adopted, M,’ = 
0, Z,’ is the aerodynamic force in the z-direction, and 
Q, is the angular velocity for the steady-state motion 
in the new trim position. For the steady-state motion 
of the missile in a vertical plane, equilibrium of forces 
in the z-direction, neglecting squares of the angle of 
attack (Fig. 2), requires 

—Z, — mg cos O; = mMRQ? 


where R, is the radius of curvature of the flight path. 
But 7; = Q,Ro, hence 

—Z, — mg cos 0, = mU\Q; or —Z’ — gcos®, = UiQ; 
and the right-hand side of Eq. (7) is zero. The problem 
is now reduced to finding a solution to the homogeneous 


equation, 
® 7 2"@) _ Mw) — M,! (? a 2) 7” 
U7 l U; U; 
M,'w =0 (8) 
with the initial conditions 
w(0) = wm, w(0) = 0 


In order to formulate the problem in a manner amenable 
to solution by analytical means, the pitching moment 
and normal force as functions of the angle of attack are 
approximated by polynomials. It is convenient 


to use cubic polynomials of the form 

M(d) = M9 + Kw? + Kod? 

Z(v) = 24) + kw? + kod? 
where # is the increment in angle of attack measured 
from the trim position. Since w is a small quantity 
compared with U, 8 = w/U, neglecting higher order 
terms. Substituting these relations into Eq. (8) and 


collecting terms, 


d3+(A + Bd 4+ Ed*)d + (C+ Dd + Fd’) 3 = 0 
(9) 








306 JOURNAL OF THE 
and the initial conditions are 
H(0) = ao = w%/Ui, 90) = 0 
where 
La E q' Se? 
bie Eh ee oe Oe a 
U1; 27,U, 
aI. 
(—Cu S Cu aa 2 Cy ) 
_ mc 
— 2k,’ Ss ee 
B= i _ q : 4 
U; ml Oa? 
E = — 3h" GS -OC, 
vin U1 7 2ml 7 Oa," 
’ Ww! M,’Za,, q' Se C 
— se ay =i U; sion l, c= Ma) -+ 
(=) Ci, Ctg) 
U; 27m 
k,'M,’ ‘Se OC 
D= —K/+——t =" “(- “) 
l l 2/, Oa,” 
(=) (1/2)(0°C;,/Qa1”)(— Cag) 
U; 2T,m 
k'M,’ ‘Sc: f =F OG.) 
aoe K.! 4 i q ( ) 
l 1 8 3! Oa, 
(<=) (1/3!) (08C,,/Oa;,*)(— Cuq) 
Uy 2T,m 


Before proceeding with the solution of Eq. (9), it is 
This 


is done in the usual manner by considering the solution 


desirable to establish the stability of the system. 


in the neighborhood of the critical points of the equa- 
As- 


suming a stable system, it is desired to find a solution 


tion and identifying the type of the critical points. 


to the equation 


I+ f(I)dI + g(d)d = 0 (10) 
subject to the initial conditions, 9(0) = ao, 8(0) = 0, 
where, in particular 

fv) =At+ BI+ EY? 

gv) = C+ DI + FY 
and without loss of generality, aj is taken < 0. Un- 


fortunately, the exact solution to this equation cannot 
be expressed in terms of the known functions of mathe- 
OorB=E= 


However, a solution to any desired de- 


matics (except in the cases A = B = E = 
D = F= 0). 
gree of accuracy can be obtained by the methods of 
numerical integration. This procedure is tedious when 
done by hand and is best accomplished by means of 
automatic computing machines. The latter method is 
expensive and the additional cost cannot always be 
justified. It is therefore desirable to find an approxi- 
mate solution in terms of the tabulated functions of 
mathematics. 
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An iteration or perturbation procedure would appear 
to be in order. However, because of the nature of thy 
nonlinearities assumed, the first-order approximatioy 
is not sufficient. Furthermore, since g(3)d is of th 
third power in J and /(#) is of the second power in 9 
the differential equations for the second and _highe 
approximations become increasingly complex, and thy 
numerical application of the results to any specific ex 
ample becomes even more tedious than numerical] 
integration, while the results are less accurate. Fo; 
these reasons, this approach has been abandoned. 

A step-by-step approximating procedure has beey 
resorted to. The equation is linearized over a closed 
interval of the amplitude of #, the linearization being 
based on a time average of the response over the in 
terval. The constants used in the linearization ar 
changed for each step. In the following section, this 
method is developed and explained. 
mathematical rigor and _ relies 


The analysis is of 
necessity lacking in 
heavily on physical intuition and analogy to the well 
known linear oscillator. 

The solution to the equation 


A) + nd + w,"3 = 0 (11 
where yn and w, are constants, subject to the initial 


conditions, 


(0) a, (0) = 0 
is well known to be 
I = ay(w,/we ">" cos (wat — 2) (12 
where 
w = [w,? — (n/2)2]'7 
GQ = tan— [(»/2)/w] 


and the first derivative is given by 


od = —apw,(w,/wie °"~” sin wl (13 


If values of w,” and are selected to replace g(#) and 


f(d), respectively, in Eq. (10), the approximate solu- 


tion can be written as Eq. (11), which will be valid for 
sufficiently small values of time. It is noted that, in 
the theory of small oscillations, w,? = C and y = A, 
but this approximation is unsatisfactory for the magnt- 
tude of values of 3 under consideration here. 

Before selecting a value of w,”, it is observed that in 
the linearized case, the damping factor, n, has a second- 
order effect on the frequency, w. For example, i 
(n/2)/w, = 0.6--i.e., the system is 60 per cent criti 
cally damped—w = 0.S0w,—i.e., the frequency is ‘0 
per cent of the undamped frequency. If 
(n/2)/w, = 0.2—1.e., the system is 20 per cent critically 
damped—w = 0.98w,—1i.e., the frequency is approxi 
mately 9S per cent of the undamped natural frequency. 
In the cases that arise in practice, the missile is con 
siderably less than 60 per cent critically damped; ™ 
fact, it is usually of the order of 25 per cent or less. 


natural 


Since it is desired to find # as a function of time, it would 
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gem reasonable to select a value of w, corresponding 
to the frequency of the undamped system. 

Accordingly, consider for the moment /(3) = 7 = 0. 
Now Eqs. (10) and (11) become, respectively, 


i+ g(d)3 = 0, BO) = a BO) = 0 (14) 
§+ 0,20 = 0, 80) = a, d(0) = 0 (15) 
Bq. (14) can be solved as follows: 
bd = (d/dd) (87/2) = —g(d)d 
Integrating with respect to d, 
(16) 


6? 79 

v- . "4 
= — | g(x)x dx = G(ao) — G(e) 

9 


If (8) is a polynomial in 3, G(ao) — G(#) will also be a 
polynomial in 3, and, hence, 3? = P(d), where P(d) is 
again a polynomial in #. Writing 

[P(a)]” 


3 = d8/di = (17) 


7) 
| dx /{P(x)]'”? 


which defines 3 as a function of f. 
Similarly, Eq. (15) can be solved, and its solution 1s 


7% dx 
Jan Ola — x) 


Letting the interval under consideration be ay < 3 < J; 


and integrating, 


t= (1S) 


{ = (19) 


0, then from Eq. (18), 


ri dx 
a fee 


where ay < 
ty = 


and from Eq. (19), 


us 
; ov dx >) 
t = rs : = 
Jao Wyle” —~ X*) Wy» 


Equating /; = ¢,;’ and solving for w,, 


(1/2) — sin—!(0)/ ao) 





emi 


it in 
ynd- 
if 
riti 
NU 
If 
ally 
XI 
cv 
on 
in 
SS. 


ild 


oO, = » 
/ - dx 
n - P(x)]| 


and the approximate solution to Eq. (12) can be writ 
ten: 
(20) 


Jj = 


ay COS Wy 
for which 

J} = aw, Sin w,f (21) 
Consider now Eq. (10) with /(#) # 0. Forming an 


energy balance, 


i *d *d 
9 + | g(x)x dx = - / fd) 3 dd 
= ao a 


. 


Where the integral on the right-hand side represents 


CHARACTERISTICS 


AND DYNAMIC RESPONSE 307 


the energy dissipated during the time # varies from 
Letting yg; represent this integral, 


“y t 
/ f(d)I dd = [019 ar (22) 
SJ ao J/0 


Similarly, for the linearized system, 


ve W»" : , "0 
| - (3 Q0*) = n ddd 
Es.’ = 1% I v(r) dr 

Jd 


a to d. 


Fais 


(23) 


Substituting the approximate solution to Eq. (14) 
as given by Eqs. (20) and (21) into Eqs. (22) and 
(23), 
Eais. = O07Wy [ f(ap cos x) sin? x dx 
+0 
a Pe ' w,f sin 2w,f 
Kais = A wW,N -_ ; 


Now at the end of time /;, equating the energy dissipated 


in the two systems--i.e., Haj, = Dajs.’—and solving 
for n, 
/ f(ap cos x) sin*x dx 
— 0 
(w,f:/2) — (sin 2w,f,/4) 


The approximate solution to Eq. (10) valid in the time 
interval 0 < ¢< f; for which ay < 3 < 3; 1s, then, 


J = ay(w,/wye -"~ ‘cos (wt — 2) (24) 
where 
sin! (2;/ a) 
Wy of) ix 
[ ax 
: # [ P(x)] 
ft, = sin! (8)/ao)/a, 
wnt aye (25) 
f(ap cos x) sin? x dx 
-’0 
n = 
(w,t,/2) — (sin 2w,t,/4) 
WwW = lw, aes (n 9)2] 
Q = tan |(n/2)/w 


Differentiating Ey. (24), the value of J is found to be 


ov = — apw,(w,/we cee)? sin wl (26) 
It is observed from Eqs. (16) for the linearized system 
that, as 7 > 0, w > w, and 2 + 0; hence, 3 > — aqw, 


The 


quantity apw,| represents the maximum velocity of the 


sin w,f, which agrees identically with Eq. (21). 


undamped system. 
For the nonlinear case with zero damping, the angu- 
lar velocity is given by Eq. (17) as 


J = [P(9)]’ 


However, letting 7 > 0 in Eq. (26), the angular velocity 


1S 
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o = — aoe, sin Wyl 


which in general will not equal [P(#)].'* The maxi- 


mum value of 3 for the undamped nonlinear system is 
obtained when # is equal to zero in Eq. (17), 


' ‘ I/o 
Umax. si [G(ao) : 
For the linear system, the maximum value of # is 
‘ ' 
Cmax. = | AoWn| 


Hence, if the quantity — aqw, in Eq. (26) is replaced by 
- ve ° ° . é 
[G(ao)|’*, the approximate expression for 3 becomes 


I = [G(ao)]'"(wn/w)e"™ sin wt (27) 
which in the limit as 7 ~ 0 becomes 
i = [G(ao)]'”? sin w,f 


It should be noted that the expression for 3 given by 
Ey. (27) is not the time derivative of the approximate 
solution for 3 given by Eq. (24), but it has the same 
general character and gives the proper value of J,,,, 
in the limit as 7 > 0. Because the approximate value 
of 3(t) obtained in this way is not the time derivative 
of the approximate value af 3(/), the functions J(t) and 
v(t) should be treated as unrelated functions. 

Since g(#) is, in general, nonconstant, it is clear that 
the appropriate value of w, will be a function of the 
amplitude of the oscillation. If the system is not con- 
servative—i.e., {(8) # O--this amplitude will change 
with time, and after time /. its magnitude is given by 


| F | | 
v; = | Xo( Wy, @)| é 


-(n/2)t 


The general procedure outlined above can be continued 
for another interval 3,’ < 3 < Jd, using the reduced 
value of the oscillation amplitude to determine the 
value of w, and modifying the procedure to account for 
the different initial conditions for the second interval, 
since J(t2) will, in general, not be zero. This solution 
will be valid for t2 < ¢< fy. 
process, the solution can be continued, and, by selecting 
the intervals of time sufficiently small, the approximate 


Following this step-by-step 


solution can presumably be made to differ from the exact 
solution by as little as desired. Furthermore, the ap- 
proximate solution for the undamped case given by 
Eqs. (24) and (27) by letting » ~ 0 and the exact solu- 
tion in the undamped case given by Eqs. (17) and (18) 
agree at ¢ = 0) and at ¢ = f, irrespective of the magni- 
tude of f. This suggests, by analogy to the damped 
linear system, that in the case of small damping the 
approximate solutions should be reasonably good for 
fairly large values of fo. 

It is convenient to select f2 such that 3; = 0, since 
this simplifies the initial conditions for the next in- 
terval. With this choice of f:, the approximate solu- 
tion can be written, for 0 << f< fe, ao < 8 < 0, 


Wy =e ‘ 
BF = a ( y e °"% cos (wt — 2) (28) 
w 
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Fic. 3. General character of pitching moment vs. angle of attack 
for missile ‘‘A.”’ 
; ry, {On ee 
vd = [G(ao)] ‘( Je /*)" sin wf (29 
w 
where 
w/2 
W, = ~ (30 
0 6d 
1/o 
- ay 
} w/2 
y = (a cos v) sin? y dy (3 
Tv 0 
9 9 9}! ) 
w = [w,”? — (n/2)?] (32 
Q = tan (n/2/w) (33 
and when 
l = ty = |[(r/2) + Q)/w 


= . 1/9 Wn (n/2){[ (w/2 
I(t.) = [Glao) | e 
Ww 


o = me (n/2)t[ (#/2)4+-2)/@} 
l — (y 


For the next interval, it is convenient to choose fs such 
that J varies from J(f2) to 0. This gives the new values 


of w, and n, 


T / f dx 
Wy, == / We 
2/ 0 [ P(x) ] 
- eo — 
mn = f(a; cos y) sin? y dy 
T 0 


where 6, is the amplitude of the peak of the undamped 
oscillation on the opposte side of zero from a;. Hence 


G(6,;) = Gla). If g(3)d is an odd function of #, then 


61 =e Gi. 

Letting r = ¢ — tf, where ¢ = 0 corresponds to 3 = 
a, 3 = 0, the second interval of the oscillation is given 
by 

ny —(m/2)r o: 
v= 56 e SIN wT 
Ww 
F F Wn, (m/2)7 
v = D(te) e cos (a7 + Q)) 
Wy] 
where 
» 1 
oO = lw, (m1 )) 2] 
Q, = tan! [(m/2)/o] 
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The process is then continued for a third interval such 
that & varies from #(t;) to 0; hence, this interval is 
exactly similar to the first. The fourth interval is simi- 
lar to the second, etc. 

The procedure outlined above can be followed for any 
reasonable functions /(#) and provided the 
amplitude of the oscillations considered is such that the 


It is found in practice that the 


g(a), 


system is stable. 
pitching moment vs. angle of attack, which determines 
the function g(#)3, can be satisfactorily represented by 
a cubic polynomial. This has the advantage of lead- 
ing to elliptic integrals of the first kind for the deter- 
mination of w,, and these functions are tabulated. If 
higher order polynomials are used, the integrals will 
have to be computed. 

Seven representative numerical examples have been 
worked out by the method outlined above. These 
were taken from actual missile configurations and in- 
clude high- and low-altitude operation at different Mach 
Numbers and angles of attack. The initial conditions 
were chosen to correspond to step-function control- 
surface deflections necessary for turns of the order of 
from 1 to 5 g’s. The general character of the moment 
and lift curves vs. angle of attack, together with their 
mathematical approximations by cubic polynomials, 
iS presented in Figs. 3, 4, and 5. 

The seven examples can be described by the values of 
the coefficients A, B, C, D, E, and F and the initial dis 
placement ay. These are presented in Table 2. 

The details of the calculations for these cases are pre- 
sented in Appendix (I). The results of the analysis in 
comparison with numerical integrations obtained by 
IBM punched card techniques and the classical linear 
ized solution are presented in two ways. The responses 
in both 3 and # are plotted against time and are pre- 
sented alternately with the appropriate plots of the re- 
storing moment —g(#)3 and the damping function 
(8) against displacement in Figs. 6a to 12b. Phase 
space diagrams on which the trajectories of the solu 
tions are plotted are presented in Figs. 13 to 19. 


(IV) CONCLUSION 


The response to a sudden angle of attack in the case 
of a nonlinear moment curve is found to be qualitatively 
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General character of pitching moment vs. angle of attack 
for missile ‘‘B.”’ 
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General character of lift vs. angle of attack for represen- 
tative missiles 


Fic. 5 


similar to that for a linear moment curve, though quan- 
titatively the difference is marked. This difference in 
the examples treated, which are considered to be rep- 
resentative, indicates that the classical linearized solu- 
tion in such cases leads to large errors. The approxi- 
mate solution derived gives results that are consider- 


TABLE 2 
Values of the Coefficients in the Nonlinear Equation of Motion [Eq. (39)] and Initial Displacement for the Seven Examples Considered 


Damping Function 
Coefficients 


A B E 
Missile ‘‘A’”’ 
Example I 10.0 0 17.0 
Example II 1.32 0 2.24 
Example III 9.0 3.5 0 
Example IV 1.3 0.5 0 
Missile _ 
Example V 2.65 1.0 0 


Exampie VI 5.1 2.0 0 
Example VII 0.37 0 0 


Restoring Moment 


Coefficients Initial 


( D I Displacement 
120 0 17,000 —0.10 
14 0 1,950 —0.10 
700 1,300 — 74,800 —().08 
80 150 —8, 600 —().08 
60 300 0 —0.10 
172 — 1,580 0 —0.10 
11 — 102 0 —0.05 
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Fic. 16. Phase space diagram for Example IV. 


ably more representative of the actual behavior and ap- 
pears to be adequate for design purposes. If extreme 
accuracy is required, the methods of numerical inte- 


gration should be used. 


APPENDIX (1) 


Application of Approximate Method of Solution to 

Specific Examples 

In the application of the approximate method de- 
veloped to actual cases it is assumed that f(#) = 
A + Bd + Ev’ and that g/d) = C+ Dd + Fd’. The 
determination of the effective damping constant is 
straightforward by application of Eq. (31), and the re- 
sult is 


n = A+ (4B/3mr) a9 + (1/4) ao? 
The algebra in the determination of w, varies some- 
The 
different cases encountered in the examples cited in 


what depending on the values of C, D, and F. 


this paper are treated separately below. 
Case I.—-(Examples I and II): D = 0; 
ay < 0; hence, 


C,F> 0, 
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29) = C+ Fer? (I.1) 


From Eq. (16), 


ra) = ¢ = 2 f (Cx + Fx*) dx 


F ; 2 (1.2) 
P(d) = - (ag- — #*) F +a’?+ 
: 
| 
Let ¢ = &8/ao and a? = ay?/[ao? + (2C/F)], then 
0 <a? <1. Now Eq. (1.2) becomes 
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oF 2C 
Pia) = av?Pi(t) = — (a + © )o — ¢*)(1+a%? 


and 


“i dx 1 ‘ 
[ [P(x) - [(F/2) ao? -+- cx 


~1 dt 
o [1 — &) (1 + a%?)]” 


Let ¢ = COS ¢, 
“= dx l 
J [P(x)]* (C+ Fav?) 
ats dg K(y) 
f (1 — k? sin? g)'”? ws (C + Fay?)’ 
where 


y = sin-'k = sin ( : = » 
2 Fay” + 2C 


and A(w) is the complete elliptic integral of the first 
kind. From Eq. (30), 


On, = (x 2) (C + Faw?) 3 K (yp) 


lA 


The approximate solution for the first interval 0 
t< T), a < 3 < Ois given by 


Wno (n,/2)t : 
g = a es™*" COS (wot — Q) (1.3) 
Wo 


F Wg (n,/2)t . 
v= Im, e SIN aol (1.4) 
wy 


where the constants are given by the following equa- 


tions on setting 2 = 0: 


ny = (x/2) (C + Fa;”) */K (y;) (1.5) 
¥; = sin— [Fa,?/(2Fa;2 + 2C)]' (1.6) 
n = A + (4B/3mr)a; + (£/4)a;? (1.7) 
wo: = [wr — (nm 2)2]'/7 (I.8) 
Q;, = tan! [(;/2)/a;] (1.9) 
du, = [Ca;? + (2D/3)a;2 + (F/2)a;*]'* (1.10) 
Now 7} is the value of / when 3 = 0; hence, 
Ti = [(/2) + Q%]/on (1.11) 
3(7:) = dye7 /?™ (1.12) 
a, = —ae (ye/2) Ts (f.i3) 
Since g(3)3 = Cd + FY, which is an odd function of 


J, the amplitude of oscillation for the undamped sys- 
tem is symmetrical with respect to 8 = 0. The solu- 
tion for the second interval, 7; < ¢ < TJ», #(T) > 
i> 0,0< 8 < a, is given by the following equations, 
where r = ¢ — 7): 


Wn —(m/2)r . 
v= a é SIN WT 
@, 


(1.14) 


Ow 
— 
Ww 


DYNAMIC RESPONSE 


; eee Wn, ade ; Es 
3 = (Tr) (**) e wine OO8 (ay7 + Q)) (I.15) 


® 


The constants are given by Eqs. (1.5) through (1.9) on 


setting 7 = 1. Now 7> is the value of ¢ when 3 = 0; 
hence, 

T2 = T1 + $[(#/2) — Q]/o} (1.16) 

a = ae ” a se. (1.17) 


For the third interval, 7. < t < 73, as > 3 > O, the 
solution is given by the following equations, where 
r=t—-— T>2: 


[ (2/2) . \ 
dF = a»|w,,/we)e °™~’" cos (wer — Qe) (1.18) 


(2/2 


3 = Iar,(Wn, We) * sin WoT (1.19) 


the constants being given by Eqs. (1.5) through (1.10) 
on setting 7 = 2. Now 7; is the value of ft when 3 = 
0; hence, 


Ts = T2 + }[(#/2) + Q2]/ar} (1.20) 
3(T3) = due '* 2)(Ts Ts) (1.21) 
aa=- - Qe (m2/2)(7 T2) (1.22) 


For the fourth interval, 7; < ¢ < 74, 3(T;) < 3 < 0, 
0 > 3 > ay, the solution is given by the following equa- 
tions, where r = ¢t — T3: 


F = a3(wp,/w3)e °™ “T sin wer (1.23) 


3} = 3(7T3) (w,,/ws)e\"”~’"” cos (w3r +3) (1.24) 


the constants being given by Eqs. (1.5) through (1.9) 
on setting 7 = 3. Now 7, is the value of ¢ when 3 = 
0; hence, 

Ts = Ts + t[(w/2) + Q3]/ws} (1.25) 


3/2)(Ts—Ts sii 
a4 = ase" : (1.26) 


The solution for the next four intervals is obtained in 
exactly the same way as for the first four. 
Case I]._-(Examples III and IV): C, D > 0; F, 


ay <0; hence, 
g(30) = C+ Dd + FX 


From Eq. (16), 


P(3v) = -—2 I (Cx + Dx? + Fx’) dx 
P(3) = = (3 — ay) (8? + dV? + cD + d) (1.27) 
where 
b = ay + (4D/3F) 
c= ap + (4D/3F)a + (2C/F) 
d = a’? + (4D/3F) avy? + (2C/F) ao (1.28) 


Factoring, 


P(d) = (—F/2) (8 — ao) (& — By) (8 — Yo) (8 do) 
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where (po, Yo, 69 are the roots of the cubic, 


e+ bP +cd+d = 0 


Since the system is stable under the given initial condt- 
tions, Bo, Yo, 59 are real, By < ay < 0 < 59 < yo, and the 


oscillation in # is in the interval ay < 3 < dy). Now 


“0 dx l 
| [P(x)}"* ~ (—F/2)" ‘~ 
"0 dx 
/ I(x — ao)(x — Bo) (wv — Yo) (% — 60) 


» 


ss aoa F 2) (ay Yo!) (Bo - 60) 


yo dg 
J0 (1 — Ro? sin? ¢) 


(1.30) 


where 


(ap — 50) (Bo — Yo) | 
Ww = sin-! ky = sin”! : 
(Bo — Oo) lao — Yon) 


. | (Bo — 5p) |“? 1.39 
0 = sin = (1.35) 
Bo( ay — 06) 


and F {W, go} is the incomplete elliptic integral of the 
first kind of amplitude, gp. 
By Equation (30), 


(1 1) ((—F 2) (ao - Yo) (Bo 60) | 
F} Wo, vo} 


For the first interval, 0 << ¢< 7), ag < 3 < O, the solu 
tion is given by Eqs. (J.3) and (1.4); the constants are 
determined by Eqs. (1.7) through (1.10) on setting 7 = 
T, and 3(7)) are ob 
Since 


0 and Eggs. (1.31) through (1.33). 
tained from Eqs. (I.11) and (1.12), respectively. 
g(#)d is not an odd function of 3, the amplitude of os- 
cillation of the undamped system is not symmetrical with 
respect to ) = 0. Hence, for the second interval 0 < 
0d < 6, where 6; is the lesser of the two positive roots 


of the equation, 


° + byw? Sa Cw + dy = >) = 


(3 By) 6:) (& — 1) 


bo < a <0< 61 ~ Vy 


The quantities 5, c, and d are obtained by substituting 
a = ae ™™ for ay in Eqs. (1.28). Hence, for the 
second interval, 7; < ¢< 7», d(7;) > J > 0, O0O< v< 
62; letting 7 = ¢ 7\, the solution is given by Eqs. 


(1.14) and (1.15) where, now, 


(1/4) [(—F/2) (ay 


Kivi) — Fh, ¢1) 


4) (21 6) 


» 


AERONAUTICAL 


(1.29) 


SCIENCES MAY, 1952 
¢1, Yi are determined from Eqs. (1.31) and (1.32) upon 
replacing ao, Bo, Yo, 60 by ai, B1, ¥1, 41, respectively, 
and Eqs. (1.7) through (1.10) on setting 7 1 apply 
7, is obtained from Eq. (1.16), and 6 is now given by 
do = be \"/* rents (1.35 
For the third interval, 72< ¢< 73, 6: > 3 > 0, the soly 
tion is given by Eq. (I.18) on replacing a2 by é& and Ey 
(1.19). 
and replacing a. with 6) in Eqs. (1.7) through (1.10 


The constants are obtained on setting 7 = 


w,, is obtained from Eq. (1.34) on replacing a, 8), y;, § 


Yi and g with az, Bs, y2, 6, Y2, and ¢ge, respectiveh 
y2 and ge being determined from Eqs. (1.31) and (1,32 
upon replacing a, Bo, 
tively. 7; 
(1.21). For the fourth interval, 7; < ¢< 74, d(7 

v< 0,0 > 3 > ay, the solution is given by Eqs. (1.23 


Yo 69 with Q2, Bo, 2» Od», respec 


and &(7;) are given by Eqs. (1.20) and 


and (1.24). a; is the numerically smaller of the two 
negative roots of the cubic 
> + ded? + c33 + dz = 0 = (8 — ay) X 

(Jj - fe} (v - 


and d; are obtained from Eq. (1.28) on r 
n2/2)(T3— T2) 


where bs, Cs, 
placing ay by 6; = de The constants are 
determined by setting 7 = 3 in Eqs. (1.7) through (1.10 
and by replacing a, Bo, 59, Yo With as, Bs, 63, Ys, respec 
tively, in Eqs. (1.31) through (1.33). 74 and ay are 


given by Eqs. (1.25) and (1.26). 


Case III.—F = 0, g(#) C+ DI. Now, 
P(8) = 2 / (Cx + Dx?) dx 
= (—2D/3) (8 — a) (8 — Bo) (V by (1.36 
where 
ao + (8C 2D) 
By 
2 
Sao 
yr ; 
( ay + (3C/2D)} § 
1.37 
- ao + (3C/2D) |] 
00) 


j bay 
l l - 
( | ay + (SC = | 


Since the system is stable under the chosen initial 
conditions, By) and 6 are real. 

Case Illa. IfC, D> 0, Bo < a <9 
6) > O, the solution is obtained in the same way as 1 
Case II with the exceptions that: Eq. (1.33) is replaced 


(Example V): 


by 
(r/4) [(2D/3) (60 Bo) 
a { ’ I »” 
Fi Wo, Pos 
Eq. (1.31) is replaced by 
Wo = sin! [(69 — ao)/(b0 — Bo (1.39 


and Eq. (1.34) is replaced by 





D 
in | 


Eq 


Eq 


all 


st 


tre 
lin 
of 
de 
da 
tv 
wl 
tit 


in 





) upon 
‘tively, 
apply 
n by 


(1.35 
€ Soly- 
id Eq. 


{= ) 


vely 
(1.32 


CSpec 


and 


mn re 
S are 
I.10 


spec- 


» 
18) 


tial 


; in 


“eC | 


»C) 











AERODYNAMIC 


(r/4) [(2D/3) (6: — Bi))” 
big ath + (1.40) 
; (A (yy) — Fivs, C1} 
Case J11b.—(Examples VI and VII): If C > 0, 
D <0, Bo) > & > 0, aw < O, the solution is obtained 
in the same way as Case IT with the exceptions that: 


Eq. (1.33) is replaced by 


ay) 


(r/4) [—(2D/3) (Bo 


i Fi Wo Go} 


1 (3 _ “3 i I 
E & 
By — ao 


Eq. (1.31) is replaced by 

Yo = sin 
Eq. (1.32) is replaced by 
sin 


Ar 


go = 


and Eq. (1.34) is replaced by 
»1 
(wr /4) (2D 3) (8; — ay) | 


K(y) 


(1.44) 


APPENDIX (IT) 


Effect of Nonlinear Pitching Moment Curves on the 
Determination of Stability Derivatives from Flight- 
Test Data 
One of the techniques used for the determination of 

stability derivatives from flight-test data is the analysis 

of the transient response to a step-function input in con- 

trol-surface deflection. It is customary to assume a 

linear system. From a measurement of the amplitude 

of the oscillation vs. time, the logarithmic decrement is 

The zeros of the oscillation determine the 

The combination of these 


determined. 
damped natural frequency. 
two then defines the undamped natural frequency, 
which to all practical purposes determines the quan- 
uty 


M,,’ = (q’ Sc I.) (dCyy da) 


If the moment curve is a nonlinear function of the 
angle of attack, the ‘‘frequency”’ of the oscillation varies 
with the amplitude. In this case, the usual procedure 
of measuring the zeros of the oscillation and averaging 
them will not determine the value of dC,,/da at the 
The 


comparison of the classical linearized solution with the 


trim position but some sort of average value. 


numerical integration for the seven examples considered 
in this paper clearly indicates that the value of dCy,/da 
obtained in the usual way will not be very meaningful 
if the nonlinearity is at all severe. 

Alternatively, an attempt to determine the nonlinear 
moment curve, using the approximate solution for the 
response, may be made. The procedure to be outlined 
requires that angle of attack, normal acceleration in 
the pitch plane, and pitch rate be measured as functions 
of time and assumes, of course, a knowledge of Mach 
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Number, velocity, and dynamic pressure in order that 
the moment curve may be reduced to coefficient form. 
In addition, it is necessary to assume the general shape 
of the moment curve. (This is analogous to the as- 
sumption of linearity usually made.) 


Consider first a moment curve of the form 


—g(v)3 = —(C + Fd’)d; C, F>0 


Assume also that the damping function is constant, 
f(3}) = yn. From the records of angle of attack vs. time, 
normal acceleration vs. time, and pitch rate vs. time, 
the damping constant can be determined. The initial 
amplitude, ao, of the oscillation in angle of attack is also 
measured. Designating 7) as the first zero of 3 and 


using Eq. (I.11), 


where 2; = tan! [(n/2)/wo]. If the damping is small 
compared to critical damping, 2) may be determined as 


Q) = (n/2)/w,,, and w,, can then be determined from 


T ag \° n/2 (n/2)? 
rig ad [—e 


n, ®,, and ae are determined and can be substituted in 


Eq. (1.5), 


Wn, = (r/2) (Fay? + C) K (Yo) 


where 
Yo = sin—! [Fao?/(2Fao? + 2C 
which determines a relation between F and C. 
Now the normal accelerometer measures the quan- 
tity 


n=w— U6 - g cos O; 


where 6 is the pitch rate. Hence, from the normal ac- 


celeration and pitch rate records, the quantity 


v= (n/U) + 6+ (2/U) cos O; 


is determined as a function of time. Designating 7» 
as the first zero of J and proceeding in the same manner 
as before, using now Eq. (1.16) instead of Eq. (1.11), 
T 
= il 
1 ( Ts = 7) 


etn) ~tts +27 
1(72 — 7)) (T2 — 71) . y 


and a second relation between F and C is determined 


Wn; 


from the equations 


Wy, = (w/2) (Fa? 4 C)'?*/Kyy 


y= ' [ Fay?/(2Fa;? + 2C 


sin 


n/2)1 


ay = age 


Following this same procedure for the subsequent 
zeros of the oscillations in 3 and J, additional relations 
between C and F can be obtained. If C is plotted 
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against F, these various curves will intersect each 
other. Since the equations used to obtain the relations 
between C and F are only approximate, it cannot be 
expected that all of the curves will intersect in a com- 
mon point. However, depending on the accuracy of 
the data, they should define a small closed region in the 
F, C plane, and any point (/, C) in this region may be 
used for the approximate values of F and C. 

As a second example of the application of the approxi- 
mate solution to the determination of the moment 


curve, assume the curve to be of the form 
—g(v)3 = —(CV + Dd), C,D>O0 


and assume that the damping is constant. 

From the envelope of the angle of attack vs. time 
record, ay and 6) are determined as the initial values of 
the amplitude on opposite sides of 3 = 0. Proceeding 
as before, using now Eq. (1.38), 

1 ‘ 11/9 
(7/4) [(2D/3) (60 — Bo)] 
+r) 
J } Wo, vo} 


Wo — sin I [ (do > 
(2 = =) 
go = Sin 
Bo( a — 60) 


Observing that 


) 
Wn, 


a) (do = Bo) |’ ; 


ao) (8 Bo) (8 — bo) = 


-2D/. 3C., a 
a ") + oD ve — a+ oD ay 


and comparing coefficients Gy avéo/(ao — 4) and 
3C/2D = a + Bo — 4. This determines D and C. 
Repeating the procedure for the other zeros of 3 and J, 
additional determination of C and D may be made. 
An average of the different values may then be used as 
the correct values of Cand D. 

Obviously, a similar procedure can be used for other 


re) = — (ov 


assumptions of the moment curve shape. 

It is to be emphasized that the above procedure ap- 
plied to a nonlinear system, since it is based on an ap- 
proximate solution, cannot be expected to give results 
so accurate as those obtained for a linear system where 
the exact solution is known. However, it should yield 
a more useful result for a nonlinear system than that 
obtained from an assumption of linearity, since it ac- 
counts, in part at least, for the nonlinearity. Further- 
more, this method has the advantage in actual applica 
tion that large amplitude oscillations may be used; 
hence, the effect of instrument errors is reduced. In 
the application of the proposed method, there are sev- 


eral sources of error which are to be guarded against. 
The first case discussed, where the moment curve was 
assumed to be symmetric with respect to the origin, 
applies to the case of a symmetric missile trimmed at 
zero angle of attack. Consequently, the abscissa of the 
various instrument records is accurately known (after 
correction of § and wW records for effect of gravity, which 
is small). Furthermore, after 7 is determined, the 
envelope can be faired through the peaks of the oscil 
lation and the zero should bisect the envelope. This 
affords a check. 
be in the amplitudes of oscillation. 


Hence, the major source of error wil] 
Since the records 
from all three instruments will exhibit the same damp- 
ing, all three should be used in the determination of » 
The determination of J(f) requires the sum of the two 
instrument readings, u(t) and 6(t); hence, it cannot 
be expected to be so accurate as that of d(f). 

The case of a moment curve that is not symmetrical 
about the origin arises where a symmetrical missile js 
trimmed at a finite angle of attack. In this case, the 
oscillation has different amplitudes on the opposite 
sides from zero. The abscissa of the oscillations cannot 
be measured in this case and must be obtained from a1 
extrapolation back from the steady state following the 
transient. This is, of course, an extremely inaccurate 
procedure and will undoubtedly introduce large errors 

An attempt has been made to apply the above method 
to the numerical solution of a case with a moment curve 


of the type 
o(3) } = C¥ + Dv 


Errors of 10 per cent were assumed in the measured 
quantities, and it was found that such errors could re- 
verse the curvature of the moment curve, though the 
slope of the moment curve at trim was virtually un 
changed. Hence, it would appear that, in order to de 
termine a nonlinear moment curve from flight-test data, 
extremely accurate instrumentation will be necessary. 
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Aerodynamics of Blasts’ 


L. TING? anno H. F. LUDLOFFt 


New York Unwersity 


ABSTRACT 


The problem of a blast passing over the surface of an arbitrary 
flat structure is investigated. The pressure and density fields 
in the whole domain behind the advancing blast are obtained in 
explicit analytic form. As a simple example, the constant pres- 
sure lines (isobars) and the constant density lines (isopyenics ) 
are plotted for a strong blast of pressure ratio 7.308, passing over 
a wedge surface. Comparison with the interferograms taken in 
Bleakney’s laboratory shows surprising agreement The “slip 
stream”’ connected with the Mach reflection of the shock comes 
out asa thin strip, in which the density gradient is steep but finite 

To obtain an analytic solution, no cone field transtormations 
are used as previously, because they restrict the generality of the 
results. Instead, (a) a Lorentz transformation of the coordinates 
is introduced, (b) the solution is given in terms of a Possio in 
tegral, (c) a ‘“‘reflection”’ of the boundary condition is carried out, 
(d)a “lost boundary condition” is taken into consideration. 


SYMBOLS 


1p, A,;, Av, Ay = constants defined in Eq. (5.1) 

a, a = constants defined in Eq. (5.1) and (6.4 

BA = sound speeds of the undisturbed flow behind 
and in front of the shock, respectively 

D, D = differential operators defined in Eqs. (3.16) 
and (4.5) 

f = the wall surface 

WV, M = Mach Numbers of the undisturbed flow be 
hind and in front of the shock 

P; Pe = the pressures of the undisturbed flow behind 
and in front of the shock 

p = the disturbed pressure 

Vv = the reflected shock front (Mach circle 

S = the “separation line’ 

1 = the triple point 

t = time 

l = time in the system of Lorentz coordinates 

U, l = the velocity of the undisturbed flow behind 
and in front of the shock 

u = the x-component of the disturbed velocity 

= the y-component of the disturbed velocity 

(<,.¥ = coordinates fixed on the undisturbed flow 
behind the shock 

(z", 9) coordinates fixed on the wall 

(x, 7 = Lorentz coordinates 

Y = 1.4 for air 

6 = angle between the normal of the diffracted 
shock and the x-axis 

Ao, Ar, Az, As = constants defined in Eq. (5.1) 
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Ao, Ai, Az, As = constants defined in Eq. (6.4) 

R, po = densities of the undisturbed flow behind and 
in front of the shock 

p = the disturbed density 

y¥ = the form of the diffracted shock 

Q, 24, Qe = constants defined in Eq. (3.15) 


(I) INTRODUCTION 


I THE FOLLOWING, the pressure and density fields 
that develop behind blasts advancing over arbitrary 
flat surfaces have been investigated. Computations 
along similar lines have been carried out previously by 
Bargmann! (who assumed the flow field to be irrota 
tional) and by Lighthill* (who took the rotationality of 
the flow into account). The method to be presented 
here is more direct, and the results obtained are more 
general than in the previous investigations. Shocks 
passing over flat surfaces of arbitrary shape can be dealt 
with in such a manner that explicit analytic expressions 
result for the pressure and density fields in the whole 
domain behind the advancing shock front. In this way, 
comparison with shock-tube experiments can be readily 
carried out, and interesting details about the nature of 
the ‘‘slip stream’”’ occurring with the Mach reflection of 
the shock can be determined. [See Sections (V, 2) and 
(VI). 

Assume that an originally plane shock front of given 
intensity (strong or weak) advances over the surface 
f(x’) of a given flat structure (or a given thin airfoil) 
into still air of density pp and pressure Py (see Fig. 1). 
The velocity of the shock may be Uo, the speed of the 
air behind the shock is (Uy) — LU’), and its density and 
pressure there are p and p. The length of the wall dis 
turbance (or the chord of the airfoil) may be /. If the 
inclination of the body surface f’(x’) with regard to the 
direction of propagation of the shock is small, the shock 
front will end perpendicularly on the surface at any 
instant so that a curved shock front results, as well as 
a shock configuration that may be interpreted as Mach* 
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the body surface, so that it differs from the stationary 
Ug-U eat : ; ary 
G. — S| pressure distribution. Hence, the field in domain J 
may be treated as a time-dependent perturbation away ; 
from the state of (relative) rest existing in domain | R 
‘ A linearized theory can be derived by a formal perturba- the ¢ 
tion procedure consisting essentially of an expansioy the 
in terms of a parameter ¢, which can be interpreted as turb 
thickness ratio of the wall disturbance (or of the air. ¥ 
foil). equi 
\y ~ ‘ 
oer t 3 (II) THE LINEARIZED DIFFERENTIAL Eouartions 
J ; 
a< a . ° ‘ : 
Fixing the coordinates in the undisturbed flow behind Phe 
the shock, the x-axis coinciding with the plane wall and 
the y-axis being parallel to the undisturbed shock front. 
y the origin is so chosen that, at ¢ = 0, the shock front hits 
the leading edge of the bump. Then, the differential 
equations determining the two-dimensional unsteady ro- 
tational flow behind the shock are T 
ales ‘ mal 
the continuity equation: 
; an 
Pr t+ (pu)z + (pv), = 0 Fig 
: ; . Du 
the equations of motion: = — -)p;; 
(2.1 : 
Dy I 
SP 
Dt p 
Fic. 2. Domains of disturbance. | 
, : ; ) _ anc 
the adiabatic relation: (p/p’) = 0 
Dt 
reflection. Because of the curvature of the shock, the resulting from the assumption that both viscosity and | 
flow behind the diffracted part of the shock will be rota- heat conduction are negligible. 
tional and, for sufficiently strong blasts, supersonic with Now expand p, p, u, and v in terms of « 
regard to the body surface. Figs. 2a and 2b represent 
ce ems ° : = P pvoix. y - 6 p(y wy | 
conditions behind the advancing shock front depending P P+ ep (x, yt +E p(x, yD +... = 
on whether the airflow behind the shock is subsonic or pp = R + e€-p()(x, y, t) + e+ p(x, y, t) + 
. - ‘ - ls \ 
supersonic relative to the obstacle. In either case, the r (2.2 
° . ° ° ° ° ° =_ Ae - ’ { 2. 2 . ’ 
time history of the diffraction is depicted in an x-y-t- " EU WO Tr U(X, Wt) +... 
space, so that every cross section of the figure represents y = ¢-y(0(x, y, tf) + &-v@(x, y,f) +... 
the domain of disturbance at a given instant /, after : 
the incident shock has hit the leading edge of the ob- Inserting Eqs. .s-8) into Eqs. (2.1) and equating the 
stacle and a disturbance has spread with sound speed muaienentte he ae ponent at «, one eae, wh 
fe as : : approximation, 
throughout domain II up to the circular reflected shock. 
The coordinate system, x, y, ¢, is to be fixed in the un- po?) + R-u,) + R-v, = 0 (2.3 
disturbed flow behind the shock—.e., in the flow oc- R-u,) = —p, (24 
curring in domain I. With regard to this system, the R:-o,9) = —p,O (2.5 ( 
leading edge and the whole body surface f(x’) are mov- p) = c?-p,“ (2.6 


ing to the left, and the shock is moving to the right. 


; ; ; , Eliminating from the above equations three out of the 
In the subsonic case (a), the leading edge of the disturb- 6 


four unknown functions yields 


ance moves slower, and in the supersonic case (b) faster, 
than the perturbation produced by the shock diffrac- ee Oo” 1 0° bis “ 
tion. Therefore, in the latter case, domain III has Op = (= ay? gt ~) Pe eer 
been disturbed by the leading edge only, but it has not ; 
been reached by the shock disturbance. Hence, here and 

Ackeret’s formulas of stationary supersonic flow apply. (0/ot) (Og) = 0 (2.8 


Conditions in domain II are clearly nonstationary. where g denotes either u, v,“” or p™. 
The pressure in the vicinity of the body surface must be It is now necessary to find two initial conditions and 
expected to influence the speed and inclination of the one boundary condition for p“ to uniquely satisfy the | 


shock front, and this, in turn, will affect the pressure on wave equation [Eq. (2.7) ]. 
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AERODYNAMICS 


THE BOUNDARY AND INITIAL CONDITIONS 


(III) 


1) On the Incident Shock Front 


Relative to the undisturbed flow behind the shock 
the air in front of the shock and the wall are moving with 
the constant velocity —(U, U), while the undis- 
turbed shock front moves with the velocity + U. 

The disturbed shock front can be expressed by the 


equatic yn 


\ Lt + ew (y, 1) + Oe) (3.1) 
hen, the shock inclination is given by 
4 oe e-W,'(y, t) + Oe?) (3.2) 
Ov — . 
6A) + Ole’ (.3.03) 


The complete shock velocity, which is directed nor 
mal to the shock front at any point, may be split into 


an x-component U, and a y-component vs, where’ (see 


Fig. 3 
dx Ox ( .: / e ( 
S : t Ole r ey,’ (y, f Ve 
at Ol 
».4) 
and 
U's U's:tg0 U0 + Oe (3.5) 


The oblique shock conditions on the shock front 


x= Ut- e-y,\? 1 ()(€° 
are 
pq, Pog,,0 
P+ PGn” = Po F Podno™ 


1 (’) =F vy Po (3.6) 
=~“ 4 = - qu’ 4 
2 v= 2 nr 
q, = 
where 
g, = normal component of velocity of flow behind, 
relative to the shock 
Qn v — vs) sinéd + (u — Us) cosé = 
—U — eb) + eu + O(e) (3.7) 
q, = tangential component of velocity of flow be- 
hind, relative to the shock 
gq, = (v — vs) cos@ — (u — Us) sin @ = 
e-v') + O(e*) (3.8) 
g,0 = normal component of velocity of flow 
in front, relative to the shock 
Guw = —vs sind — [(Up — U) + Us] cos @ = 
—Uy — ef? + O(e*) (3.9) 
q, = tangential component of velocity of flow 


in front, relative to the shock 
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x= Utt e-W' (yt)... 








gq vs cosé+ [(Uy — l + l’s| sin 0 
—(U9 LU )-e-W,\? + Of« (3.10) 
Inserting Eq. (3.7) into Eq. (3.10) and Eq. (2.2) into 


) 


Eq. (3.6), the coefficients of like orders of « may be 


equated. The relations between the coefficients of « 
are equivalent to the customary normal shock condi 
tions. 

The relations between the coefficients of €! yield four 
equations between the first-order quantities that, in 
keeping with the approximation of the linearized theory, 
are to be fulfilled at x = Ut, y > 0. 
tions can be written in the form 


These four equa- 


c*- p(x vig > Gt 1 + Q)p 
\ Ut, y 0, t 3.1] 
0), 
u \ Ut, 4 0, ¢) p" 
; Re 
\ Ut,y >0,t 3.12) 
(2, 
(l ( Vv, (y, b) ) 
R! 
v = Ut,y >0,t) (3.13) 
(l Up, (y, 2) = v0 (x = Ut,y > 0,0) (3.14) 
where 
te (y — 1) (MW 1) 
My — 1) + 2)’ 
: 1 (3y — 1)M? + (3 — vy) (3.15) 
( 
=] > M4 = +2 
Q = — (1 — M*)/M 


By using differential Eqs. (2.3), (2.4), (2.6), and (2.7), 
), p, and y from the pre- 
In this way a boundary cond1- 


one can eliminate uu“, v 
ceding four equations. 

tion for p™ alone, to be applied at the shock front 
x = Ut, can be formulated. 


3.16) 


D,, p(x = Ut, y>0,t) = 0 


where D, , is a linear differential operator defined as 
O° - 
D, » = (Q + M + &M) >? + (1+ M? + 


32 
2MQ,)-c: + (QM? + M —M®Q)c? — (3.17) 
t re) 
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Fic. 4. The forms of f(x’) and f(x’) for a shallow bump. 


(2) On the Wall 


The boundary condition on the wall is, in first ap- 


proximation, 


v(x < Ut, y = 0,t) = (Un — U)f' [x + (Uo — Vt) 
(3.18) 


where the prime means differentiation with respect to 
the whole argument. 


lim p,“(x = Ut, y > 0, t) 


y—>O0O 


Similarly, Eq. (3.19) implies 


lim p,(x < Ut, y = 0, t) = —R(Uo — U)?-f"(Uot) 
x—> Ut 
(3.21) 
One may note that the two limits for p,‘” are equal 
if 


[4/(y + 1)]-U-Uy = (Uo — U)? (3.22) 


i.e., for a special shock strength P/P») = 7.3076 or M = 
0.5100. In general, the above two limits for p,“ 
are unequal. This indicates that p,‘” has a singu- 
larity at x = Ut, y = 0, as defined by Eqs. (3.20) and 
(3.21). 
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In general, f’(x’) is a bounded, sectionally continuoys 
function of x’ and has one-sided derivatives (see Fig, 
4). Since it is possible to find a continuous function 
that is close to f’(x’), one may assume that /’(x’) jg 
continuous and has continuous derivatives. This 
must be remembered if later on one tries to derive ex. 
pression (4.11). 

Using the differential Eq. (2.5), a boundary condi. 
tion for p“? at y = 0 can be derived from Eq. (3.18). 


py (x < Ut, y = 0, t) = —R(Uy — U)?f"[x + 


(Uo — U)t 3.19] 


(3) The Point of Intersection of the Shock Front and the 
Wall 

In deriving the preceding two boundary conditions 
[Egs. (3.16) and (3.19) |, the special relation 


v(x = Ut, y = 0, t) = (Uo — U)f'(Uot) (3.18 


which can be applied at the intersection point of both 
surfaces, has not been really used. Additional informa- 
tion will be obtained if Eq. (3.18a) is combined, on the 
one hand, with the shock relations (3.13) and (3.14 
and, on the other hand, with Eq. (3.18), which holds on 


the wall alone. From Eq. (3.14) it follows that 


vy,“ (0, t) = —[1/(Uo — V)]v (x = Ut, 0, t) = —f' (Ut 
Vy (0, f) = —Uo-f"(Uot) 
Differentiating Eq. (3.13) with respect to y, one ob- 
tains 
Q: py (x = Ut, y > 0, t) = —R(Us — U) dy (y, 8 
which implies 
RI hy ( Uo — U) - = 
= I ( l ot) 
Qs ; 
(3.20 


! 
=—  RUUo-f"(Udt) 


4) Onthe Reflected Shock (Mach Circle 


Since the disturbance that spreads after the shock hit 
the leading edge is small, the boundary line separating 
the disturbed region II from the undisturbed region | 
(or from region III in case this region exists) is, eXx- 
cept for terms of 0(€), the Mach circle; more explicitly, 
the curve 


r= c+ e-q(¢, t) + O(e*) 


where r = Vx? + y? and g = tg —'y/x (see Fig. 3). 
However, in a consistent linearized scheme boundary 
conditions have to be prescribed at the undisturbed 
Hence, one has to re- 


boundary—.e., here at r = ct. 


quire that on the outer side of the Mach circle—e., 
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pb = () where domain II borders on 


—— = (ct)*, 1 
domain I, and p\? = domain II 


borders on domain III. 
that, under these conditions, p“ varies continuously 
across the Mach circle; this means that, to the order 
., there is no discontinuity in p“ across the Mach circle 


p Ackeret ) where 


It is shown in Appendix (I) 


i.e., the reflected shock is of o(€). 

Actually, it is not necessary to state the boundary 
condition in such detail. One may simply prescribe a 
condition at infinity. 


p™ — 0, as V x? + y? > © (3.23) 

Then, the form of the solution of the wave equation, 
as assumed in Section (IV, 2), can be shown auto- 
matically to vanish in the whole domain I and to ful- 
fill the above-stated conditions on the Mach circle. 
In this way, the domain of disturbance is, only formally, 
extended, the boundary now consisting of the full 
length of the incident shock, the x-axis from — © to 
x = Ut, and the circular are at infinity. The shifting 
of the circular boundary to infinity is an artifice that 
makes the solution procedure developed in Section (IV) 


more lucid. 


5) Initial Conditions 


Obviously, the two initial conditions are 


pY(x< Ut,y => 0,t< 0) = O (3.24) 
bY(x< Ut,y >0,¢t< 0) = 0 (3.25) 
These conditions are actually implied by Eq. (3.19). 


(IV) DERIVATION OF THE ANALYTIC SOLUTION* 


FOR p"! 


1) The Lorentz Transformation 

The underlined equations [| Eqs. (2.7), (3.16), (3.19) to 
(3.21), and (3.23) to (3.25)] summarize the formulation 
of the mathematical problem. The fact that the com- 
plicated boundary condition (3.16) is prescribed on the 
plane x = lt may suggest the introduction of new inde- 
pendent variables (%, 7, 7) related to the old variables 


(x, y, ¢t) by the well-known Lorentz transforma- 
tion, 
= (x — Ut)/V1— M*, F=~y) (41 
= 
i = (ct — Mx)/V1 — M? f 
The plane ¢ = 0 corresponds to the plane x = Ut, and 
the wave equation 
Deg’? + pgs? — par? = 0 (4.2) 


remains unchanged. 


* The authors wish to express their appreciation to C. S. Gard- 
ner, Institute for Mathematics and Mechanics, New York Uni- 


versity, for his valuable suggestions 
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flow 


The boundary conditions become 


p — 0, as V rt? + y? —> a (4.3) 
py (& < 0, § = 0,7) = Re®Aof”[G(Xov + 7) (4.4) 
D,; 7p = 0,5 > 0,2) =0 (4.5) 
and 
lim ps? (#, 0, 2) = Re*Aof” (al) (4.6) 
g—> 0 
lim ps(?(0, 9, 7) = Ref” (al) (4.7) 
i—~>O0 
where 
: Us = Ml - c 
a = ; > No = l— WV - 
cV1—M ( l 
{ (Uo — l i } UoM 
Ao=-— _ -=— = 
; c v+1 « 
and 
“1 1 Oo O° O° 
D; i = = ae eh Se aa coe 
M,? OF OFO/ OF 
The initial conditions are now 
p) = p;™ =0 (for 7 < 0) (4.8) 
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bo 


(2) The Possio Integral 


In general, the solution of such boundary initial value problems as the one defined by Eqs. (4.2) to (4.8) can be 
solved in terms of what have been called ‘‘temporary sources” spread over a certain area in the x-t-plane character. 
izing the motion of the body (here, the motion of the wall relative to the air behind the shock). Such a solution 


can be written as Possio integral,*® 


. l -_ ps (é, 0, 7) 
pY(z, y, 1) = — Lf : - dé dr (4.9 
7 JR V @ — r)? — @ — €)* — ¥° 


in which the numerator p;"” of the integrand represents the “source strength’’ and the denominator represents the 
‘pseudo distance’’ between “‘source point” £, 7 and point x, y, ¢. The integration area R in the & — r-plane is con 
fined by the hyperbola 


f—7 = VG@— & + F' 


and the straight line r = 0. Looking at Fig. 5, it is obvious that the method is applicable only if p;‘") is pre 
scribed on the entire plane ¥ = 0. 

Actually, however, p;‘” is given for the left half of the plane ¥ = 0 (@ < 0), while it is unknown in the right half 
plane (> 0). Hence, the next step is to find an equation for ps“? (% > 0, ¥ = 0, 7), which will replace condition 
(4.5) prescribed on the plane @ = 0 and also produce the correct type of singularity for p;‘” at the origin as defined 
by Eqs. (4.6) and (4.7). 

Inserting the proper integration limits and using condition (4.4), the Possio integral (4.9) becomes 


Ww 


= Rc*Ao * ’ "9 ij | a T > hoe) J 
pu (2 <0; ¥, t) — / dr / dé = - - 
r 70 J2t— WV E—1)*—3 Yu — 7)? — & — 2)" — ¥ 


| "7— #24 “ET W/(F—7)?-5 ps (E > 0, 0, 7) 
dr dé : 1.10 
r Jo J0 V (i— 7)? — (% -— €)? tai 


Substituting Eq. (4.10) into Eq. (4.5) will provide the desired equation for the unknown function ps"? (£ > 0, 0, 7 


(3) Determination of py\(~ > 0, 0, r) 
By carrying out the proper differentiations with regard to ¢ and /, one obtains 


_ : ; : —Re®Apa?-H(—X) 7-3 "0 f° [a(r + RoE) 
D;,ip(& = 0, ¥ > 0,2) = dr i— 
PV AL e r)?=%f?2 : 2 1) - 


VG C= =F — 


*t 


Tv V 
| *V (F=7)?=52 D, ,ps(& > 0, 0, 7) | y Ao(r) + 2M-Ay'(r 
dr dé = — dr = (4.11 
7 Jo 0 Vi-r?-#-F So Vii -— 7)? — 
where 


H(x) = (1/M,?) — 2Mx + 2» 


Ai) = px(& = OF, F = 0,2) — p,O(E = O- F = O, 2) 
Aj’ (2) = (d dt) Ay (2) 
Ao(E) = Pog (4 = OF, F = 0,2) — fog (E = O-, F = 0,2) 


Since p is a solution of the wave equation, it is seen that the identity 
ps(& = 0+, ¥ = 0,2) + ps(# = O-, F = 0,2) = 2p,(& = 0, F = OF, 2) (4.12) 
holds true, which represents a kind of ‘‘mean value theorem’’ at the singularity. [For proof, see Appendix (2).| 
This relation and the property of the singularity defined by Eqs. (4.6) and (4.7) yield 
Ai(?) = 2Re?(u — Ao)f” (al) (4.13) 
or 
ps ( = OF, ¥ = 0,2) = Re?(Qu — Ao) f"(al) (4.14) 
According to Eq. (4.5), expression (4.11) is required to vanish. This will be the case, if 
D, i[ps(& > 0, ¥ = 0,2)] = —Re®Aga?-H(—Xo)-f" [aE — Xok)] (4.15) 


and 
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Ao(?) = —2MA'(t) = —4MRe?(u — Ao)af’’’ (a) 
or 
psx“ (x = 0 a y = 0, t) = ARe?| Ay(Xo + $17) —_ 1 Mu | f’"" (at) (4.16) 
Hence, the problem is reduced to obtaining a function ps" (% > 0, ¥ = 0, 7), which satisfies differential Eq. (4.15) 
and boundary conditions (4.14) and (4.16). It is shown in Appendix (3) that the solution has the form 
py?(F > 0,5 = 0,0) = Re Aif" [ai — Ak)] + Aof”’ [ai — Xet)] + Asf"[a(? — X3%)]! (4.17) 


where the expressions for d, \;, and A; are listed in the next section. 
(V) THE FINAL RESULTS 


1) Expressions for p, pp, u,v, andy" 


The results obtained in the previous section can be summarized as follows: 


aed Re?A, [*-3 , ’ S” (a(t + Xo£é)] 
pM, 9,1) = — dr dé : — 
uy 0 WE — os VG — 7)? — @ — $)* — § 


2,3 a . - 7. r= ee 

Re 2A; b— 4/g24 52 £+ Vv (f—1r)?-—7f2 ‘ {"lalr os \i€) " 
» dt dé : (5.1) 
i=1 T 0 J/0 V (i — rr)? — (% — ie se : a 


‘a 


where 
Uo ” MU,\ « 
@ = , Xo = (1 — M - 
cV 1 — M \ c Up 
Ay = —(Uy — U)?/e?, hi = Yo 
Nor + 2MX + (1/M_?2) 
1, — Ay 


Xo? — 2MX + (1/Mo?) 
\y and \3 are the roots of the quadratic equation 
2? — 2MX + (1/M)? = O 
and A, and A3zare the solutions of the two simultaneous linear equations 
Az + Az + Ai + Ao = —[8/(y + 1)]/(MU b/c) 
a S . a S MU, 
hoAe + A343 + A1A1 + AvAn = 2M — Ay 
7} i «£ 


The above solution is unique and satisfies the wave equation and all the initial and boundary conditions. 
By using the transformation (4.1), the disturbance pressure p“’ (x, y, ¢) can be obtained from p" (2, J, 7). 
From differential Eqs. (2.4) to (2.6) and the boundary conditions (3.11) to (3.14) the following expressions are 
obtained: 
The density variation: 
p(x, y, ft) = (1/c?)p™ (x, y, t) + W(x, y) 
= (1/c?)p (x, y, t) + (Q/c?)p™ [x, y, ¢ = (x/U)] (5.2) 


The x-component of the disturbance velocity : 


"% Q x 
u(x, y, tf) = — / dr:-p,) (x, y, 7) — p (x. y¢{ = ) (5.3) 
. R Jx/t : Re i l 


The form of the shock front: 


0, 
x= Ut — e: — ; / P(x = Ur, y, r) dr (5.4) 
R(Uy — U) So 3 


The y-component of the disturbance velocity: 
: I : 


» 


; as 7 
v(x, y, t) = — bV(x, y, r) dr + = | bY (Ur, y, 1) dr (5.5) 
. R | l f . R Jo Pr ¢ 
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(2) The Separation Line Between Rotational and Irrotational Flow—Its Physical Significance 
From Appendix (1) it can be seen that (on the shock, outside 7 = cf) 
p(x, y, ¢ = x/U) = Ofor x? + y? > x?/M? 

hence, also, 

% x . ° 

--pO |x, y, tf = —]) = OforxV 1 — M’?< My 

l 

Taking this into account, it follows from Eq. (5.2) that the lines of constant density (isopyenics) coincide with the 


lines of constant pressure (isobars) for x/y < \//V 1 — M?, but deviate from the latter for x/y > M/V 1 — jp. 
i.e., inside the area OTB of Figs. 6 or 7. 

In other words, looking at Fig. 7, the flow inside OTB is rotational, outside, it is irrotational. This is in aecord 
ance with Kelvin’s theorem, since no particles cross the separation line OT. Since the reflected circular shock js 
of o(e), in a first-order theory particles can acquire vorticity not by crossing the reflected shock but only by cross 
ing the curved portion of the incident shock TB. These particles stay inside the area swept over by TB—ie., 
inside the area OTB. 


The form of this separation line can be represented by (see Fig. 6) 


yp = €- S (6, 7) + O(e?) 


where 
yv=x-V1l—-M—y-M, oc=x-M+yV1l—M 
SO, 7) = I [Vv 1 M?-u (x, y, r) — M-v(x, y, r)] dr 
J0 
ef | 
y=oVl1—- VW 
Since 
lim p(x, y, x/U) = 0 
<—y(M/v/1 —M2)= 0 


there is no discontinuity of density —1.e., no slipstream in the mathematical sense, at least not in first approxima 
tion (perhaps in higher approximation). However, it will turn out that in a physical sense there is a slipstream 
of finite width, as explained in Section (V1). 


(VI) APPLICATION OF THE THEORY 


a 


The shape of the wall (or of the airfoil) y = ¢€-/(x’), 0< x’ < 1, may be approximated by a polynomial 
f(x") = Gx’ + cox’? +... Get (6.1) 
The accuracy of the approximation increases with the order of the polynomial. While the case of higher poly- 
nomials (a circular are airfoil) is treated elsewhere,® one may here consider the simplest example, 


f(x") = x" (6.2) 


representing a corner of inclination e. Then, 
f'(x’) = 1 for x’ > O and f’(x’) = O for x’ < 0 
Furthermore, 


f(x") = © (for x’ ~ 0)) 


ss ” , , j (6.3) 
f 7" (") dx’ = 1 f ee 


Using Eqs. (6.3), the pressure integral (5.1) can be evaluated and yields 


Re? | : Thy , ' 
pY(x, v,t) = - Ao: N(x, y, t, — bo) + > AN (x, y, t, Ad? (6.4) 
Tra i=1 
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AERODYNAMICS OF BLASTS 


where 
; 1—AM NM — x! : 
N(x, y, t, A) = - cosh! 7 (for A > 1) 
Vr? - 1 V(t — Ax)? + (2 — Dy? 
= (1 — AM) (VBE — x? — y?/t — Ax) ~(forA = 1) 
1—AM Ni — x : 
= — -cos~! (for A < 1), 
V1l-—--» V (¢ — Ax)? + (A? — Ly’ 
and 
% — M ell Uo 
Xo = —— a=aVvl-M= 
| a Ao M Cc 
4 + M ot 
Xi = a (for ¢ = 1, 2, 3) 
1+ A.M 
Substituting Eq. (6.4) into Eqs. (5.2), (5.3), (5.4), and (5.5), one obtains, respectively, the density variation: 
R hy 
pY(x, 7, f) = - AoN(x, y, t — 0) + DS AiN(x, y, t, Ad) | + W(x, 9) (6.5) 
T i l 


where 


QR x = x 
W(x, y) = -— — | Ao- («. y,t = a - ws) oo 7: A,N (x. y, t = ; | 


7 i=l, 


and the form of the incident shock: 


x= Utt eY(y, 1) + Oe) (6.6) 
where 
‘ on ae , 2 , ‘cee 
yH(y, 2) = - AopCtN(x = Ut, v, t, — Xo) + > A.€ -t- NCU, y, t, Xs) — 
‘ 1 Mra . i=1 
2,3 
y:Aol.(Ut, y, t, — 0) — y- D> AL(Ut, y, t, d,) 
t l 
with 


9 


t — Ax Viti— x? — y' 


L(x, y. t, \) = tan 
M—x y 


As a numerical example, the pressure and density fields have been computed for P/py = 7.3076 or (Up — U)/c = 
1.21 and M = 0.5100. The isopyenics and isobars, as well as the curved shock front, are plotted as shown in Fig. 
It is interesting to note how strongly and suddenly each isopycnic splits off from the corresponding isobar after 

Within a thin strip along OT, there is 
This may be interpreted as a slipstream 


‘. 
traversing the separation line OT between rotational and irrotational flow. 
a steep density gradient, while there is no appreciable pressure gradient. 

of finite width. This slipstream, which must be expected with every Mach reflexion, can be clearly seen in Bleak- 
ney’s’ interferograms; also, the characteristic shape of the isopyenics in the vicinity of OT has been recently* 


verified experimentally in all details. 


APPENDIX (1) 2, a and b—namely x? + y? = ct’, x < Ut—will, also 
; in Lorentz coordinates, be a conical surface—namely 
Proof That Solution (4.10) for p\(x, y, t) Vanishes at ee = ae R — , 
Infinity and Fulfills the Conditions on the Mach Circle, © + Y= ,2< 0. (See Fig. 9,aandb.) Any point 
as § in Secti A —— . : : . , : 
tated in Section (III, 4 #, 9, 7, lying on this conical surface will be influenced 
Since /"(x") = 0 for x’ < 0, the integrand of the only by sources in the §—7-plane which lie inside the 
Hea 3 nen SE ne Der. ‘ i J 
Possio integral (4.10) vanishes when Af + ‘Ss 0. Ex area R bounded by the r-axis and by the hyperbola 
essing it in other words: ‘‘Sources of intensity p;‘?”’ 4 . ‘ e 3 
. . ; Ps r? — £% — 2ir + 2 = 0, which passes through the 
are spread only over that part of the &-7-plane which is 
to the right of the line A~ + 7 = 0. (See Fig. 8, a, b, 
and c.) Now, the conical surface appearing in Figs. 


origin (€ = 0, 7 = 0). Since in that case the integra- 


tion area R does not include sources with p5(& > 0), 
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the second integral in Eq. (4.10) does not come into play, 


and 


. Rc?Ao 
p(x, ¥, f) = x 
T 


[/ f"[a(XoeE + 7)] 
L_. VU er - & - = 


The same reasoning applies for points outside the Mach 
cone, for which «* + yr — i? >Oorx? 4 gy — cF> 0. 
The corresponding hyperbola is shifted to the left. 

(a) If X\» > 1—i.e., (Up — U)/e < 1 (see Fig. 8a) 
the bounding line (L.E.), A + 7 = 0, is relatively steep 
and will not intersect the hyperbola of any point on or 


didr (4.10a) 


Str 


outside the cone, because 


dr dé L.E \ > l 
but 
dr i 
2 = l 
dé hyp Vi + yy? 
Hence, inside R the integrand is identically zero. 
Thus, p"' 0. 
(b) If Xv» < 1—ie., (Uy — U)/c > 1—the bounding 


line is less steep than in (a) and will or will not intersect 
the hyperbola depending on the location of point 
An intersection will 


%, j, 1). See Figs. 8, b and c. 


+ ¥? > Xo, which is now smaller than 1. 
The line 3/7 V 1/X? — 
points @, 9, ? which furnish intersection from those that 
do not, hits the Mach circle at the point of tangency of 
the Mach line, as indicated in Fig. 9b. For any point 
lying between C and 7, the hyperbola will be located 
How- 


occur if &/V & 
1, which separates the 


as indicated in Fig. Sb—.e., here again, p© = o. 
ever, for points lying between F and C, the hyperbola 
will be located as in Fig. Sc—i.e., now the ‘source in- 
tensity’ is different from zero in part of the integra- 
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tion area R. However, the integrand depends only on 
ps (E < 0) and not on p;(— > 0). Only the latter 
expression is related to the shock diffraction by Eq. 
(4.15); the former is determined by the boundary con- 
dition on the body surface alone, so that the disturbance 
pressure agrees with that given by the Ackeret theory 
for steady supersonic flow, and Eq. (4.10a) yields 
Rc?Ay cee ‘ 7 = 
f’' [ai + Xo % — 5-VI 


py = — 
_ do” 


av | 
For points Z, J, / inside the cone, the hyperbola will 
again be shifted in the direction of positive & and will 
enclose sources with ps"(é > 0, 0, 7)—1.e., p will then 
be influenced by the shock diffraction. 
It is evident that p“ = O for 


1+ hy-% —7-V1 —12< 0 


i.e., to the left of EC, which is tangent to the Mach cone. 

Furthermore, it is obvious that the hyperbola cor 
responding to any point in domain I, with large co 
ordinates & or ¥, will be shifted considerably in the 
direction of negative — or negative 7. In other words: 
p\?(«) = 0, as stipulated in Section (IIT, 4). 

From the preceding discussion it can be seen that in 
either case (a) or (b), if the point under consideration 
moves continuously from the outside to the inside 
across the Mach circle of Fig. 9a or 9b, the respective 
hyperbola moves continuously into the domain which 
is covered by sources of intensity py"(€ > 0). Since 
the source intensity is everywhere finite, it follows that 
p will vary continuously across the Mach circle. 
Thus, to the order e’, there is no discontinuity of pres 
sure; the reflected shock is of 0(e). 

From the uniqueness of the solution follows that re 
quiring the solution to vanish at infinity is equivalent 
with imposing the boundary condition on the Mach 


circle. 


APPENDIX (2) 


Proof of the Identity (4.12 
If ps (%, 0, 7) has a discontinuity at Z = 


Possio integral (4.10) becomes 


I ——_ Sa 
0, 5,2) = — | dt | 
Tv 0 0 


pY(F = 


Integrating by parts, one obtains 


pM = O, 7,2) = + 


ps(é) > 0, 0, 7) 


V (i —ry’?’—#— ¥ 


= (0 and a right- and left-hand derivative with respect to t, then the 


ps(—#’ < 0,0, 7) _, 
dé 


2 
me 
a) 

— 
2 S 


é g=Vii 


1 el —F - P 
y| dr Ex > 6,0, r) cos * —— = 
T 0 V il = rT)? ar x t=0+ 


(Equation continued on page 328) 
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¢/ 
s 


Vi-7r)? — - 


ra 
— 
| 
Sn, 
—~H 
o 
~ 
¢ 
> 
% 
V 
= 
x< 


bs(—t’ < 0, 0, 7) cos! 


¢ *0 Zr 5 
cos~! ; = dé — / psz°?(—é#’ <0, 0, 7) cos! = dt 
VGi-—r)? —¥ J+vi-r 7? Vidi-—7r)?-—-7, 
, 1 “t—F*F 
pM (E = 0, 7, 2) = -3f [py(x = Ot, O, 7) + p3?(x = 0-, 0, r)] dr — 
Z /J0 


l "Fr —J / (f—1r)?— 9? é 
dt p(E > 0, 0, 7) cos! <—- 
7 JO J0 Vi-7r)? -— - 


70 ¢ - 
/ p56) (—#’ < 0, 0, r)cos* : : dt’ | 
+ V (i —1)?—§? V(t{— r)? — 7 | 


Differentiating both sides with respect to y, and letting y approach zero, the terms due to the second integral ap- I 

proach zero. Therefore, one obtains Eq. (4.12), si 

| ~ 

lim ps" (x = 0, 9,0) = [ps Hs = H+. 0, {) + p3(& = 0-, 0, t) | ure 

eh | - regi 

the 

cor 

mes 

7 

an¢ 

APPENDIX (3) This leads directly to Ma 

Wit 
Derivation of the Solution of Eqs. (4.14) to (4.16) P'2(?) = Re*A2f" (at) 
rs ?) = Rc?A3f" (ai) 


The partial differential equation (4.15) for p{(# > 
0, ¥ = 0, 2) is always hyperbolic, since the discriminant The constants A» and A; are determined uniquely by 


the two simultaneous equations 





(2M)? — (4/Mo?) = 2(y + 1) (1 — WM) X " = 
; : ; Az + Az — [H(—Xo)/A(X) Ao = Qu Ay 
= " S _ , _ a 1 
(: 7 9 uw) *e hohe + A343 — ApAol[ F1(—Xo)/H(Xo)] = 1 
—Ao(X + 4M) + 4Mu : 
for all values of bo ij from 0 to l. . -_ a a en 
The solution of the differential equation (4.15) is 
i J 0 
™ A 7 — - _ — - _ es 7A 
p;“ ir > 0, y= 0, t) = [s(7 — Ao) +- I (ft — A3t) + 2 A3 ( 
“ Ft is) “wT - = =\7 
ad H(X%») * oJ [a(t — ot) | REFERENCES : 
0) 
' Bargmann, V., N.D.R.C., Report AMP 108.2R, 1945 

where \» and 2; are the roots (distinct, real and positive) ? Lighthill, M. J., Proc. Roy. Soc., Series A, Vol. 198, pp. 454 

of the quadratic equation, #70. 

*’ Courant, R., and Friedrichs, K. O., Supersonic Flow and 
H(X) = dX? — 2M + (1/M,2) = 0 Shockwaves, pp. 326, 341; Interscience Publishers, Inc., New ia 

; York, 1948. 
. . : ° ‘ Ts . i . *) 7 2 ‘on x6 th: 
T's and [; are two arbitrary functions to be determined Lamb, H., Hydrodynamics, 4th Ed., pp. 7, 8; Cambridg i 
sate. University Press, 1916. di 
by the boundary conditions (4. ‘ (4.16). ; : 
' haan ons (4.14) and 4 ») 5 See, e.g., Gardner, C. S., and Ludloff, H. F., Journal of the sic 
r(?f) + 13(2) + Re2Aol Ht =F H (Xo) 1" (ai) = Aeronautical Sciences, Vol. 17, No. 1, pp. 47-59, January, 1950 th 
‘ x 9 Ww =F 6 Ferdman, S., Master’s Thesis, New York University, College n 
Re?(2u — Ao)f" (al) of Engineering, 1951 
os ' s ™ - i ’ Bleakney, W., and Taub, A. H., Interaction of Shock Waves, th 
./ ,/ 4 2, nat Lsrrti a _. : - ° ? ’ J , 
Nols’ (2) + Asl's’(2) + ApaRe Ao [H| No) /H1(Ao) | f (al) = — Review of Modern Physics, Vol. 21, No. 4, pp. 584-605, 1948 is 
— ARc*[Ao(ro + 4M) — 4Mu]f’’’ (at) 5’ White, D. R., Doctoral Thesis, Princeton University, 1951 
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The Measurement of Downwash and 
Sidewash Behind a Rectangular Wing at 
a Mach Number of 1.6 


THEODORE DAVIS* 


United Aircraft Corporation 


ABSTRACT 


Downwash and sidewash angles were measured with a conical 
probe in vertical planes 0.5, 1.0, and 1.5 chord lengths downstream 
of a rectangular wing of biconvex section. Finely spaced meas 
urements were taken to determine the downwash field in the 
region of steepest gradients, with the result that the rolling up of 
the vortex sheet was well defined. The experimental results are 
compared with the theoretical solutions calculated from the 
measured wing loading. 

The Mach Number and dynamic pressure in the empty tunnel 
and behind the wing were also determined. The free-stream 
Mach Number was 1.60, and the Reynolds Number, based on 


wing chord, was 0.58 X 108 


List OF SYMBOLS 


b = span of wing including image; span of vortex line 
q = dynamic pressure 
Vl = Mach Number 
= surface distance along wing in x direction 
tl = stream velocity in x direction 
u = velocity perturbation in x direction 


= velocity perturbation in y direction 


= velocity perturbation in z direction 


2, 9, = coordinates of point at which downwash is wanted 

x’, y’, s’ = coordinates of point expressed in per cent semispan 
of wing; asx’ = x/(b/2) 

a = wing angle of attack 

8 = ~V/M? — 1 

r = vorticity £u ds 


INTRODUCTION 


 gnsetenagy OF THE FLOW FIELD behind the lifting 
surfaces of a supersonic airplane or missile is 
important in handling the various stability problems 
that arise. Several references are available which pre- 
dict, by means of linearized theory, the downwash and 
sidewash behind wings of various plan forms. Each of 
these theories is subject to the limitations of small 
angles of attack and a nonviscous fluid, with the result 
that the downwash region directly behind the wing tip 
is shown as a discontinuity. 

In order to provide basic experimental data on the 
subject, a program was established at the U.A.C. Re- 
search Department to measure the downwash and side 
wash angles at various stations behind a straight rec- 


tangular wing. 

Received December 12, 1950. Revised and received Novem 
ber 15, 1951 

* Aeronautical Research Engineer, Research Department 


The complete results of these investigations, which 
are too extensive for publication herein, are presented 
in reference 1. Representative results obtained from 
the test program and a discussion of the principal as- 
pects of the instrumentation problems are presented 


in this paper. 


DESCRIPTION OF TEST EQUIPMENT 


Tunnel 


The test program was conducted in the U.A.C. Re- 
search Department 8- by 9'/,-in. supersonic wind tunnel 
This 
wind tunnel is a vacuum-type, continuous-operation 
channel powered by two turboblowers driven by two 
Condensation control at this 


at a nominal free-stream Mach Number of 1.60. 


300-hp. electric motors. 
Mach Number was obtained by recirculating enough 
of the tunnel flow to raise the total temperature to 


approximately 200°F. 


Models 


The model and probe were mounted on the sidewall 
of the wind tunnel as shown in Fig. 1. The rectangular 
wing was 6 per cent thick with a 2.0-in. chord and a 
2.7-in. semispan. The section was biconvex with the 
maximum thickness located at 50 per cent of the chord. 
The angle of attack was controlled to within 0.10° by 
means of a screw-driven protractor keyed to a shaft 


attached to the model. 








Tunnel installation of wing and conical probe 
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Fic. 2. Sketch of four-orifice conical probe. 


Conical Probe 


The conical probe shown in Figs. | and 2 had a total 
apex angle of 20° and was equipped with four 0.015- 
in. pressure orifices spaced 90° apart in a plane 0.25 
in. from the apex of the cone. The leads from the 
orifices were connected to water U-tube manometers in 
diametric pairs and were also connected to a mercury 
manometer. The water manometer was used to indi- 
cate the pressure difference across the cone in the down- 
wash and sidewash planes, and the mercury manometer 
was used to measure the average cone surface pressure. 
The cone was calibrated at Mach Numbers of 1.50, 
1.60, and 1.70 over a range of yaw angles from —5 
to +10 It was found that the slope of the calibra 
tion curves was essentially constant and independent 
of Mach Number over the range of interest. 


TEST PROCEDURE AND ANALYSIS OF DATA 


A survey of the flow in the empty tunnel at all the 
probe locations was made with both conical and total 
head probes. The wing was then installed, and the 
survey was repeated for wing angles of attack of 0°, 3°, 
and 6°. Finely spaced downwash measurements were 
taken in the region directly behind the wing tip. 

The downwash and sidewash angles were determined 
from the probe readings and the probe calibration 
curves. The data were corrected for tare flow condi- 
tions by subtracting the angles for empty tunnel flow 
from those measured behind the wing. 

The Mach Number at each test point was deter- 
mined from a composite curve of the ratio of pitot 
pressure to cone surface pressure. With the Mach 
Number known, it was possible to determine the static 
pressure and the dynamic pressure. 


DISCUSSION OF RESULTS 


Conical Probe Calibration 


Since total head readings were desired for use with 
the average cone pressure in the calculation of Mach 
Number, the feasibility of incorporating a total head 
tube into the cone was investigated. It was found that, 
for the 20° cone angle with the static orifices located 
close to the apex, the presence of the total head orifice 
caused a significant change in average cone surface 
pressure. * 


* Conical probes that can be used to measure the true surface 
pressure and pitot pressure simultaneously have since been de- 


veloped. 
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The validity of the average cone pressure and pitot 
pressure measurements may be questioned when the 
readings are taken in a stream having downwash and 
sidewash components. The effect of yaw on the read. 
ing of a pitot probe was checked in the supersonic wind 
tunnel, and it was found that 10° of yaw produced 4 
change of 4 per cent in the measured pitot pressure. 
Sinular calibrations of the cone indicated that the 
maximum change in average cone pressure was 2.85 
per cent of the pressure measured with the probe aligned 
with the air stream. Since both changes at angle oj 
attack were in the negative direction, the net change 
in the ratio of the pitot and static pressures was 1.3 
per cent, which resulted in a difference of less than 0.00] 
in the calculated Mach Number. Subject to  thes¢ 
limitations, the use of the total head pressure and the 
average cone pressure without correcting for the down 
wash and sidewash angles was justified for the measure 
ment of Mach Number. 


Downwash and Sidewash 


The downwash data that best indicate the suitability 
of the test techniques were obtained in the x-y plane of 
Fig. | at z = O where the steepest gradients were en 
countered. The downwash angles measured in this 
plane at stations 0.5, 1.0, and 1.5 chord lengths down 
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Fic. 3. Downwash distribution behind a rectangular wing. 
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Sidewash distribution behind a rectangular wing. 
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Comparison of theoretical and experimental downwash 
distribution behind rectangular wing 














%:4°% M 4 y= guds (FT®/SEC 
—— _ : ——————— 
NG 50 
———— im, 
__~ a = 
ooc 
——e a ae 
™ ~ 
TIP —e 
~ 20 oc 
Wine es 7s an . 
™. 


DISTANCE FROM TEST SECTION CENTERLINE, Z-INCHES 








EE 
. 
TRAILING EDGE 
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;. 7. Change in Mach Number and dynamic pressure caused 


by wing at angle of attack. 
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stream of the trailing edge are shown in Fig. 3. In 
order to define the transition region between downwash 
and upwash, the spacing between probe settings was 
reduced from '/, in. to as little as '/3. in. Even with 
points so closely spaced, the absolute maximum values 
of the downwash and upwash angles were not estab- 
lished because there were regions less than a !/j¢ in. 
wide where the indicated angle varied as much as 7‘ 
or 8°. 

The sidewash angles downstream of the wing are 
plotted against distance above and below the plane of 
the wing at the 92 per cent span station in Fig. 4. 
The closely spaced probe settings used for downwash 
were not repeated for these tests, but it is believed that 
sidewash gradients steeper than those shown did exist. 
It should be noted that the angle scale at the 1.0 
chord station is twice that for the '/. and 1!/» chord- 
length stations. In general, the sidewash angles are in 
agreement with what would be expected from considera- 
tion of the vortex pattern. 


Theoretical Considerations 


Most of the theoretical solutions for the downwash 
are based on the assumption of an ideal uniform dis- 
tribution of lift on the wing. 
lift distribution is not in agreement with the theoretical 


Since the experimental 


loading as a consequence of nonuniform velocity dis- 
tribution in the wind tunnel or interference flow from a 
body, the experimental values of downwash often 
A method of 
calculating downwash from a known wing loading was 


differ from those predicted from theory. 


developed by Smith in reference 2, and this method was 
used for the comparison with the experimental results 
presented in this paper. Reference 2 is a further de- 
velopment of Schlichting’s work on this subject re- 
ported in reference 3. 

The wing loading was determined experimentally 
from a field of pressure taps incorporated into the wing 
This loading was then applied to the determination of 
the vorticity, as follows: On the basis of a linear 
approximation, the measured pressure coefficient, C,, 
may be expressed in terms of the free-stream velocity, 
uo, and the perturbation velocity, u, as 


C, = 2[1 — (u/u)] (1) 
Solving Eq. (1) for u, 


u = m(l1 — (C,/2)] (2) 


P 

The experimental pressure coefficients were sub- 
stituted in Eq. (2), and the values of u at various chord 
stations were calculated. These velocities were then 
plotted against the chord, s, and the various line inte- 
The 
path of integration was from the leading edge, over and 
through the wing, and back to the leading edge. The 
values of I for an angle of attack of 4° are plotted on 
Lines of constant 


grals, fu ds = I, were evaluated graphically. 


a layout of the wing in Fig. 5. 
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vorticity are drawn through points having equal values 
or I’. 
lines were curved, and it was necessary to approximate 


Since the lift distribution was not constant; the 
» the 


them with straight-line segments parallel to the leading 

- . ss) 
edge. A total of 28 horseshoe vortices were used jp 
obtaining the solution for the downwash. The down. 


wash was then calculated from the equation 


PPV x’? — p(y’ — 1)? _ 


ne Vx"? — By’ + 1)2 
mb x’(y’ — 1) 


x'(y’ + 1) 


° 
[ee 


9 
\d) 


A grid, giving the values inside the brackets for varioys 
points relative to the vortex line, was prepared for use 
with this problem. 

An important advantage of this theoretical treatment 
is that it considers the effect of the mirror image of the 
wing at the tunnel wall and thus accounts for the root 
Mach cone. The method also gives a solution based on 
a measured, rather than an assumed, wing loading, 

The theoretical downwash in the plane of the wing 
one chord length downstream from the trailing edge for 
angles of attack of 2° and 4° was calculated by the 
wing loading method and is plotted in Fig. 6. The 
lines of constant vorticity were carried straight to the 
wing tip because insufficient pressure measurements 
were available to justify fairing the data inboard of the 
tip. As a consequence of this assumed distribution of 
vorticity, the center of the vortex is shown directly 
downstream of the tip instead of slightly inboard as 
indicated by the experimental data also shown in Fig, 
6 for angles of attack of 3° and 6 The variations of 
tneoretical downwash inboard of the tip are attributed 
to variations in wing loading. The pressure measure- 
ments indicated a dip in the spanwise loading curve at 
45 per cent span which lead to an upwash in the caleu- 
lated flow distribution behind this region. Unfortu- 
nately, it was impossible to check the points less than 
an inch from the wall because of the limitations im- 
posed by the dimensions of the downwash probe. 


Mach Number and Dynamic Pressure 


The change in Mach Number and dynamic pressure 
distribution with the wing installed at 3° angle of 
attack is shown in Fig. 7, where a minus sign indicates 
It is 
readily apparent from this plot that the primary dis 
The points along 


a decrease from the empty tunnel condition. 


turbance is in the wake of the wing. 
the horizontal tunnel centerline are actually slightly 
above the wake centerline and, hence, do not show the 
maximum variations existing in the wake. 


CONCLUSIONS 


On the basis of the overall results, some conclusions 
about the nature of supersonic downwash may be drawn. 
The flow field behind a wing is greatly influenced by the 


(Continued on page 340) 
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Supersonic Wing-Body-Tail Interference’ 


GEORGE MORIKAWAt 
Hughes Aircraft Company 


ABSTRACT 


The study of wing-body-tail interference is generally compli 
cated by nonlinear effects, particularly in the description of the 
flow field behind lifting surfaces. In this paper, these difficulties 
are avoided by restricting the study to angle of attack near zero 
for in-line configurations in which both wing and tail lie nomi- 
nally in the same horizontal plane. Comparison of interference 
effects on the tail due to the two /imiting representations of the 
flow field behind lifting wings—namely, the vortex sheet and 
shows the same qualitative behavior for the limiting 


line vortex 
Thus, the indications are that 


cases of the design parameters. 
this analysis, which essentially yields the maximum interference, 
may be applied to predict the behavior of configurations in which 
the tail is not in the same plane as the wing—i.e., where the 
maximum interference occurs at nonzero angle of attack. 


NOTATION 


= radius of body 


A.R. = aspect ratio 

i = \ I 

(x) = spanwise lift distribution 
L = lift 


interference lift on tail-body for configuration and flow 


field given by (C) in Fig. 3 


m = ratio of tangent of leading edge angle (measured from 
free-stream direction ) to tangent of Mach angle 
Vo = free-stream Mach Number 


= free-stream dynamic pressure 
= lifting surface semispan measured from body center- 


line 
W = free-stream velocity 
x = spanwise coordinate 
y = vertical coordinate normal to plane of wing or tail 


= x + vy, cross-flow plane in complex variable notation 


= line vortex position in first quadrant of s-piane 


a = angle of attack 

a = angle of attack of body 

8 = V M2 l 

r = line vortex strength 

t = transformed plane (see Fig. 7 

” = efficiency 

d = a/s, diameter-span ratio 

Ne = a/xo, diameter-line vortex span ratio 

rm = ctn~! V WM? — 1, Mach angle measured from free- 
stream direction 

¢ = velocity potential in cross-flow plane 

¢ = stream function in cross-flow plane 

Subscripts 
W = pertains to wing 
i = pertains to tail 
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INTRODUCTION 


Eine. general STUDY OF supersonic wing-body-tail 
interference is inherently one involving nonlinear 
effects even for relatively small angles of attack, sav a 
< 10°. For finite lifting systems, these effects are due 
primarily, directly or indirectly, to viscosity. Con- 
sider a typical lifting system like that shown in Fig. 1, 
consisting of a slender body, a forward lifting system 
(wing), and an aft lifting system (tail). For certain 
relative positions of body nose, wing, and tail, each 
component part of the total configuration can be stud- 
ied separately, and the indications are that a satisfac- 
tory correspondence exists between experiment and 
linearized theory. However, in many cases, aerody- 
namic interference is unavoidable, and, hence, it 1s es 
sential to minimize any adverse effects. From momen- 
tum considerations, it is evident that any finite lifting 
system will induce a flow field that will affect a lifting 
system placed downstream of it. Behind the sharp 
trailing edge of a wing, but in its vicinity, experiments 
indicate that a vortex sheet well represents the flow 
field, but farther downstream a continuous rolling-up 
of the sheet generally takes place. This is a nonlinear 
process that is difficult to describe; 
attempts, with due reference to previous work, are found 
Even the description of the two- 


the most recent 


in references | and 2. 
dimensional motion of line vortices (or fully rolled-up 
vortices) in the presence of a circular cylinder is dif- 
ficult if more than two vortices, symmetrically placed 
with respect to the free-stream direction, are involved 


The present investigation avoids these difficulties by 
restricting the study to angle of attack a near zero, the 
analysis being pertinent then for in-line configurations, 
in which both wing and tail lie nominally in the same 
the validity of linearized 
Then, with the 


horizontal plane (within 
theory) passing through the body axis. 
possible exception of the viscous wake interference, sta 


tionary linearized theory is applicable. 


FORMULATION OF PROBLEM 


In order to make the aerodynamics of the configura 
tion shown in Fig. 1 mathematically tenable, yet hold- 
ing to reasonable assumptions, some essential simplifi- 
cations can be made. The principal condition that 1s 
imposed is that the configuration and flow field can be 
separated into component parts and analyzed sepa- 
rately —a general procedure that is inherent in the analy- 
ses of practically all physical problems but which must 
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always be performed with great care. The configura. 
tion to be studied is shown in Fig. 2a, where the line 
vortices schematically represent the flow disturbance 
created by the wing. Then, in many cases, the cop. 
figuration can be separated into the components showy 
in Fig. 2b. 

The separation of nose and wing-body is feasible 
only if the flow field itself can be separated in a simple 
manner—for example, if complete information is known 
in the cross-flow plane at the separation. Such is the 
case if the nose disturbance can reasonably be assumed 
to be at — ~ with respect to the wing-body, so that the 
flow field in the cross-flow plane at the separation js 
that for flow past an infinitely long cylinder. Simi- 
larly, the wing-body and tail-body can be separated, 
The aerodynamics of the wing-body without afterbody 
is given in reference 3, and some progress has been 
made in the difficult analysis of the flow field around the 
afterbody in reference 2. In the present analysis for 
a — 0, the flow disturbance created by the wing is 
represented by either an unwarped vortex sheet or 
line (fully rolled-up) vortices that remain in the plane 
of the wing. Then the wing-body-tail interference 
problem is reduced to the study of the tail-body with 
given flow discontinuities. By linearization, this prob- 
lem can be reduced further by superposition, as shown 
in Fig. 3. For the lifting problem—i.e., body angle of 
attack, ay = arp (tail angle of attack) 
3 reduces to the problem given in reference 3; the 


part (B) of Fig. 


wing-body-tail interference problem finally reduces to 
that given by part (C) of Fig. 3. 

The present study is restricted to the aerodynamics 
of a configuration consisting of a semiinfinitely long 
circular cylinder, a delta wing, and a delta tail with no 
afterbody, although part of the analysis—namely, that 
given by slender body theory—is directly applicable to 
plan formsof arbitrary-shaped leading edges and straight 
trailing edges normal to the free-stream direction. The 
analysis can easily be generalized to cylindrical bodies 
of other shapes and wing and tail plan forms with 
straight edges in a manner similar to that described in 
reference 3. 


DELTA TAIL-Bopy Lirtr DuE TO DELTA WING-Bopy 
INTERFERENCE 


(1) Parametric Representation 


Taking the same general point of view as expressed 
in reference 3, certain limiting cases of the aerodynamic 
parameters are considered, from which some insight 
for the more difficult analysis over the entire range of 
the parameters is obtained. Consider the simplified 
configuration at angle of attack ao, ay, and az, as 
shown in Fig. 4. The basic parameters for this family 
of delta wing-, delta tail-body combinations are 8 = 
ctn p = VM? — 1, dw, Ar, My, and my, where yu 1s the 
Mach angle, \ = a/s is the diameter-span ratio, m 1s 
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SUPERSONIC WING- 
TABLE | 
Limiting Cases 
Case I I] II] IV 


0 {) x 
0 0 


the ratio of the tangents of the leading edge angle to 
the Mach angle, and the subscripts correspond to wing 
and tail, respectively. Again, asin reference 3, the basic 
aerodynamic variable is taken as the ratio of specified 
lifts, so that there is no Mach Number dependence ex- 
The limiting cases are 


plicitly i.e., 8 iS eliminated. 
given by the limiting values of the parameter m = 0, 
»o, where m = O implies slender body theory and m = 


» implies strip theory. Then there are four possible 
cases given by my = O with m; = 0, © and my = 


For a more familiar geometrical 


0 


with my = O, «= 
interpretation of the parameter m, note that, for a delta 
wing alone, m = BA.R./4; BA.R. is called the effec 
tive aspect ratio. Of the four limiting cases, the two 
cases given by my = 0 with m; = 0, © are of more in- 
terest, since (1) for my = © with m, = O, the inter 
ference lift approaches zero like 1/my as my — ©, and 
2) for my = © with my; = ©, the interference lift 
approaches zero like 1/(my-my) aS My, Mr —> 
However, the last two cases are by no means trivial, 
later. 


and some remarks will be made about them 
Table 1 earmarks the four limiting cases. 


2) Representation of Interference Lift 

The interference lift 1, can be presented in several 
different forms (lift ratios), each of which is useful in 
the various considerations necessary in the design 
analysis of efficient stable configurations. The lifts 
pertinent to a useful representation are easily visualized 
by the superposition scheme in Fig. 5, where the body 
is eliminated to simplify the discussion; Ly, Lr, and 
L, generally include body lift. Then the interference 
lift L; may be represented by the following lift ratios: 
(—L,/L, (1 — nr) (la) 
where 7; is defined as the tail efficiency; 
( —L, Lw) (1 Nw) (1b) 

where ny- is defined as the wing efficiency ; 


-L,/(Lw + Lr) (1 — ») (le) 


where 7 is defined as the total lifting surface efficiency; 


( L; ) (1 -— Nw) (1 — Yr) 
= (id) 
a L, Ly 


(ld) is usually made di 


and 


The lift ratio given by Eq. 
mensionless by expressing either L; or Ly as a lift coef 
ficient based on the body cross-sectional area.* It is 

lo the author’s knowledge, A. E. Puckett and I. Naiman first 


introduced this ratio and have used it effectively in the design of 


high-speed configurations 
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B) EQUIVALENT PROBLEM 
Case I with vortex sheet and equivalent problem for 
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A) CROSS- FLOW PLANE 


B) TRANSFORMED PLANE 


Fic. 7. Cross-flow plane for Case I with line vortices 
TABLE 2 
Identically Zero Lift Ratios 

Case I II III I\ 
(1 — n7 0 ) 
(1 — oF 0 0 
(l—7 0 0 0 
(1 — nu 

Ly 


noted also that (1 — ny) given by Eq. (1b) is the only 
i.e.,0 < (1 — gw) S 1 
In addition, 


true efficiency ratio this is 
evident from momentum coiisiderations. 
for the line vortex representation of the wing disturb- 
ance, (1 — my) [and, hence, (1 — nw)/Lr] can be ex 
pressed independent of the vortex strength 
pendent of the wing loading and geometry. 


i.e., inde- 
Thus, 
these lift ratios afford a somewhat more general inter 
pretation of interference lift. 

It is of interest to note from Table 2 that certain lift 
ratios are identically zero (for all Aw and A;) for some 
of the four limiting cases considered. 

Table 2 clearly shows that limiting Case I gives the 


maximum nonzero information for all lift ratios and that 
(1 — nw)/L7 (1 — nr)/Ly 


vields the maximum nonzero information in the com 
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parison of the limiting cases. Most of the present 
analysis is presented in terms of (1 — nr) and (1 — 


nw), both having simple physical interpretations. 


(3) Case I (mw = mr = 0) with Vortex Sheet 


Unless otherwise stated, the lifting problem with ao 
= ay = ar is considered. Referring to Fig. 6, this 
limiting case, assuming an unwarped vortex sheet, has 
a particularly simple solution (slender body theory), 
and comparison with other limiting cases (in particular, 
Case I with line vortices) is especially instructive. For 
Ar < Ap (tail span greater than wing span), the dif- 
ference between the total lift of the equivalent prob- 
lem, and the total lift without interference is clearly 
the lift on the cross-hatched area in Fig. 6b. Then the 
interference lift is exactly the wing lift Ly. Similarly 
for \y > Aw, the interference lift is the tail lift Ly. 


AERONAUTICAL 


SCIENCES MAY, 
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(4) Case I with Line Vortices 


Since slender-body theory is applicable for this case, 
the analysis is reduced to a two-dimensional potential] 
problem in the cross-flow plane (cf. reference 4). They 
the problem for interference lift given by Part (C) of 
Fig. 3 becomes that shown in Fig. 7a, where the usual 
complex variable notation is used—i.e., 2 = x + jy, 
The vortex strength is denoted by I (positive counter. 
clockwise), and the boundary condition is that the 
velocity ¢, normal to the surface is zero on body and 
tail. Then, by symmetry about the y-axis, the problem 
for the right half-plane is analogous to that of finding a 
Green’s function. This is most simply accomplished 
by transforming the tail-body into a circle, as shown in 
Fig. 7b, and using the known reflection properties on 
the circle in order to satisfy the boundary conditions, 
The transformation z — ¢ is easily seen to be 


Thus, 
a” | a*\ (. | 
Ly —lLy - l = a) & . - a8 9 5 a $ + 3) 
= — = —2za’* | ——— . - . s 
qa qa Aw 
(for ke <= Aw) and the solution for the complex potential in the ¢- 
ee (2a) plane is 
4] —Lr , 1 — Ar” “ 
= = —2ra’" T 4 Co)lC + (1/e)] 
qa qa Ar g+w= - ; In +. a = (4a) 
(for Ar > Aw) ) 2a [¢ =a fo) }(¢ = to) $ 
and Writing A = [s + (a*/s)] for convenience and trans- 
forming back to the cross-flowplane, the argument of 
Ar 7h = hy" the log becomes 
(1 — gr) = 2 
a ad ot 1\ 
(for Ar < Aw) a) Cc = fo) (« + J (: -_ ‘) - (« is *) 
$o $ $o 
= | (for \r > Aw) | - (4b) 
_ _ , ‘ ‘ , = ‘ (: —_ ala -+- £0) (: sd ‘) + (< = *) 
rhe other lift ratios are obtained simply [ef. Eq. (1) ]. Fo § 60 
Then on the tail or body ( = x,a <x < s; ors = ae’), the denominator of the argument is the complex con- 
jugate of the numerator, and the potential on the tail is 
/ 
Im V1 — (1/A®)[20 + (a?/z0) |? | 
(5a 


r 
¢o(x, 0; 2) = — tan7! 1 a? | a?\2} 
7 i=> (» +-+R@i — =, (« + ) 
A? x A? Z 


where R/ and /m denote the real and imaginary parts of the radical expression. If the vortices are in the plane 


of the wing (2) = Xo), then for x9 < s the imaginary part is identically zero and for xo > s the real part is zero. 


Thus, in the \ notation, * 


g(X, Ar; Ao) = O (for Ar < “ 
r \ _ (Ar? — Ao?2)(1 — Ao2Ar?) ; (5b) 
= —tan™! - = xe (for Ar > Xo) 
T No (A? = Ar*)(1 = do7A7) 
Similarly, on the body (z = a”), 
G(X, Ar; Ao) = 0 (for \7 = ”| 
r 1 /(Ar? — Xo2)(1 — Ao?Ar?) , (6) 
= — tan =n ceraoaee (for Ar > Xo) 
+ No (1 + Ar*)? — Ar*(x7/a’") 


* Note that here and in the integration [cf. Eq. (8a)]\ denotes the running variable a/x. 
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It should be noted that Eqs. (5b) and (6) state that the velocity potential ¢ is identically zero downstream of the 


cross-he atched area shown in Fig. 8 
the line vortices lie in the plane of the tail.* 
cf, Fig. 6b) shows an analogous behavior. 


potential, since 


== 2 | | dx = = 
q a q a! 


}and thereby imply that the lift is nonzero only on the cross-hatched area when 
Also a comparison with the results of Case I with the vortex sheet 
The interference lift can be expressed as an integral of the velocity 


LU + 


8 *s 
ds) dx = =a [y(t.e.) — o(le.)] dx 


where (1.e.) and (t.e.) denote the leading and trailing edges, respectively, and W is the free-stream velocity. 
Then the wing portion of the interference lift, combining Eqs. (5b) and (7), is 


L /g = —4(a-T/W)[(1 — do) /Ao] (for Ar < Xo) 
aah S/a-T\ [' do _ | (A — Av?)(1 — Ag?A2) [A 
= : tan! = = — a (Sa) 
q r\ M7 (Ar? — Ao”)(1 — Ao?Ar*) | A? 
(for Xr > Xo) J 
and combining Eqs. (6) and (7), the body portion of the interference lift is 
L;| 8 fa: ") "e | | dx. ) 
= Vf tan a (for Xr < Xo) 
q W 4 re 1 — (x*/a*) | a 
wae “(a “\ (Ar? — Yo®)(L — do?Ar?) | t (Sb) 
W 0 No (1 + Ar”)? — 4Ar?(x?/a?) 
— do’ l } dx : z 
tan eer | (for Xr > Xo) 
2ro 1 — (x*/a*) \ a | 
Integrating by parts, Eq. (Sa) can be reduced to standard elliptic integrals. 
Both integrals of Eq. pe are of the form 
1 : 
[= i tan—' (B/V C* — x*) dx 
0 
where B and C are constants with respect to the integration. Integrating by parts, 
/ / 
B | VB?+C 2BV (B? + C*)(C? — 1) 
I = tan™ ; + Btan™ = _ tan= = ——pr = (Sc) 
7c? — 1 wot = 1% 2 (B? + C*)(C? — 2) + C 
The simpler of the two integrals on the right side of Eq. (8b) has a particularly simple solution, 
a i 1 — Xo? | 
tan! - |} dx = (1 — Xo) (8d) 
us 0 2ro I — 
The total interference lift, adding Eqs. (Sa) and (Sb), is 
Li g = (L|}™ + L;) qd (Se) 





* This is easily visualized by considering a two-dimensional line vortex near a plane wall and evaluating the velocity potential in 


the limit as the distance between vortex and wall approaches zero. 
elsewhere it is the usual complex variable. 


+t Note that the free-stream direction is here denoted by :; 


The interference lift is in relatively general form, since 
it is expressed in terms of vortex strength I and posi- 
tion Ay. For the particular case of a delta wing with 
My = 0, the relationship between (1) wing loading and 
geometry and (2) I’ and Xo is necessary. Recalling the 
assumption (for a — 0) that the rolling-up process of 


the vortex sheet behind the wing is such that the fully 
rolled-up vortices (line vortices) remain in the plane of 
wing and tail, it is easily seen that the lateral center of 
gravity of the vortex sheet on either side of the body re- 
mains invariant in the downstream direction. This 
implies that each line vortex on either side of the body 

































































338 JOURNAL OF THE AERONAUTICAL SCIENCES MAY, 1952 
. : . is exactly at the center of gravity position of the Vortex 5 
sheet. Thus, a relationship is found between the line i 
r° vortex position x» and the spanwise vorticity distriby. by 
- FA r tion at the trailing edge of the wing, and the vortey nd 
MMi a0 i strength is directly related to the spanwise lift evaly ing 
4 | x : ; , Y > ees . 5 
r : 6 4 a ° ated at the root chord of the wing (Kutta Joukowsk 
q b 4 law). 
LLLI_ Ly = ine. 
— f 4 rhe reflected vortex inside the body is at the InVerse 
ceaemeactiace ee i point a*/x 9. Then calculating the body portion of th 
4 4 | é : the 
“<Q total wing lift corresponding to this line vortex spacing 
te | indicates that a lifting-line representation of the Wing 
“F | body, as shown in Fig. 9, yields a good approximatio; 
to the total lift and also gives a reasonable approxima 
oe — aoe tion to the wing disturbance some distance behind thy 
Fic. 8. Case I with line vortices in plane of tail. cs die ne 
trailing edge satisfying the body boundary conditions 
Since the calculations are somewhat simpler, the vor Fo! 
a tex spacing given by the lifting-line representation js giv 
: used. Then the total lift is given by Ly 
* r 
Wh Ko L/q = [2al(a) Git (1 — Xo*)/Ao}t , 
r » N 
amie este =5 T => = (4al'/W)[(1 — Xo?)/Ao] 
' a | ' 
: t | . 
rp , 4 since, by the Kutta-Joukowski law, As 
/ 
ii > : my 
l(a)/q = 2T/U 9b 
Since the spanwise lift at the root chord of the wing for lin 
“t a o is 
Fic. 9. Lifting line representation of wing-body for mw = 0 my = 01s br 
l(a)/qa = 4a[(1 — dw?) /An (¢ M1 
: ‘ 7 — ne : the 
[ and the total lift for my = 0 is \ 
he22| Ao24] 2-26) Aloe Mn 
1.6 i the 
it 9 2 2 4 ( 
= 2ra°|(1 — Aw*)/Aw] J an 
Ga 
yf My 
1.4 + | ; : . ; 
the relationship between wing span and vortex spacing 
| . , 6 
is given by 
| b 
| 
1.2 } } + + 4 | 
(1 — An =) An 1 (1 — \o”) No (He 
T t re] 
| Te 
1.0} — — a Since the wing lift Ly and the tail lift (without inter- | 
ference) Ly for mp = 0 is given by Eq. (9d) (using the II 
(I-9_) ; js ‘ fj eee ; ; ' ‘ 
%T proper subscript for the tail), the lift ratios given by for 
8 t Eq. (1) can be expressed in terms of Ao [or Aw by Eq tit 
(9e)] and Ay. A comparison of (1 — 7) for Case I with ar 
vortex sheet and Case I with line vortices is shown in Wi 
6 t t 4 Fig. 10. A marked similarity between the two repre ou 
sentations is shown. In fact, the two cases are identical C; 
for X\y < Aw. The difference is found for A; > Aw at 
.4 t t t / (tail span < wing span), where for the vortex sheet se 
| representation the tail is completely ineffective —1e. tv 
My =0 (1 — nr) = 1-—~and for the line vortex representation | I. 
2 } eye the maximum interference* (more than 100 per cent 
LINE VORTEX ne > ' ; 
——. ineffective) in the interior region, 0 <A; < 1, occursat | 
VORTEX SHEET Av = Ao and has a minimum (1 — nr) of about 0.5 
oO . (about 20 per cent effective). 
0 = - 6 8 1.0 in 
rj ; = tbody ar 
Fic. 10. Comparison of (1 — nz) for Case I with vortex sheet Actually, the maxium occurs at the limit Ar = 1 (all body), wr 


and Case I with line vortices. where (1 — nr) is infinite, 
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(5) Case II (mw = 0, mr = &) with Line Vortex 


For mr = ©, the interference lift is easily obtained 
bv a “strip theory” calculation over the tail plan form 
and the adjacent body carries no lift. Then, integrat- 


ing over the tail area S7, 
5 


i irr 
i = wT / a, dS (10a) 
q S7 
L; 4a? (*") fl Cr = Ao)(L + dodr) 
= : n 
g mmr \W/ irr |(Ar + Xo)(1 — Addr) 


For the calculation of the lift ratios, the wing lift Ly is 
given by Eq. (9d) and the tail lift with no interference 
l 1 la 


é ; 1 — A,’ 
= = + 2 In Ar 
qa Myr Ar” 


As pointed out in Table 2, the nonzero lift ratios as 
mp > © are (1 — nr) and (1 — nyw)/Lr for this case. 
A comparison of (1 — nr) for this case and Case I with 
A marked similarity 


Ly by reference 3, 


(11) 


line vortices is shown in Fig. 11. 
is seen between these two cases, which describe the 
bracketing limits of the tail parameter m;—namely, 
The essential differences are (1) 


mr = Oand mp = &~. 
© occurs for 


interference for mp = 
i.e., when the tail span is smaller than 


the maximum 
Aw < Ar . do 
the wing span, but greater than the vortex span 


and (2) (1 — nr) = 8/r?, a finite limit at 47 = 1 for 
My = oe . 
(6) Case II (mw = ©, mr = 0) and Case IV (mw = 

mr = @) 


If the wing-body disturbance for these two cases 1s 
represented by line vortices, the nonzero lift ratios (cf. 
Table 2) are identical to the results for Cases I and II 

i.e., the lift ratios (1 — ny) and (1 — nyw)/L7 for Case 
III are identical to those for Case I, and (1 — ny)/Lr 
for Case IV is the same as that for Case II. These iden- 
tities are immediately evident, since the two lift ratios 
are independent of my. The representation of the 
wing disturbance by a vortex sheet has not been carried 
out for these cases, but the comparison of results for 
Case I with vortex sheet and with line vortices is indic- 
ative of the expected difference between the two repre- 
sentations. Fig. 12 shows the similarity between the 
two representations of the lift ratio (1 — nw) for Case 


I. 


APPLICATIONS AND FURTHER STUDIES 


From the few experimental results available, the 
indications are that for a — 0 the interference effects 
are well described at least qualitatively, and the viscous 
wake and other nonlinear effects modify the description 
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where a, is equal to'the inverse tangent of the ratio of 
the upwash due to the line vortices in the plane of the 


tail to the free-stream velocity IV’; for small angles, 


(i) Xo ( l a” ) 
a; = - _ — —— (10b) 
a mw \x*? — Xo" Xo?x? — at 


and 
(Ar? — do?) (1 — Ao*Ar?)) 

— — Ip |— -| Xo In | (10c) 
ho =|Ar2(1 — Xo)? Ar?2(1 — do?) f 


in a quantitative way. It appears also that, once the 
geometry of the flow disturbance behind a wing-body 
(wing not necessarily in the horizontal plane) is suf- 
ficiently well understood, then the present analysis 
can be applied to tail-body configurations where the 
tail is not in the horizontal plane. Another important 
problem is that of obtaining the dependence (non- 
linear) of the interference effect on the angle of attack 
a. Some progress in this direction has been made® 
with planar wing-tail configurations having no body. 
Further work beyond the present analysis for inter- 
mediate values of the parameter m, does not appear 
warranted at the present time, until further experi- 


mental work indicates such a need. A reasonable ap- 





T = 


Ao=8 
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(I-9y) 

















Comparison of (1 — nr) for Case I and Case II with 
line vortices 
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do?.6 ho? 8 











a = : 
ce) = 4 6 8 1.0 
AT 
My=O my =O LINE VORTEX VORTEX SHEET -— — 
Fic. 12. Comparison of (1 nw) for Case I with vortex sheet 


and Case I with line vortices 


proximate analysis for immediate values of my js es 


sential and appears possible. 
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(Continued from page 332) 


vortex originating at the tip of the wing at angle of 
attack. The core of this vortex trails out behind the 
tip and bends slightly inward and downward. The 
resulting effect on the flow behind the wing is strongest 
directly behind the wing near the vortex core, and the 
disturbances show only a slight reduction in intensity 
as much as 1.5 chord lengths downstream of the wing. 
It has been demonstrated that these effects can be 


measured by the use of a conical probe having four 


static pressure orifices and a total head orifice in the 
tip. 
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Heat Transfer to Constant-Property Laminar 

Boundary-Layer Flows with Power-Function 

Free-Stream Velocity and Wall-Temperature 
Variation’ 


SOLOMON LEVY? 
Unwersity of California at Berkeley 


ABSTRACT 


Numerical computations have been performed for the bound 
ary-layer form of the energy equation for incompressible flows 
with power-function variation of free-stream velocity (#4 = 
cx™) and of wall temperature (7) = Ax”), the pertinent solutions 
of the momentum equation in this case being those of Hartree 
The numerical computations given herein are to some extent a 
repetition of those given by Schuh and by Chapman and Rubesin, 
the object of the present computations being the resolution of dis 
crepancies appearing in the previous solutions and an extension of 
their range. IBM machine calculations were employed in the 
finite difference calculation presently utilized, the results thereof 
covering a range of wall-temperature function exponents y from 

2.5 to 4.0 and of Prandtl Numbers from 0.70 to 20.0 for four 
values of m(4, 1, 0, —0.0904). The accuracy of the numerical 
computations is examined in detail, and the accuracy of the com 
puted functions at the wall, which determine the heat-transfer 
rate, is estimated to be within 2 per cent 

Examination of the results reveals that the results of Schuh 
for the flat plate are in error. For the range of the calculations, 
it was found that the local heat-transfer coefficient can, with the 
exception of large negative y values, be expressed within +5 
per cent as 

(hx/k) 


V ux/v = Bim, y)o7 


where the exponent A of the Prandtl Number @ varies from 0.254 


to 0.367 for 0.0904 <m <4 and where the function B(m, y) 


can be approximated by the equation 


0 104 
} 0.205) 


NOMENCLATURE 


1, C, D = constant of proportionality unless otherwise de 
fined 
B(m, y) = function of m and y 
Gn, b, = coefficients defined by Eq. (14) 
Cn, d, = coefficients defined by Eq. (17) 


specific heat at constant pressure, B.t.u./lb. °F. 
step interval in the finite difference solution 
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si = extrapolation error defined by Eq. (24) 

EE, = extrapolation error defined by Eq. (25) 

f(n) solution to Eq. (7) (Hartree profiles 

k thermal conductivity, B.t.u. ft/hr. ft.? °F 

m = the exponent in the statement ™ Cx™ 

\ total number of steps in the finite difference solu 
tion 

p = the pressure in the boundary layer, lbs. per sq.ft 

g = the local heat rate, B.t.u./hr. ft 

O = the total heat rate, B.t.u./hr. ft.? 

R = extrapolation error defined by Eq. (26 

r = the exponent controlling the heat rate at the sur- 
face 

S = rounding-off error defined by Eq. (28 

7 = temperature in boundary layer, °F 

T; temperature at the edge of the boundary layer, 

F 

Ti = temperature of the surface, °F 

qs = Rounding-off error defined by Eq. (27 

u = velocity component of flow in x direction, ft. pet 
SCC 

u velocity at the edge of the boundary layer, ft. per 
sec 

1 = velocity perpendicular to body surface, ft. per sec 

x = distance along the body surface, ft 

y = distance normal to the body surface, ft 

B = velocity distribution parameter, 2m/(m + | 

Y = the exponent controlling the temperature distribu 
tion along the surface (7) = Ax? 

n independent variable occurring in the Hartree pro 
files 

4 the dimensionless temperature, 7/7 

d = the exponent determining the effect of the Prandtl 
Number 

bu = absolute viscosity, lb. sec. /ft.? 

v = the kinematic viscosity, ft.?/sec 

p = density of the flowing fluid, lb. sec.?/ft.‘ 

a = the Prandtl Number, uc,/k dimensionless 

y¥ = the stream function 


INTRODUCTION 


p | ‘HE CALCULATION OF THE HEAT TRANSFER to the 
laminar boundary layer is of importance in many 


engineering applications. This type of flow occurs in 
the aerodynamic heating of bodies in flight, in the cool 
ing of gas-turbine blades, along airfoil surfaces, and 
along rocket motor nozzles. In such instances, the 
character of the variation of the wall temperature often 


has an important influence upon the rate of heat trans 
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fer, and, in consequence, solutions of the momentum | 
and energy equation leading to an accounting of the 
effect of variable wall temperature are of significance. 
In certain circumstances, such as the stagnation point I 
region, the solutions can be used directly for design 
calculations; for others in which an arbitrary variation 
of free-stream velocity and wall temperature or heat 
rate occur, approximate methods based on these solu- B 
tions have been proposed. ! con 
All boundary-layer solutions in which variation of sim 
the wall temperature can be accounted for deal with E 
incompressible constant property flows with free-stream ert® 
velocity a power function of the distance from the stag- and 
nation point. These are the ‘“‘wedge flows” (mu = cx’) Fic. 1. Flow boundary layer about an arbitrary body Wh 
for which the solution of the momentum equation has nec 
been given by Hartree.’ Exact solutions of the corre- where the y-axis is taken perpendicular to the surface fact 
sponding energy equation are only possible for the case (Fig. 1). tior 
in which the wall temperature is a power function of a Bie bt , : ; For 
E ; : tt rhe dissipation term is omitted from the energy 
distance from the stagnation point (7) = Ax’). : . * ; : ‘ Ps. 
: ; ; ‘ ' : : equation, Eq. (3). For use with flows in which fric 
Solutions of the energy equation as given in this re hk ; a ae : eae, ee 
. é . tional heating is significant, the solution of Eq. (3 
port have already been presented by Schuh*® and by . | 
: : ; : Sg “is to be added to the appropriate adiabatic wall solution k 
Chapman and Rubesin.‘ Schuh considered the flat- ps ae te ie 
to satisfy the prevailing boundary conditions. The gral 
plate flow and two wedge flows (m = 1/2, m = 1) and “ . hae 
; z : ? temperature 7 appearing in Eq. (3) represents a tem- 
obtained his results by the method of successive ap- ; : 
: ; : - ; : perature difference between the local and free-stream 
proximation. Chapman and Rubesin considered flat- , 
ae : : temperature and, in consequence, partially accounts 
plate flow only but considered variation of properties , 
: ; : for the term, u(dp/dx), which appears in that equation 
to the extent that up = const., under which circum- ; ; 
: ; ; : when the state temperature is the dependent variable. 
stances the solutions apply for the imcompressible case, pa is 
; - ae # panies The boundary conditions are 
as well as for the case of compressible flow. This solu- : 
‘ ye . ‘ Py » y “yay 7 > > vr « ‘ > ¢ ; ‘T° “7 IS 5! 
tion was obtained by use of the Runge-Kutta integration T = T(x) go y = 0 ‘ 
‘ sf P ° = . ‘ : 0 
method. Certain discrepancies exist within Schuh’s re- Fas, OO u= ee , 
ra : fut 
sults and between the resultsof Schuh and Chapman and with uw, = cx” 
: a“ ‘ , the 
Rubesin. The present solution was undertaken to re- : 
P a P . n m . > ‘ 2e Variz ac det 
solve these discrepancies; to investigate further the ef- Introduction of the Hartree variables, 
. a . . ; at i 
fect of the Prandtl Number fo S case, since the pre- / ; 
ect f the Pra itl Nu n yer for thi case, since the pre n = Vil + m)/2)(u1/vx) y rl 
viously available solutions considered only Prandtl _ 
Numbers of about 0.70; and to present solutions for y= Vv [2 ‘(1 + m)|vxu, f(n) (5 Q = 
two additional “wedge flows” (m = 4 and m = : : . 
where y is the stream function defined by 
— 0.0904). : 
u = Oy/oy v = —OoOy/Oox (6 
ANALYSIS , , ae . , 
yields the following form for Eqs. (2) and (3) through 4 
- % 01 
Derivation of Equations Eq. (1): 7 
. . . tlk “rr 719) (> 
In the laminar boundary layer of a body in flow at J +H t+ PU —s*) = 0 \/ the 
high Reynolds Number, the velocity increases from zero Lor or 9 oT Ire 
at the wall to the free-stream velocity u;.. The bound- yn? +f > ak + 6 > = 0 (8 I 
. ° ° o e m x 
ary-layer forms of the continuity, momentum, and ” ” soli 
energy equations for an incompressible constant prop- where 8 = 2m/(m + 1) and o = uC,/k is the Prandtl rie 
erty two-dimensional flow are Number. tha 
Eq. (7) has been solved by Hartree for several values me 
(Ou/dx) + (dv/dyv) = 0 | ee ak 
of 8, and these solutions are employed for the solution not 
Ou n Ou Op O*u of the energy equation. This equation is changed to an zer 
u v = 2 a ; 
, Ox ” Oy Ox a ov? (<) ordinary differential equation by introducing the tem . 
r r r perature function 6 = 7°/7), where @ = @(n) and To = * 
. 3) 3) 0” Pet ae 
pC, (« +y = k — (3) Ax’. Thus, 
Ox Oy Oy? a 
* The primes indicate the order of differentiation with respect a 
where x is measured along the surface of the body and to 7. Ch 
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HEAT TRANSFER TO 
d°0 e dé 7 
- + Je — o(2 — B)rf'6 = 0 (9) 
dn dn 


The boundary conditions are 


dé= 1 n 
6=0 

Before the general method of solution of Eq. (9) is 
considered, a few particular cases that yield rather 
simple numerical solutions are discussed. 

Eq. (9) has been solved by Pohlhausen® and by Eck- 
ert® for the case of an isothermal wvall, so that y = 0 
and the solution can be obtained by quadratures. 
When y is not zero, numerical solution of Eq. (9) is 
necessary except for the particular case in which the 
factor (2 — B)yisequal to —1. In this case the equa- 
tion is exact and solution by quadratures is possible. 
For this special case the solution is 

d=e (10) 

For this case (d0/dn)n» = 0, and the temperature 


gradient at the wall, which is 


(°°) be (*) l 
= Ax’ 

Ov y=0 dn m= V2 — B vx 
(07) i (“*) x! 
Ov /y—0 ; dn/ no V (2 — £B) 


is seen to be zero for all values except the indeterminate 
pont atx = 0. At that point the temperature of the 
fluid immediately adjacent to the wall is infinite, and 
the heat introduced impulsively at this point can be 
determined from the heat flow in the boundary layer 


atany point x. This is 


ie eal "oT @ 
V0 = pC,mTo a dy == 
J0 T 9 uy , 
ApC, V vC(2 — 8B) / Of’ dy (11) 
0 


It can be shown that for y > — 1/(2 — £8) the heat 
flow at the wall is always into the boundary layer, while 
for y < — 1/2 — @ the heat flow is into the wall despite 
the fact that the wall temperature always exceeds the 
free-stream temperature. 

For values of y other than noted above, numerical 
solution of Eq. (9) is necessary, and this has been car- 
ried out by the method of differences. Reduction of 
that equation to finite difference form, with the incre- 
ments in @ beitiy distant e apart and subscript m de- 
noting the number of the increment beginning with 
zero at the plate, produces the equation 
a, — 20,1 + 6,2 . 6, — O,-1 

; T+ Fn = 
e e 
yo(2 — B)fraO,1 = 0 (12) 


This yields, for the temperature at increment n, 


BOUNDARY-LAYER 


FLOWS 3438 


2+ coef, + oye(2 — B)fn-1'’ 
6, = . 0, l 
Ll + o@ Ja 
I 
6, (13) 
l 4 oe} » 1 
If the factors of Eq. (13) are defined as 
2 + oe l + oye?(2 B)f,, a 


G,.4 * " 
l + oe} » 1 (1 4) 
6... = 1/(1 + osf,_s) 


n 
then Eq. (13) acquires the form 


= An—19n-1 — Dn—1On-2 (15) 


6 


n 
In virtue of the boundary conditions at the surface 
a = 1, 


6. = a0, — b, 
0, = Abo — bob, etc. 


so that all 6, values are expressible in terms of 6;. This 
may be expressed generally as 
0, = Ce-101 — dn. (n > 2) (16) 
where 
Cu-1 = Gn—iCn—-2 — DgrCp 
d, 1 = Gd, id, hs b,, id», 
where 
a= gy d, = by - 
(17) 


oo = 1 do = 0 


The calculations leading to 6, for large m were accom- 
plished on an electronic IBM computing machine, in 
which six decimals were carried in the machine and the 
results were rounded off to five decimals. The incre- 
ment e = An was taken in most instances as 0.1. In 
virtue of the boundary conditions, @, becomes zero for 
large n (i.e.,7 — ©). In this calculation, if the value 
of 6; was unaffected within the sixth decimal place by 
evaluation from Eq. (16) with @, = 0 for two positions 
nand n — 1, the number of increments was deemed suf- 
ficiently great. With this value of 6, the entire tem- 
perature field could be determined from Eq. (16). 

The temperature gradient at the wall gives the heat 
flux and as such is the most important result. This is 


(d6/dn),-0 = —(1 — 4)/e (18) 


and can be evaluated from the known value 4. 
Appendix (A) contains a discussion of the accuracy 
of the method and the probable precision of the results 


RESULTS 


Calculations have been carried out for a sufficient 
number of cases to describe the combined effect on the 
heat transfer of free-stream velocity and wall-tempera- 
ture variation and of changes in the Prandtl Number 
of the fluid. These cases are indicated in the following 


tabulation: 
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TABLE 1 
Temperature Gradients at the Wall for Flat-Plate Flow (8 = 0). Variable Wall Temperature 
= - Calculated——— — a Schuh Rubesin. 
y/o 0.70 2.0 10 20 0.70 0.715 
—0.50 0 0 0 0 
—0.25 er — 0.3789 — 0.6257 — 0.7668 cet we 
0.00 —0.4065 —0.5822 —0.9863 — 1.230 —0.414 —0.418 
0.25 —0.4989 —0.7130 —1.210 —1.513 Mae 
0.50 — 0.5690 —0.8112 —1.377 —1,721 - 
1.00 — 0.6746 —0.9593 —1.625 —2.024 —0.576 —().691 
2.00 —0.8218 —1.165 —1.965 —2.445 ne —().845 
3.00 —0.9296 —1.316 —2.211 —2.741 aa —0.967 
4.00 —1.017 — 1.437 —2.406 —2.974 —0.828 —1.053 
(a 
TABLE 2 48 = 1.6 di 
Temperature Gradients at the Wall for Stagnation Flow (8 = <y = —2.5 to 4 
0 ‘ariable Wall Te Te > - - 
: 1.0). Variable Wall Temperature le = 0.70, 1, 5, 10 
- Calculated ~ Schuh ? ; 
y/o 0.7 1 5 10 0.7 48 = —0.199 (separation flow) 
—1.0 0 0 0 0 3 dy = —0.455 to 4 
—0.75 —0.1755 —0.2001 ae —0.4062 } Ps =. 
—0.50 ne re me —0.8094 —0.318 ¢ = 0.70, 1, 5, 10 
—0.25 —0.4093 —0.4708 ie —1.081 - ; _ 
0 —0.4879 —0.5603 —-1.011 —1.286 —0.496 The large number of cases considered prohibits com- 
0.25 —0.5535 —0.6345 -1.141 —1.451 seicall il aie a Wiis. cic 
0.50 —0.6094 —0.6979 —1.251 —1.590 plete presentation of the results for all. The most im- 
1.00 —0.7033 —0.8116 —1.482 —-1.818 —0.707 portant result, the gradient of the temperature func- 
2.00 —0.846 —0.9647 —2.15 —0.85 . . > - oP 
2’ 40 61 “ i795 190 ).894 tion at the wall, is presented for all cases in Tables 1, 2, 
3.00 —0.9567 —1.089 —1.914 —2.418 : 3,and 4. The knowledge of this magnitude enables the 
00 —1.048 —1.192 —2.086 —2.630  —1.050 : 
. 1.04 1. tua ant 1.08 calculation of the heat rate at the wall as 
‘J dé I uy FIG 
TABLE 3 go = —kTo (| ) V9 8 (19 flow 
ae . : y , 7 , C N= meas vx 
Temperature Gradient at the Wall for Wedge Flow (8 = 1.6) ~ 
. icant cctnanate if the wall temperature 7) = Ax’ is given. Should the - 
- : é e 
aanme , a s ie wall heat rate, gq) = Dx’, be given instead, the wall 
y/o ‘ o — af ‘ mo 
. re . te > rE ( , 
~2 50 0.0000 0.0000 0 0000 0.0000 temperature 7) is determined by Eq. (19), with ia 
—1.50 — 0.2687 —0.3101 —(). 5587 —0.7005 , te . , 
—0.50 —0.4413 -0.5085 —0.9303 —1.186 y=r-+[(1l — 8)/(2 — B)] in 
0.00 —0.5062 —0.5828 — 1.064 — 1.357 Pave : . see oes the 
0.25 —0.5353  —-0.6161  —1.122 — 1.432 Distribution of the temperature within the boundary : 
0.50 —(0. 5626 —0.6468 —1.176 —1.501 = 2 a eeey om ee a. eee : 
1.00 -—0.6120 —0.7031 1.275 —1.626 layer is illustrated by Figs. 4, 5, and 6, covering th pre 
2.00 —0.6975 -—0.7995 —1.442 — 1.836 following cases: 
3.00 —0.7692 —0.8808 1.581 —2.010 exe 
4.00 -—0.8315  -—0.9512  —1.701 —2.159 fo = 0.70 re} 
18 = 1.6, 1.0, 0, —0.199 ; oof 
ab! 
TABLE 4 i) B= 1 
Temperature Gradients at the Wall for Separation Flow (8 = lc = L, & 10 
—0.199). Variable Wall Temperature j 
f re iv 7 te ~ = —_ 2- vere 
. a All values given for the cases 1/(2 — 8) were 
y/o 0.7 1.0 5.0 10.0 obtained from Eq. (10) rather than from the fimte [| 
—1/2.199 0.0000 0.0000 0.0000 0.0000 difference calculation. 
—0.25 —0.1955 —0.2168 —0.3290 —0.3894 
0.00 —0.2930  —0.3227 —0.4884 —0.5806 
0.25 —0.3476 —0.3820 — 0.5760 —0.6848 DISCUSSION 
0.50 —0.3861 —0.4237  —0.6375  —0.7581 
1.00 —0.4412 —0.4835 —0.7257 —0.8629 7 . . r > adi p 
200 —0 5134 —0 5622 —0 8424 —1 002 Examination of the results for the gradient of the wt 
3.00 —0.5641 —0.6174 —0.9245 —1.099 temperature function at the wall reveals a check with 
4.00 —0.604 —(0.6608 —0.9890 = 76 . . . : s#hs 
_ _ — hc the results of Chapman and Rubesin which is within be 
the accuracy of the present calculations. There is also th 
( ‘ -ement wit > re > f ‘ ation 
(3 = 0 (flat-plate flow) agreement with the results of Schuh for the stagn gee be 
% y : ~ = > fatenlate co: > ATE eV > Vv 
y = —0.5 to 4.00 flow, but his results for the flat-plate case are evidenth ai 
a 0.70, 2, 10, 20 oe Caeee ; wi 
The effect of the parameter 8 on the wall-temperature of 
48 = | (stagnation flow) gradient as revealed in Fig. 2 is small in the region 0 < m 
> = —lto4 8 < 1 but becomes significant at higher or lower values by 
¢ = 0.70, 1, 5, 10 of 8. For interpolation over the entire range of values, A 
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Fic. 2. Calculated temperature gradient at wall for wedge 
flow with variable wall temperature. Prandtl Number = 0.70 


the representation of the quantity (hx/k) V ux vy is 
more convenient, as indicated by Fig. 3. Figs. 2 and 3 
contain a curve indicating the values of 8 and y, which 
in combination produce a wall heat rate invariable in 
the direction of flow. 

The effect of the Prandtl Number can be simply ap- 
praised by cross-plotting results for given 8 values as 
exemplified by the results for 6 = 1 of Fig. 7, which are 
represented on Fig. 8 within a narrow range by the use 
of a suitable factor, o°-*°®. Treatment of all the avail- 
able results in this way yields the following factors: 


8 1.6 1.0 0 —0.199 
d 0.367 0.355 0.327 0.254 
Thus, 
hx/k . 
7 = B(m, y)o 
V Wyx/v 


where Bo” is given by Fig. 3 for ¢ = 0.70. 

The exponent of the Prandtl Number as obtained can 
be approximated on theoretical grounds. Any analysis 
that postulates a linear velocity distribution within the 
boundary layer, such as the analysis of Bond,’ yields 
a factor o’*. Analysis presuming a constant velocity 
within the boundary layer yields a value of o'”* typical 
of all “heat conduction” solutions. Tifford,* approxi- 
mating the velocity distribution for the separation flow 
by the parabola u/u,; = 0.10n?, obtained a factor o ’* 
Approximations of this type become more valid as the 
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Fic. 4. (a) Temperature distribution for stagnation flow with 
variable surface temperature. Prandtl number = 0.70. (b) 
Temperature distribution for flat-plate flow with variable wall 
temperature. Prandtl Number = 0.70. 
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Fic. 5. (a) Temperature distribution for wedge flow (8 


0.70. (b) Temperature distribution for separation flow ({ 
—0.1988) with variable surface temperature. Prandtl Number 
= 0.70. 


1.60) with variable surface temperature. Prandtl Number = 
5 = 


thickness of the thermal boundary layer is reduced, and 
this effect is produced by increasing the value of o and 
also, as examination of the present results reveals in 
Figs. 4 and 5, by increasing the value of y. Because 
of this effect, the exponents » given for the present 
results are approximations relevant to the range of ¢ and 
They do, however, 
1 


y considered in the calculations. 
follow the trends predicted, being near '/, for the sepa- 
ration flow and near '/; in the range 0 < 6 < 1 where 
the approximation of a linear velocity variation is jus- 
tifiable and the value does increase as 6 increases be- 
yond unity, where the velocity distribution within 
the boundary layer becomes increasingly more uni- 
form. 

Further investigation of the results reveals that, 
within the range of the calculations, the function B(m, 
y) can be approximated within +5 per cent by the 
equation 


( 2m \0. 104 
Bim, y) = 0.57 + 0.205) a 
: m+ | 


9 0.37 + [0.12m/(m + 1)] 
<¥7 ‘ 

( + ) 

m+ | 


This expression, however, fails at large negative 
values. 
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APPENDIX (A)--ACCURACY OF THE NUMERICAL 


SOLUTION 


Numerical Solution 


The use of punched-card machines for solving linear 
differential equations by means of a difference method 
was illustrated by Feinstein and Schwarzschild.’ Pro. 
ceeding in an identical manner, the constants @,, 
byt, Cn—1, and d,_; were calculated on the electronic 
computer. The value of 6; was obtained by setting 
6: = O, thus giving, for 8 = O and y = 0, 


6, = 0.9296657 
If one takes 0% = O, the corresponding 4, is 
6, = 0.9296657 


Both values are seen to agree within the sixth dee; 
mal, which confirms the convergence of the tempera- 
ture profile. 

The use of the difference method for solving differen 
tial equations has been discussed in references 10 and 
11. The size of the interval and the total number of 
steps NV must be chosen to ensure negligible errors in 
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Fic.6. (a) Temperature distribution for stagnation flow with 
variable surface temperature. Prandtl Number = 1.0. _ b 
Temperature distribution for stagnation flow with variable 
surface temperature. Prandtl Number = 5.0 
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HEAT TRANSFER TO 


In this analysis the interval e was 


the final results. 
taken as 0.10 mainly because the Hartree functions are 
tabulated for each 0.10 value of the independent vari- 
The interval was changed in one particular 
case from 0.10 to 0.05, and results of the same order of 
This smaller interval, how- 


able 7. 


magnitude were obtained. 
ever, meant successive interpolations of the functions 
‘fand f’. The resulting error was considered too large 
and the interval e was taken as 0.10 for the remainder 


of the calculations. 


Accuracy of the Solution 


An absolute check of the accuracy of the calculated 


results can be obtained by means of an iterative 


method. A more accurate solution can be evaluated 
by taking the machine solution as the first approxima- 
tion in the integral form of the energy equation. A 
more accurate temperature distribution thus follows. 
Such a calculation, however, was not accomplished. 
The possible error in the machine solution was esti- 
mated instead. 

As pointed out by Feinstein and Schwarzschild, the 
total error made in one step consists of two parts 
the error of the extrapolation formula and the rounding 
off error. 

The energy equation, Eq. (9), for 7 = (nm — l)je can 
be written 


6,-1° = o(2 — B)fn—1'On1 — fo8,—1' (20 


By developing 6, and 6, -» into a Taylor series, one 


gets 
, é ,? e sre 
6, = On-1 + CO,-1 + 5) On : + =e +- 1) 
° oO. 
/ e ‘Ad e AAA 
0, ’> 6, es ed, l + >! 6, 1 a py On 1 7 (22) 
- >. 


Eliminating 6, ;’ and 6,_,’" between Eqs. (20), (21), 


and (22), there results 


where the fifth and higher values of the step value e 
have been neglected and where a,_; and 6, ; are given 
by Eq. (14). 

The dominant error in Eq. (23) is 


Ee. he = ( fae? 2)6,, "ie (24) 


The next effective error is 
E, 2° = (foe*/6)0,1'"" (25) 
Every E,_, produces errors in the following 6's, in- 
creasing roughly linearly with the number of subsequent 
steps. Hence, the error in the final @ produced by 
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Calculated temperature gradient at wall for stagnation 
flow (8 = 1). Variable wall temperature 
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Fic. 8. Effect of Prandtl Number upon temperature gradient at 
wall for stagnation flow with variable wall temperature 


1 is (V — n + 1)E,-1, where N is the total 


number of steps. 


one £, 

If one sets F,-1 = (Ex-1) average these errors can be 
summed up by the square root of the sum of their 
squares. The summation over gives, for the error 
in the final 4; produced by all the /,,_:’s, 


Rm (Ea (iN SS (26) 


The value of R was estimated for 8 Oand y¥ 0 
It was thus found that R is approximately equal to 


0.011. This value was obtained by putting V = 61 
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The temperature gradient at the wall, however, can REFERENCES 


easily be determined by stopping at N = 50. The ' Seban, R. A., Calculation Method for Two-Dimensional Lap 
extrapolating error in (d0/dy)»» is thus of the order of — inar Boundary Layers with Arbitrary Free Stream Velocity Varin 
tion and Arbitrary Wall Temperature Variation, Institute of En 
gineering Research, Report 12, Series 2, University of California, 
Berkeley, 1950. 


0.0095 or about 2 per cent. 
The rounding-off error of 6, is produced by the errors 
Of Ay—-1, bp—-1, On-1, and 0,2 appearing in Eg. (15). Since 








: : : : ° 2 Hi « BD. 2. o Equation Occurring Y er — 
the calculations were carried to six decimals, each of 9 @tttee, VD. KR On an Aquation Occurring in Falkner ang 
eee t ‘ a r ff (1/2 Skan's Approximate Treatment of the Equation of the Boundary 
lese four terms has a rounding-off error of (1/2) X Layer, Proc. Camb. Phil. Soc., Vol. 33, pp. 223-239, 1937. 
10-®. The corresponding rounding-off error of 6, be- ae ; 
I 8 8 - 3’ Schuh, H., Laminarer Waermeuebergang in Grenzschichte) 
comes bei Hohen Geschwindigkeiten, ZWB_ Forschungsbericht 20}9 
oe : - “Heat Transfer in Laminar Boundary Layers at High Vy 
i = (a, 1 + b, 1 + 6, 1 + 6, »)(1 2)10°° (27) ties.” "se ' 8h Veloci- 
- 5 
Using dn) ~ 2, bps ~ 1, 0,1 ~ 1, 0,2 ~ 1, one gets ‘Chapman, D. R., and Rubesin, M. W., Temperature ang 
Velocity Profiles in the Compressible Laminar Boundary Layer 
— me 6 s . : ‘ ae ‘ ge ’ : 
i = 2.5 X 10°* with Arbitrary Distribution of Surface Temperature, Journal of 
. : : _ the Aeronautical Sciences, Vol. 16, No. 9, pp. 547-565, Septem 
Proceeding in a manner similar to the one used for 5. 1949 Exte 
the extrapolating error, the final error in 65, produced 7 : . Ble 
I 8 3 the final error in 66) I ; 5 Pohlhausen, E., Der Waermeaustausch zwischen festen Koe | 
by all the rounding-off is pen and Flussigkeiten mit kleiner Reibung und Kleiner Waer Robert 
. ‘ia ee Vs meleitung, Z.£.a.N.M., Bd. I, S. 115, 1921. Aerop 
S = T,(N3/3)'” (28) ' Ohi 
6 Eckert, E., The Computation of the Heat Transfer in the Lan Januar 
For 8 = 0 and y = 0, S is approximately equal to inar Boundary Layer of Bodies in Flow, V.D.1.-Forschungsheft 
7 X 10~*. Thus the total error in 6; is 0.0117, and the No. 416, Berlin, 1942. "9p 
error in (d6/dn)n—o is about 0.010. This last estimate 7 Bond, R., Heat Transfer to a Laminar Boundary Layer with P 
checks the discrepancy noted in (d6/dy)n» for B = Q Non-Uniform Free Stream Velocity and Non-Uniform Wall Ten Decen 
ss . : , , : > on oineer: spare . 0. Series 9 far 
and y = 0. Indeed, the computed value of —0.48S see go A sy —" Report 10, Series 2, ae 
‘ ore . ‘ ss Iniversity of California, Berkeley, 1950. Wi 
is 0.008 higher than the value —0.496 derived by Eck- ; i ie ate 
ert * Tifford, A. N., Some Prandtl Number Effects on the Transfer F thr 
. 0 
p ? of Heat, Engineering Experiment Station News, The Ohio State i 
It is seen that the extrapolation error is by far the ;... ’. 7 Algo Golan 
f ph: : 7 ; University, February, 1950. iad 
larger one of the two. This extrapolation error in- _ ; ; : as 
cans Wi 0 1 -_? ® Feinstein, L., and Schwarzschild, M., Automatic Integration of The 
pene ‘ oe — eT oe oo aoe Linear Second Order Differential Equations by Means of Punched from | 
For y = 0, the error in (d6/dy)»-o must be less than 2 Curd Machines, The Review of Scientific Instruments, Vol. 12 sultan 
per cent, while for y < 0 this error may reach 3 to 4 No.8, p. 405, August, 1941. trifug: 
per cent as 6 converges rather slowly for negative values Boole, G., Treatise on the Calculus of Finite Differences, New sae 
of y. Higher values of Prandtl Number will tend to York, Stechert (1946 reprint of 1872 edition). 
decrease the error in the temperature gradient at the '! Jordon, C., Calculus of Finite Differences, New York, Chelsea 
wall. (1947). 
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Extension of ‘‘The Dynamic Effects in Rotor 
Blade Bending”’ 


Robert A. Rogers 
Aeroproducts Division, General Motors Corporation, Vandalia, 

Ohio 
January 9, 1952 

ONSIDERING THE PAPER by A. H. Flax and L. Goland ap- 
C pearing in the JOURNAL OF THE AERONAUTICAL SCIENCES, 
December, 1951, it is interesting to extend this analysis to the case 
of a rotating propeller where the blades may be near resonance 
with the periodic aerodynamic forcing function caused by the 
air stream entering the propeller disc at an angle to the direction 
of thrust. Results are identical to those obtained by Flax and 
Goland, although the boundary conditions of the problem are dif- 
ferent. 

The fundamental equation of blade bending may be derived 
from consideration of Fig. 1, where f represents the blade re- 
sultant transverse load per unit length and where the axial cen- 
trifugal load is represented by p. The equation is, after dif- 


ferentiating twice with respect tor, 

o? __ d*y oy O*y °R 
z (21 *) 2h = fo >= 
or? dr? or or? 


where the transverse loading is composed of three parts: h, the 
aerodynamic loading considering a rigid blade; t, the loading 
and ¢#;, the increment of 


p dr (1) 


r 


due to the weight and inertia forces; 
aerodynamic loading caused by deflection of the blade. For 
propeller blades, the weight term is negligible in comparison with 
the aerodynamic and inertia loadings. The increment of aero 
dynamic loading may easily be derived by consideration of the 
increment of transverse velocity of the air relative to the blade, 
AV,, produced by the bending motion of the blade. Consider 
ing Fig. 2, 


U? = (rQ)? + V? (2) 
Aa = AV,/U (3) 


where Aq@ is the incremental angle of attack caused by bending 


deflection of the blade. The increment of aerodynamic loading 


is 
ts = —aocp/2U*Aa (4) 
= —apcp/2UAV, (5) 


Since AV, = Oy/Ot, then 


t; = —<doCp 2U0d, ol (6) 
Letting 
—a = —aocp/2|(r2)? + V2]'/? (7) 
Then 
tj =— a(r) oy ol (8) 


Letting 


*R 
F(r) = | p dr 


we arrive at the basic equation of blade bending which is the 


same as Eq. (5) in Flax and Goland’s paper. 


o? _,d*y  . OF o*y oy oy 
EI — F(r) + om + p +a = it, (9) 
or? dr? or? or? or ot 


They have solved this equation utilizing the boundary conditions 
of a rotor blade hinged at the root. Considering a propeller 
blade, the equivalent boundary conditions are 








ds = O¢;,/Or = V (r = 0) 
re) Or | EI( 0d, Or?) | =() (r = R) 
0°¢,/dr? = 0 (ow BY (10) 
F(r) = 0 (ry = R) 
T, (1) 
Mp 
a py) 
Mt 
Y 
y 














Fic. 1. Loads acting on deflected propeller blade 
AV, 
4 
Fic 2. Velocity diagram of oscillating airfoil section 


349 








350 JOURNAL OF THE 


*R 9 9 9 

a 0°; _, « Mb Op 

(a1 )e F(r) ds + p 
or? or? or? 


AERONAUTICAL 


ds + bs? (—R2Q2m 4+ 
r 


SCIENCES MAY, 1952 


und these are used in integrating Eq. (18) of Flax and Goland’s paper which is 


*R 
tRQa) | dr = tpbdetR2t dy 


0 


Integrating the first term in Eq. (11) by parts twice, the first two terms again vanish so that 


“R 92 __ O*ds 
- | El — } ¢, dr 
or? or? 


0 


"R 026, \? 
= El ( ) dr 9 
0 or? "| 


Similarly, the second and third terms of Eq. (11) are combined and integrated by parts so that, considering the boundary conditions of 


Eq. (10), 


- 5 OO Od, “R O¢, \? 
—- F(r) os + Pp d, |dr = F(r) dr 
0 or? or 0 or 
Eqs. (12) and (13) are substituted into Eq. (11), and the result is identical to that of Eq. (24) determined by the authors. Thus, fo; 
propeller blades, the complete expression for deflection under load ¢,, is 
*R 
t,o, dr 


0 
« 


es 


pa “ ie a 
1g ¥(! ) ; 


and the static deflection or bending moment may be multiplied by an 


n 


A.F. = 1/¥V [1 


to obtain the dynamic deflection or bending moment, respectively. 


Attachment of Rotor Blades by Laminated 
Elastic Straps 


Kurt H. Hohenemser 

Chief of Aerodynamics, Helicopter Engineering Division, McDonnell 
Aircraft Corporation, St. Louis 

January 23, 1952 


ibe ATTACHMENT OF A ROTOR BLADE by two laminated elastic 
straps offset from the blade feathering axis in the plane of ro 
tation first used in Doblhoff’s experimental pressure jet 
helicopter WN 342 and is also incorporated in the French pres 


a substitute for a 


Was 
sure jet helicopter ‘‘Ariel."’ The straps are 
flapping hinge and for a thrust bearing otherwise necessary to 
transfer the centrifugal force of the blade to the hub. In com 
bination with a floating hub the straps are a convenient, simple, 
and trouble-free means of blade attachment which also allows 
the vertical blade hinge and its damper to be eliminated. 

When designing the strap attachment, care has to be taken to 
avoid mechanical instability of the chordwise blade oscillations 
in flight and on the ground, and the spacing, length, and cross 
section of the straps are powerful parameters for influencing the 
natural frequency v, of the lowest chordwise blade bending mode 

Reference 1 gives a description and analysis of the strap at- 
tachment for the French helicopter ‘‘Ariel II.’’ The following 
criterion for avoiding chordwise self-excited mechanical oscilla 


tions of rotor blades is used: 


2>w0/"> 1 (1) 
This criterion was taken from refer 


Reference 2 is based 


where w is the rotor speed 
ence 2, but it is, unfortunately, in error. 
on reference 3, which gives a correct analysis of the problems of 
self-excited chordwise rotor-blade oscillations including blade 
damping using the frequency response method. The results of 
reference 3 are in agreement with those of reference 4. In refer 
ence 2, however, an erroneous transition was made from the case 
with finite blade damping to the limiting case of infinitely small 


blade damping. Actually, the center of the instability range in 


| " EI (Ce) a + = F(r) (= ) dr 
a or? 4 or? | 


@,\? 
992. 2 

R226? 

oo” 


amplification factor 


the case of undamped chordwise blade oscillations is given by th 
relation 


V; - vp, = WwW 


where v, is the frequency of a horizontal pylon mode. In flight 
this frequency is identical with the so-called pendulum frequency 
of the rotor craft and is quite low. Therefore the center of me 
chanical instability is avoided if 


) 
Wmas vi 2 


for the highest rotor speed obtainable in flight. This is certainly 
a safe design criterion for the strap attachment of rotor blades 
In case of the French helicopter ‘‘Ariel II,”’ this criterion is ol 
viously not fulfilled. In this case the landing gear has to be 


sufficiently stiff, so that on the ground 


Wmar. S< Vb TT VI 


where vy, is the lowest frequency of a landing gear mode witha 


horizontal component at the location of the rotor. In the air 


the rotor speed must never drop below 
Wmin , Vy T 


where v, is the frequency of the pendulum mode, otherwise mé 
chanical instability of the chordwise blade mode will occur in 
flight. When designing the strap attachment the erroneous cri 
terion (1) must be replaced by the criterion (2) or by the criteria 
(3) and (4) for the center of the instability range, and an adequate 


- 


margin Of wnar, and wmi,n, Must be kept with respect to this centet 


in order to avoid mechanical instability 
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The Shock Strength in a Two-Dimensional velocity of the undisturbed fluid)t, with apex at r = 0,¢ = 0 

Nonsteady Flow In order to obtain an approximate solution uniformly valid on 
the Mach cone, Lighthill’s technique" ? is applied 

Ly Ting Expanding both the independent variable r Cot, and the 

Daniel Guggenheim School of Aeronautics, New York University velocity potential g into power series of a small parameter e, 

December 28, 1951 r— Ct = R+ eo %(r, 0) + r@(r, 0) +... (1) 
SSUME THAT A SMALL DISTURBANCE Starts at origin (7 = 0); Ar, 6, t) = ep Mr, 0, R) + Her, 0, KR) +. - (2) 

A when ¢ = 0, the linearized theory will give a good approxi and substituting them into the basic equations (in polar coordi 

mation to the flow in and outside the Mach cone, r = C)(sonic nates) for unsteady two-dimensional potential motion, we have 


0 a. aa a Coattt 4. coh 2: oes = ¢2{ [9y7.(1 1 v/C 2o.()/C] 
T Re, i €¢ e“¢ ev¢ Tese) @ Cilla (y + \Pr ( a¢ 0} X 
yt ror 28? drOR | rOR 
(¢pr™ + evrr® +...) + €[2r, — (y + 1)eRp™/Co — 2¢-/Co + a quadratic form in ¢,", ge, 7,, and 
ro opr + o(e2) + additional terms with at most first derivative with respect to R{ (3 


The first order equation is a simple wave equation for ¢'”, (y + 1)? 
which is exactly the same as the linearized solution of ¢(t, 7, @) a ae 2 (F + B*) (R 0) | 
with ¢ replaced by (r — R)/Cy. Near the Mach cone (R — 0), (7 
yg” and its first derivative behave in general as r(2 mi Ss (42 F) (R<0) 
. Co? 
yg) = D(r,@) + E(r, 0)R + O(R (RS 0)) ’ 
gr? = E(r, 0) + A(r, 0)(—R)*? + OCR) (R<O \ (4) The shock conditions are Ag, = 0 and Ag, = [2/(y + 1 . 
+ 1 
ge’) = E(r, 0) + Bir, 0)R + OCR) (R> 0) (C?/qn) — q,| where 
Since ¢gpr‘’ tends to infinity as R — 0, it is necessary that the a te (20, + o? + evr?) (4 1)/2 
coefficient of grr", grr’, in Eq. (3) vanishes at R = 0 in order 
to have finite gp”, gp. This condition together with the g = (ee PO S)/(TVM SF + Ser 
i onditi op?) = rp) = 7 D= - Te 
initial condition, Agr sates Oat ¢ = 0, and the result ¢. = (a + oS + eS VSP + Sor 
Ag?(R = 0) = O, to be verified later, leads to 
~, ] Substituting Eqs. (1) and (2) into the first shock condition, we 
oy al Wr. @) - ° ‘ es . 
pt) ae = . E(<. 60) dz 4 : have Ag’? = 0 and Ag = 0, which was used in deriving Eqs 
2Co 0 Ci (5). With Eqs. (1), (2) (3), (5), (6), and (7), the second shock 
: | *) (5) condition leads to 
Agr? = —— 2/2 1A%s, 0) + B%z, 6)| dz oC 
ee AV A? — F+ BV B+ F 4 Aen = 
y+ + 
Ar) = [(y + 1)2/4C,2] (A B?) 
me »( F B (8 
where 
. 7 The shock strength is given by 
l or 
A* = NA dZ  /t*A? dz Ap 1 2C., 
2 ie ) ” 2 22 , 2 r 
~ Jf? J? — - —~— i AW A? I BV B?+ F4 Agr (9 
; Pp Co? Cy Y t | 
and B? is related to B in a similar manner, and where A( a , ; es wer i 
. , “ee . . where F is the solution of Eq. (8) satisfying the initial condition 
indicates the change of a quantity from R > 0 to R < 0. rE 
‘ze ()atr = U. 


The last equation in Eqs. (5) indicates that Ar‘ is positive. 
Following Lighthill’s argument,? the shock front can be repre The result is applied to study the reflected shock in a Mach 


sented by the equation reflection raised when a plane shock wave, traveling parallel to a 





St, 7,0) = Col r+es%(r, 0) + & s® + Ole = (0 (6) 
with s) = 7) and r2(R < 0) > s® > r2(R > 0). For con- ae Shock 6) (¢) 
venience, an unknown function F(7r, @) is introduced as defined by ee Shock Oe?) 


— ~ --—- Shock ole*) 
PRESSURE PR, —== Expansion wave 0(*) 
DENSITY 2B -- Slip stream ole4 
SONIC VEL. G > MG 

















4, 
F / ~ 
— \O- 











plane wall, meets a shallow bump, y = « f(x), in it. 
maximum slope of the bump, as expansion parameter, the first 


order solution’ for ‘‘supersonic flow’ (the wall is moving with 


supersonic velocity M Cy relative to the undisturbed flow behind 
the shock) gives 
; p 
iim — 


r—>Cot poCo? 


M — cos 0 


ee mos gf (O* HCO)" "ReW/ Cot — 1) 


(10) 


W(0)r-”{ReWV Cot — 1) 
for 8 > 0, = m — arcsec Mo, and for 6 < Hj, the right side of Eq. 
(10) will have an additional term M)?8~'f’[rv(@)] where B = 
/ M2 — 1, v(@) = cos@ — Bsin @ + Mo, and H(@) remains posi- 
tive for the entire range of 9@and Mp. 

Since the flow is irrotational outside the triangle O7J, a po- 
tential function y can be introduced to describe the flow near the 
Mach circle including the triple point 7 (@ = arccos M), where 
both p“ and its first derivatives are zero. Comparing Eqs. 
(10) and (4), we have A = W(0)r-'2, B = 0, E = Cy)M,28 “If! 
[rv(0)], D = —M.C B-'f[rv(@)] for 6 < 6, and A = W(6)r-'/2, 
B= E =D = 0for@>4,. 

Solving Eqs. (8) and (9), we have the strength of the reflected 
shock 

Ap Ap | 3(y + 1l)y 


we)? if W> ol 


Po p C,? at 16 (11) 
= 0 ifW< of 
The shock front is represented by the equation 
+ 1)Mo 
r-—-Cit= t =: = -- 1 fe sro + O(c?) (for 0 < 6) 
(12) 
r — Cot = &(y + 1)?W%(0)r/4 + O(e*) (for @ > @;) 
(13) 


The possible shock patterns at an instant ¢ > 0 for f’(0*) greater 
and smaller than zero are shown in Figs. 1 and 2, respectively. 

At point L, @ = 6, where the Mach circle is intersected by a 
first-order shock or expansion wave, the theory fails.? 

When the flow is region III vanishes and Eqs. 
(11) and (13) are true for the entire range of 6. The result agrees 
with Tan’s.* 

For all cases, at the triple point, W = 0, so the strength of re- 
flected shock is O(e?). As a consequence, a three-shock con- 
figuration with a slipstream does not exist there up to the second 
order—i.e., the strength of the slipstream is O(e?). This state- 
ment can be proved more rigorously by using the second-order 
adiabatic equation and the shock equations across the refracted 
shock. TJ. 


“ 


subsonic,” 
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The Reverse-Flow Theorem for Nonstationary 
Flows 


A. H. Flax 
Cornell Aeronautical Laboratory, Inc., Buffalo, N.Y. 
January 24, 1952 


F@ STEADY FLOWS, reverse-flow theorems valid within the 
limits of linearized theory have been presented in various 


With e, the 


forms and for various cases by von Karman,! Hayes,? Munk 
Harmon,® Brown,® and the writer. The most gener: fi 
‘ seneral forms of 

these theorems have been given by Ursell and Ward? and ind 
° e ° : de- 

pendently by the writer® using different methods of attack, Py 
e 5 . . *ur- 

ther work on the subject has disclosed that there is a rever 
ie TSe- 
flow theorem for nonsteady flows as well. 

stated (for harmonic oscillation ) as 


SS wx, y, t)w(x, y, t) ds = S SF vl x, y, t)w(x, y, t) ds 


where is the instantaneous value of lifting pressure and w is th 
° . e a ec 
instantaneous value of prescribed vertical velocity at any point 


The theorem May be 


on the wing plan form. and w are the same quantities for a 
flow around the same plan form with the direction of the free. 
p and p and w and w are so de- 
fined as to have the same sense—i.e., p and p are both POsitive 
for lift upward and w and w are both negative when the local 
angles of attack with respect to respective free-stream directions 


stream velocity, uo, reversed. 


are positive. 

This result can be obtained formally from simple basic con 
siderations. Consider a thin wing of any plan form harmonically 
oscillating in an incompressible flow. The boundary of the wing 


is described by a function 


s = 2,(x, y, t) = 20x, ye 


The velocity potential for such a wing satisfies the equation 
V*@ = 0 except on the wing itself and on the vortex sheet trailing 


behind the wing. The boundary conditions are 


U,(020/0x) 


twl 


+ 1w3Z9 l 


Wo = Ogd)/Os = 
where @ = goe'™ and w = woe The pressure at any point (to 
within an arbitrary additive function of time), according to the 
linearized Bernoulli equation, is given by 
+ twp 2) 


7 = p/p = Ut + (0¢/Ot) = Uou 


Here the velocity perturbation components are given by 


Od Od Od 
“z= . y= » w= ) 
Ox Ox Oz 
The flow is irrotational, and, in particular, (Ou/Oz) — (Ow/dx) = 


0. wo is the unperturbed free-stream velocity in the x direction 
Over that part of the plane zs = 0 occupied by the wing, there isa 
pressure discontinuity Az between 7 on the upper and lower sides 
of the plane, but everywhere else in the fluid Aw must be zero 
At the trailing edge, the Kutta condition requires that Ar = 0, 
in order that z go smoothly into its value in the wake. Since u 
and ¢ satisfy Laplace’s equation, 7 satisfies Laplace’s equation 
as well. 

Now, if there be a flow over the same plan form with the direc 
tion of the free-stream velocity reversed, there are associated 
with this flow functions ¢ and z, which have properties and in- 
posed conditions similar to those of ¢and x. For the reverse flow 
case, however, the wake trails from what was in the first case the 
leading edge, the Kutta condition being applied there. It should 
be noted in each case that across the plane z = 0, Ad and Au differ 
from zero only on the wing and in the wake. In each case, there 
is a singularity in « and x of the type 1/+/d at the leading edge, d 
being the distance from the edge.®~* 
The volume integral 


SS fore — $V %r) dv 


may be written, where V, the volume of integration, is taken overall 
space outside the wing. It is to be noted that V*¢ = Oand V%p =0 
everywhere outside the wing except in the wake; as x approaches 
+o in the wake of the direct flow, ¢ and x go to zero; while, as 
x approaches — ~ in the wake of the reverse flow, ¢ and x go to 
zero. Further, at the leading edge, where the derivatives of 
¢ are singular, w is zero as required by the Kutta condition in 
the reverse flow, and similarly for + and ¢. Thus, Green's 


theorem may be applied to this volume integral, using as the 
to x = 


= —a 


boundary of the fluid the plane z = 0, from x = 
+o, Then, 
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(1rV*o — oV2x) dv = f Cg »” ds 
Vv Ay, A- On 

Or 

Sf. ¢ 

Ay, Ae 


On 
where the area of integration, A, is the entire plane z = 0 from x 
_otox = +. The subscript + denotes the positive side of 
the plane, while the subscript — denotes the negative side. The 
» direction is positive for the outward normal of the bounding 


ds =0 (4) 


plane, s = 0. Further, 
oo Op Od dg 
(%*) xs as" on dl he 
$+ = —¢-; i -— —7 


because of the antisymmetry of the lift problem. Thus, 


1 in ds — rr) = ds = 0 
A 02 A OZ 


where only values on the + side of the plane z = 0 need now be 


(5) 


considered. 


Or Ou . 3 
= —U + lw (6) 
Oz Oz Oz 


It should be noted here that the sign of the first term is negative 
because, with x defined as positive in the direction of the unre- 


versed flow, u = 0¢/Ox has a sign convention opposite to the 
usual one. If the senses of x and w are kept the same for the two 
flow directions, the signs must be chosen as in this note. By the 
irrotationality condition ‘ 
Ou/Os = OW/OX 
Also, 
0¢/0s = Ww 
so that 
(Or/Oz) = —uo(Ow/Ox) + tww (7) 
Thus, 
— — 
Op ww , . 
7 ds — ri) Uo + ipww | ds =0 (8) 
e Oz SF ae Ox 
The integral 
oe 
Uo ff - (Ow /Ox) ds 
may be integrated by parts on x; since atx = —~, @ = O, 
while at x = + ©, w = 0, there results 
2 fe 
u | fo — ds = —uU ew ds = 
J JSA OX A Ox 
— Ug uwds (9) 
A 
Substituting this result in Eq. (8), 
(10) 


n ~~ 
Od : 
1 ds — (uu + twh) wds = 0 
4 O2 


But according to Eq. (2), the term in brackets in the second in 


tegral is simply z, so that 


S S rw ds = S S ww ds 


where the integrations are over the wing plan form only, since 
This is the desired re- 


(11) 


and zw are both zero except on the wing. 
sult that reduces to the theorems of Ursell and Ward’ and the 


writer® for steady flows, 


S Si w ds 


(12) 


= S S uw ds 
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The present result is also true for compressible nonsteady flows. 
The extension of the proof will be given in a later paper. For 
the special case of two-dimensional incompressible nonsteady 
flows, Schwartz® has given formulas for the lift and moment ex- 
plicitly in terms of w(x, ¢). His results are exactly those ob- 
tained by the direct application of Eq. (11) to the case he con- 
sidered. 

An adjoint variational principle of the type discussed by the 
writer for steady flows’ may also be derived for unsteady flows 
by methods entirely analogous to those for the steady-flow 


case. 
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Comments on ‘‘Aerodynamic Coefficients of an 
Oscillating Airfoil in Two-Dimensional 
Subsonic Flow’’ 


Henry E. Fettis 
Mathematician, Flight Research Laboratory, Wright Air Develop- 
ment Center, Wright-Patterson A.F. Base, Dayton, Ohio 


January 7, 1952 


ys THE SUBJECT PAPER,' the authors give values for lift and 
moment coefficients for a two-dimensional oscillating air 
foil moving at subsonic velocity. A cursory examination of the 
values presented indicates that they are subject to considerable 
question, particularly at higher values of the reduced frequency, 
w. The authors themselves have noted that the values agree 
poorly with those previously computed by the Dietze method, and 
it is difficult to believe that the latter results are greatly in error, 
since they have been derived independently? by a different 
method, the results agreeing except for differences that are 
clearly of second order in magnitude. Since the publication of 
reference 2, additional calculations have been made‘ which fur 
ther confirm the correctness of the Dietze values. 

In Figs. la, and 1b, a plot is shown of the lift coefficient L, as 
defined in reference 3, which is related to the quantities kg’ and 
k,” in reference 1 according to the equation 


Ly, = (1/w*) (k,’ + ik”) (1) 


Fig. la indicates fairly good agreement with Dietze down to a 
value of (1/w) equal to approximately 1.5. For 1/w < 1.5, 
a marked deviation in the real part takes place which is difficult 
to account for in view of the fact that the Dietze values for this 


quantity approach the known value of | as 1/w — 0. 








354 JOURNAL OF THE 


AERONAUTICAL 














Fic. la. Real part of Li vs. 1fw; B= 0.7. -Dietze; A 
Timman 
/ 2 3 ¥ s é > oe i, 
4), ° = =e z ' . ay = 7 . 
me 
° 
a 
s a" 
ve 
Fic. lb. Imaginary part of Li; vs. 1/w; B = 0.7 Dietze ; 
Timman. 


A further difficulty lies in the fact that the various coefficients 
are not mutually in agreement. This may be demonstrated with 
the aid of the following recently demonstrated relation: 

k,= iw(m, + k,) (2) 
A spot check of the values in Table 1 of reference 1 indicates that 
for the first few values of w at each Mach Number, the above rela- 
tion is satisfied with a reasonable degree of accuracy, the dif- 
As w and M 
increase, values on the right-hand side and left-hand side of Eq. 


ferences being possibly due to rounding-off errors. 
(2) become widely separated. On the other hand, the values of 
Dietze check reasonably well by Eq. (2), with discrepancies in 
The 
puted in reference 5 satisfy Eq. (2) to the fourth or fifth digit, in 


only the third or fourth significant figure. values com- 


general. 
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Comments on ‘‘Uniformly Loaded Semi- 
Infinite Wedge-Shaped Plates’’ 


M.L. Williams, Jr. 
Assistant Professor, Guggenheim Aeronautical Laboratory, 


California Institute of Technology, Pasadena, Calif 
February 11, 19592 


- A. J. A. MORGAN’S RECENT NOTE (Readers’ Forum, JOURNAL 

OF THE AERONAUTICAL SCIENCES, Vol. 18, No. 12, p. 845 
December, 1951), it is stated that ‘‘the deflections of a uniformly 
loaded semi-infinite wedge-shaped plate of constant thickness 
are explicitly determined for the case where one edge is free and 
the other is clamped.”’ Then, in one of the concluding para- 
graphs, it is stated that, for the boundary conditions considered, 
“no stress singularity occurs at the corner for any wedge angle 
0 < 28 < =.” 

While both statements are correct, it is believed that the man- 
ner of presentation leaves an impression that the stresses at the 
corner are always nonsingular. The problem is discussed in ref- 
erences 1 and 2, where it is shown that any attempt to prescribe 
arbitrary boundary conditions other than those compatible with 
Eq. (17) of Morgan’s Readers’ Forum item will indeed give rise to 
a mathematically infinite stress at the corner for certain opening 
angles. In particular, while Eq. (17) is the same form of solution 
given in reference 1, it may be shown that, if a free edge is desired 
along a circumferential boundary, for example, the stresses at the 
vertex are unbounded for opening angles greater than approxi- 
mately 90° 

It should be emphasized then that the corner stresses for this 
uniformly loaded cantilever plate (and actually even under more 
general loadings) will only be nonsingular in certain cases. The 
special solution, Eq. (17), is therefore thought to have somewhat 
limited practical applicability. 
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Surface Stress Singularities from Various 


Illinois Institute of Technology, June 


The Superposition of High-Intensity Periodic 
Waves on a Steady Flow 


J.C. Truman and N. J. Lipstein 

General Engineering Laboratory, 
Schenectady, N. Y. 

January 24, 19592 


Genera! Electric Company 


| eerecnet IN OSCILLATORY FLOW PHENOMENA has prompted all 
investigation of experimental methods for superposing high- 
intensity audiofrequency pressure oscillations on a steady flow 
Pressure oscillations that are treated in the field of acoustics are 
usually assumed to be periodic and to have small amplitude, pet 
Shock waves, 


on the other hand, may involve large pressure ratios but are in the 


mitting the use of the linearized wave equation 
usual case nonperiodic. Finite amplitude effects, which are found 
to change the shape of periodic waves of high intensity, are dis 
cussed in the present note 

It is readily seen! that a wave of finite amplitude cannot be 
propagated in a gaseous medium without change in wave shape 
From a consideration of the exact equation of motion, Fay? has 
shown that a transfer of energy from components of lower fre- 
quency to those of higher frequency takes place due to the non 
linearity of the system and the interaction of harmonics. The 
effect of this transfer is counteracted by that of viscosity, which 
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attenuates higher frequency components more than those of 


lower frequency. Hence, for any arbitrary intensity and fre 
quency, there must exist an equilibrium wave form, having har 
monic content in such relationship that the decrease in amplitude 
of any component due to viscosity is exactly compensated by the 

Any wave approaches this equilib 
Conditions for stability are found to 


increase due to nonlinearity. 
rium form as it progresses 
depend on intensity, however, so that the equilibrium form is not 
permanent but changes as intensity decreases. Experimental 
evidence of finite amplitude effects, substantiating Fay’s theo- 
retical analysis, has been reported by Allen and Rudnick.é 
Data obtained by the present authors for high-intensity waves 
superposed on a steady flow also indicate the presence of such 
effects.‘ 

Iwo different ‘‘flow-chopping”’ devices were used to modulate 
. A gate 


the airflow in a smooth pipe 3°/, in. in inside diameter 


modulator, consisting of a rotating wheel with slots equally 
spaced around its periphery, was mounted so that the flow exit 
from the pipe was alternately opened and closed. A sleeve modu- 
lator was formed by two coaxial cylinders having aligning port 
configurations in their walls. As the outer cylinder was rotated 
about the inner one, which formed a section of the pipe contain 
ing the air stream, the ports were alternately opened and closed 


A Massa 


microphone was used to obtain pressure responses at six stations 


Both modulators were belt-driven by electric motors. 


along the pipe 

Fig. 1 shows oscillograph traces obtained at stations 1 through 
6, located at distances upstream from the gate modulator equal 
to 65, 57, 50, 31, 24, and 16 in., respectively. The change in wave 
a triangular wave near the modulator to a sharp 
Here the velocity of the 


shape from 


sawtooth at station 1 is clearly seen. 


air stream was 315 ft. per sec., and the signal frequency was 275 


’ 


cycles per sec. The average sound intensity was 167 db., cor 


responding to a pressure amplitude slightly less than 1 Ib. per 


sq. in. Fig. 2 illustrates the wave shapes at stations 1, 2, and 3 


located 24, 16, and 9 in., respectively, upstream from the sleeve 


modulator. Air velocity was 380 ft. per sec., signal frequency 
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was 340 cycles per sec., and average intensity was 166 db. Here 
also, the change to the sawtooth wave form is noted. 

This work was part of an investigation performed at the re 
quest of the Propulsion Group, Flight Research Laboratory, 


Research Division of the Wright Air Development Center. 
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Recovery Factors for the Laminar ‘‘Separation’’ 
and Stagnation Flows* 


S. Levy and R. A. Seban 
Research Engineer and Associate Professor, respectively, Institute 
of Engineering Research, University of California at Berkeley 


January 30, 1952 


INTRODUCTION 


|, penton WALL TEMPERATURES for the two-dimensional in- 


compressible ‘‘wedge flows’’ corresponding to the Hartree! 
solution of the momentum equation have been computed by 
Eckert and Drewitz? and Schuh’ and can be inferred from cer- 
tain calculations of Tifford.4 Some of these results are shown 


in Table 1 in terms of the recovery factor 


r= 1 — [((7, — Taw)/(U:?/2GJIe, 
TABLE | 
Recovery Factors for ‘Wedge Flows’ at a Prandtl Number of 
0.70 . 
B 1.0 0 —().14 —0.20 
Eckert-Drewitz 0.791 0.833 0.838 
Schuh 0.840 0.840 
Tifford 0.830 0.835 
Present results 0.815 ace 0.762 


The theoretical results now reported were obtained to estab 
lish the value of the recovery factor for the velocity profile in the 
boundary layer corresponding to the value 8 = 0.20, which 
yields zero shear at the wall and, in consequence, corresponds to 
conditions at separation of the boundary layer. In connection 
with this work, the recovery factor for the case 8 = 1.0 (stag 


nation flow ) was also determined 


NOMENCLATURE 


b(n) = function defined by Eq. (7) 

B,, = coefficient obtained from Eq. (9) 
c(n) = function defined by Eq. (7) 

ty = proportionality constant (UU) = Cx 
j = coefficient obtained from Eq. (9 


¢ = 


= step interval in the finite difference solution 


F, = extrapolation error 

f(n = solution to Eq. (3) (Hartree profiles 
2(n) = function defined by Eq. (5) 

G = gravitational constant (32.2 ft./sec.? 
J = conversion factor (778 ft.lbs. per B.t.u 


* This work was carried out at the University of California under the 


sponsorship of the the Air Materiel Command (Air Force Contract No. 33 


0038-12941 
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Fic. 1. Temperature field in the insulated laminar boundary 
layer for ‘‘separation’’ and stagnation flows with Prandtl Num- 
ber of 0.70. 


k = thermal conductivity (B.t.u./hr. ft.2 °F. /ft.) 

m = exponent in the statement (U; = Cx”) 

M, = local free-stream Mach Number 

p = pressure in the boundary layer (Ibs. per sq.ft.) 

Pn» In = Coefficients defined by Eq. (13) 

r = recovery factor 

= temperature in the boundary layer (°F.) 

a = temperature of the insulated surface (°F.) 

7, = stagnation temperature (°F.) 

Ti = temperature at the edge of the boundary layer (°F.) 

u = velocity component of flow in x direction (ft. per 
sec.) 

U; = velocity at the edge of the boundary layer (ft. per 
sec. ) 

v = velocity perpendicular to body surface (ft. per sec.) 

x = distance along the body surface (ft. ) 

y = distance normal to the body surface (ft.) 

8 = velocity distribution parameter 

n = independent variable occurring in the Hartree pro- 
files 

0 dimensionless temperature [(7 T1)/(U?/2GJc,) | 

m = fbsolute viscosity (Ibs. sec. /ft.?) 

v = kinematic viscosity (ft.?/sec. ) 

o = Prandtl Number, uc,/k, dimensionless 

y the stream function 


ANALYSIS 


The momentum and energy equations upon which the solutions 
considered here are based are written in the following form: 


Ou Ou Op O7u 
plu tv + =u (1) 
Ox oy Ox oy? 


oT oT k OT nw fouy? 
plu tv = + (2) 
Ox oy Cp OY? Cy \Oy 


Eckert and Drewitz have presented a consideration that indi- 
cates that compressible flows may be defined by these forms 
of the equations for Mach Numbers in the range 0 < J/, 
0.9. 


Transformation of Eqs. (1) and (2) through Hartree variables 


n=V 11 + m)/2\(Ui/vx)y 
Y= V (2/11 + m)] vx U1 f(n) 
gives, for U; = Cx™ and B = 2m/(m + 1), 
fo" bap + al f'?) =(Q (3) 
0’’ + fod’ 2Bof'd + 2of’’? = 0 (4) 


with u/U; = f’(n) and An) = (T T1)/(U1?/2GJc,). 


* The primes indicate the order of differentiation with re spect to n 
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The boundary conditions are 
f’ =0 6’ = 0 n=y=0 
f' =1 6=0 n=yo> 


Introduction of the function of g(7) defined from 
@ = ge ‘1/2) So" of an , 


yields, for the energy equation, 


g’’ = b(n)g + c(n) P 
where 
b(n) = 2Baf’ + (1/2)ef’ + (1/4)o?f? , 
c(n) = —2Qaf’’%e" 2) fy" af dn 


The boundary conditions become 


g(0) = 0; 


gin > ~) = 0 


Reduction of Eq. (6) to finite difference form, with the increments 
in g being distant e apart and subscript » denoting the number of 
the increment beginning with zero at the plate, produces th: 
equation 


nor = Bagn — gna tC, g 
where 
B, = 2 + e*b(n) 9 
Cr = e*c(n) 


The calculations leading to g, for large n were accomplished 
on a calculator in which five decimals were carried. The incre. 
ment e = An was taken as 0.1. 

The boundary condition g’(0) = 0 is approximated by a 
parabolic curve in the same manner as by Eckert and Drewitz 


Thus, 


2 Ci 
$— B, 4— B, 


7 
ll 


Thus g and all other g, can be expressed in terms of go or 


S 


Ln = Pr¥o In é 


where 
Pn = Brapna Pn—2 13 
qn = B, 19n-1 — Gn-2 + Cans 
with p, = 2/(4 B,),, fbb = 1,4 = G/(4 — B,), an 
go = 0 


In virtue of the boundary condition as 7 ~ ©, g, becomes zer 
for extremely large n. In the calculation if the value of go was 
unaffected within the fifth decimal by evaluation from Eq. (12 
and with g, = 0 for two positions m and n 1, the number of 
increments was deemed sufficient. 


The recovery factor is obtained as 
r= 00) = 2(0) = go l4 
and can be evaluated from the known value of go. 
RESULTS 
Calculations have been carried out for the following cases: 


B= 1.0 (stagnation flow) 
8 = —0.20 (separation flow) 


for a Prandtl Number a of 0.7. 

The distribution of the temperature within the boundary layet 
is shown in Fig. 1, and the computed values of the recovery factor 
are shown in Table 1: r = 0.815 for B = 1.0; r = 0.762 for 8 = 
—().20. 
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DISCUSSION Mi, Reynolds Number R, specific-heat ratio y, and Hilton’s tem- 
The present result for the stagnation flow as revealed from perature recovery coefficient Co. 
Table 1 lies between the value of Eckert and Drewitz and that 


The uniform trend of the recovery factor 


In order to obviate the use of a method of successive approxi- 
mations for the determination of Cr and thus reduce the time 
required for computation, an equation was sought which would 
present Cr as an explicit function of 14; and R. Assuming the 
fluid to be air with y = 1.40 and Ce = 0.88, the equation 


inferred from Tifford. 
from the stagnation and flat plate flows would by the present re- 


sult involve an increase of only 2 per cent in the recovery factor. 


9) Small increases in the recovery factor in this range are indicated 
by some experimental results and not by others, but the small in- Cr = 0.0631 e—9-172 Mi R—0.182 
‘ crease may account for a masking of the effect by experimental ; ‘ . } ; 
(6) yields values for Cr which agree with those obtained by the parent 
error. é nate : ; 
. A P . é , > vet > < > ¢ @ 5 
The decrease in the recovery factor for the separation flow is ex- pry within 4 per cent in the ranges 1 < Mi < 3, 4(10°) < 
08. 


pected in terms of the movement of the region of maximum dis- 
sipation away from the wall. The reduction is apparently rapid 
and is not reflected in a contrast of the Eckert and Drewitz re- Pe 
sults for 8 = —O0.14 and 6B = 0. Such a rapid decrease suggests 
the utilization of the effect for the indication of flow separation, 
and low values of the recovery factor have been observed in the Correction for Aspect Ratio of the Indicial Lift 
crements separation region in experimental results for air flow over a 1:3 Functions of Wagner and Kussner 
elliptic cylinder. 
Huces the The accuracy of the solution here presented has been examined R. H. Scanlan 

in detail, by the methods indicated in reference 5. By develop- National Research Fellow, Ecole Nationale Supérieure de |’Aéro- 

nautique, Paris, France 


February 8, 1952 


umber of 


ing gn41 and g, intoa Taylor series, the interpolation error can be 








\ > 
determined 
E, = (e1/12) (d‘g/dn‘) 
9 , ' — HAS ESTABLISHED the correction o(k) to be added 
: 7M. ° : ic , a . ‘ . ‘ > corre: 7 y > > _ ’ . ° e ° ° 7 
This error — ler t nba isa i one in the to Theodorsen’s two-dimensional circulation function C(k) 
Eckert-Drewitz solution and partially explains the different value when the transition is made from a two- to a three-dimensional 
m ad ine = 1.0. The over: . > se i i : . “ ricnase 
' plished obtained for B + | rhe overall error of the present solution is sinusoidal flutter study. Recently, Ashley? has produced a simi- 
= estimated to be below '/2 per cent. lar correction term og(k) to be used in the case of a swept wing, 
a se where Reissner’s correction is inapplicable. 
ed by a REFERENCES 
Dr wit D.R.O ) Fall 1S) ' With either of these corrections ¢, we have the basis for the 
Cwilz Hartree, D. R., On an Equation Occurring in Falkner and Skan’s Approxi , : . a a — . a 
mate Treatment of the Equation of the Boundary Layer, Proc. Camb. Phil. Soc., correction of t he well-known indic ial lift functions k(s) of Wagner 
Vol. 33, pp. 223-239, 1937 and k.(s) of Kiissner for aspect ratio. It must be understood that, 
l() ‘Eckert, E., and Drewitz, O., Calculation of the Temperature Distribution for use in correcting these latter functions, the functions ¢, which 
Ids “ rs oe son pag ner of Unheated Bodies in a Flow at High Velocity, are valid only for harmonic motion, must be prepared separately 
ulti. Forseh., Vo 9, 1942 a 2 Piss 
Schuh, H., Laminar Heat Transfer in Boundary Layers at High Velaci- for each component (flexion, torsion, etc.) of each mode of wing 
1] lies, FB 2019 (1944): ZWB, R. and T. 810 (1947); MAP Volkenrode response envisaged and, with each mode, for a complete range of 
‘Tifford, A. N., Some Prandil Nuniber Effects on the Transfer of Heat, reduced frequency k = bw/v: 0 Sk S ~. Account must, of 
engines: rE . . St: News. The O State U “rsitv. Fe oa ° . 
; ag Experiment Station News, The Ohio State University, Febru course, be taken of the geometry of the specific wing studied 
ary, Yo 
Levy, S., Heal Transfer to Constant Property Laminar Boundary Layer Let the complex o be written 
12 Flous with Power Function Free Stream Velocity and Wall Temperature Vari — b ; 
ation, University of California, Institute of Engineering Research, Series o(k) = on(k) + toi(k) (1) 
$1, Issue 5, 195 . ° 
eons, BOSS whence the corrected C(k) is 
C(k) + of k) = F(R) + orn(k) + 7[G(k) + a7(k)] (2) 
= > —— — . + . . a . . 
) ; Now Garrick,’ for example, has established the relations 
1), and 
ce 2 * FS)... 
1es zero Comments on ‘Turbulent Boundary-Layer ki(s) = , sin ks dk 
Phe 7 J0 
to Characteristics...’ “ (3) 
q. (12 14 2 G(k) ks dk 
- = COS RS GR 
aber of Marvin L. Luther 7 JO k 
Research Department, U.S. Naval Ordnance Test Station, Inyokern ' 
iat : ita ' which relate ki(s) to C(k) = F(k) + iG(k). Garrick‘ has also 
China Lake, Calif. ; : 
January 23, 1952 pointed out the relation that must hold between h2(s) and 
(14 ki(s): 
I’ THE PAPER “Turbulent Boundary-Layer Characteristics at 1 , | a 1 
* ° * n> " . , / 5 = 5 + a — 
Supersonic Speeds—Theory and Experiment” (Journat or  *:(5) = se | ki(si) \> eae aE ds, e V se 4 
THE AERONAUTICAL SCIENCES, Vol. 17, No. 9, p. 585, September, (0 és (4) 
= . rm ‘ = . siti § > Z) 
1950), R. E. Wilson obtains the equation [Eq. (32), p. 588] ' 
] i ie 
* i 1 | 9 
sn * Vo _y¥-i!1 2 ().242 = ki(s1) . ds, (s > 2) 
(: + Cp M,? = wr Js-2 2-—-s4+5 j 
Ve 2 V Cr rae ' 9 
| The linearity of the entire theory permits the following simple 
. am : . - - . . " ~ ‘ : . 
logio (CrR) — 0.768 logio (: + Cp . us) generalizations, which define [ki(s)]an. and [ho(s)]a.R, as modi- 
laver P fied for aspect ratio for a particular wing and specified mode of 
i where at wing: 
actor : that wing 
p= o = [(y — 1)/2]M2/i1 + [(y — 1)/2]4n?} 2 ” F(k) + o(k) . 7 
i - [Ai(s) lar. = sin ks dk (5) 
which relates skin friction, Cr, to free-stream Mach Number r Jo k 
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1 j 
ki(s)lar. = [Ri(s)], | ds; + 
[ko(s) ]ar a. Ja. ES <a 
1 / 
/ (2 — .5 <2 
. V s(2 s) ‘o <'¢ = 2 (6) 
1 § | | ‘ 
= a = [Ri(s) lar \> a ds, 
(s > 2) 


The above aspect-ratio corrections are made in the same sense 
as those of Reissner—that is, for wake effects only, not account- 
ing for changes in apparent mass effects. Such apparent mass 
effects are represented, for example, in the uncorrected term 
(1/r) V s(2 — s)above. The resulting corrected terms may be 
employed for general unsteady three-dimensional analysis in the 
manner suggested, for example, by Baird and Kelley® for two 
dimensions, the usual Theodorsen lift circulation terms now being 
replaced by the following terms containing superposition integrals 
of the indicial lift functions: 


(1) Lift P for an arbitrary vertical movement of the wing 
defined by a vertical velocity function w(s) at the three-quarter 
chord point: 

P . 

= w(0) [ki(s) lar. t [ki(s Si) lar 
2rpvb 0 ds, 


w(s;) ds; 
(7) 
(2) Lift P resulting from the penetration of a gust of arbitrary 
vertical velocity distribution w(s): 


P 
= w(0) [Ro(s) lar. + [Ra’(s 


2rpub 0 


where ko'(s) = dko(s)/ds. 


x 


— Si)Jar. wisi) ds; ( 


Bisplinghoff, Isakson, and O’Brien® have already pointed out 
the application of Eq. (3) in the case of the generalization to their 
indicial counterparts of functions of harmonic oscillatory lift at 
high subsonic speeds in two-dimensional flow. 
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Design of Nozzles Having Continuous Wall 
Curvature 


J.C. Crown 

Hyperballistics Division, Naval Ordnance Laboratory, White Oak, 
Silver Spring, Md. 

February 1, 1952 


M*** NOZZLES FOR SUPERSONIC WIND TUNNELS are now de- 
signed using the Foelsch method! with a cubic initial 
curve.? Nozzle contours obtained from the Foelsch method have 
inherently a discontinuity in curvature at the ‘‘inflexion’’ point. 
Since a flexible plate cannot be bent to follow a nozzle contour 
having a discontinuity in curvature, it seemed desirable to de- 
velop a method for designing nozzles having continuous wall 
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Schematic illustrating notation 


curvature at all points. The basic formulas for designing such 
a nozzle are presented herein. * 


Symbols not shown in Fig. | are as follows: 


y = expansion angle 
c¢ = local velocity of sound 
p = density 


Along the Mach line from the inflexion point we have the char. 
acteristic condition 
where y; is the value of y in the test rhombus. Furthermore, 
dy/ds = 


pe 9 


\dt/ds = cos (0 a) 3 


ldy ds = sin (0 a 4 


Eliminating s and integrating, we have 


cos (a@ 6) " . 
&—-et= dy = Yo F cos (a — @)dy 
ws @ pe SJ 
and 
” sin (a — 6) “ — 
7 = dy = W F sin (a 6) dy (6 
0 pe J 0 
where 
fl + (Cy — 19/2] 2) 241 
F =< 
(y + 1)/2 { 
and y = ¥/Pwan; Y = Oat centerline, y = 1 at wall. The co 


ordinates &, 7 of said Mach line as functions of y can now be found 
if 0 = 6(Y) were known. This @ distribution will be further dis 
cussed below. 

Knowing the coordinates of the inflexion-point Mach line, we 
can find the coordinates of the downstream nozzle contour 
Since Mach lines of one family are straight downstream from th 


inflexion-point Mach line, it follows that 
x=€& + 1 cos (0 + a) ‘ 


y =n +/sin (6 + a@) 8 


where / follows from the equation of continuity 


L 
i Ywall ¥ = yF(1 J (9 


pe 
From the preceding equations it can be shown that for the wall 
curvature to be continuous, d@/dy must vanish at the inflexion 
point & = x, and — = &, Furthermore, for the gradient of cur 
vature to be finite, d?0/dy? - 
are satisfied by adopting a distribution along the inflexion-point 
Mach line 


© at these points. These conditions 


6= 0, sin? (4 2)y (10 


for example. The integrals (5) and (6) can now be evaluated 
with the aid of Eq. (1) and the usual expressions for vy and a 

the writer became acquainted 
Puckett 


* After the completion of this analysis 
through private channels with another method, developed by A. E 
to achieve continuous wall curvature 





The 
jt seen 
nozzle 
be see! 
10) al 
The si 
of ordi 


where 


In sl 

l 
line— 

9 

3 
upper 

(4 
tour 

(5 


or by | 


Foe 

a Paral 

March 

2 Cre 
1948 


Furt 


Marti 
Instruc 
Me 


Janua 


é b 
il 
Lees 
ackn 
Dries 
law < 
and 

by tl 
prob: 
pract 
amot 
sum] 
than 
and 

indic 
assul 
As a 


bines 





Ng such 


le char- 


> CO- 
yund 


dis- 


our 
the 


(9 


all 


ion 


ns 
int 





READERS’ 


The inflexion point is thus now a true one with zero curvature. 
it seems logical, as indicated in the remark to follow, that the 
R 


throat should also have zero curvature. This can best 


ale 
,. perhaps, by starting with a distribution such as Eq. 
10) and working backwards using the method of characteristics. 
The simplest upstream contour satisfying the imposed conditions 
of ordinate, slope, and curvature is then a quartic given by 


x # 
Yo + X, tan A ] - 
. Xx) 2x1 


Yo) cot 4 


(11) 


d ” 


where 

2(y1 (12) 
In short, the method is then as follows: 

1) Start with a @(y) distribution along inflexion-point Mach 


line—e.g., Eq. (10). 
2) Integrate to obtain £, » (coordinates of said Mach line) 


(3) Coordinates of inflexion point are found taking y = 1 as 


upper limits of integrals. 


(4) Calculate x, y (coordinates of downstream nozzle con- 
tour ) 
(5) Calculate upstream nozzle contour using either a quartic 


or by working backwards using method of characteristics. 
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Further Comments on ‘‘The Effect of Surface 
Cooling on Laminar Boundary-Layer 
Stability”’ 


Martin Bloom 
Instructor, Department of Aeronautical Engineering and Applied 
Mechanics, Polytechnic Institute of Brooklyn 


January 24,1952 


— WRITER, IN A RECENT NOTE, ! gave tworeasons for the errors 
in the stability calculations with heat transfer presented by 
Lees.2. One of these, the presence of numerical errors, has been 
acknowledged by Lees? and also seems to be implied by Van 
Driest,{ who lists differences in Prandtl Number and viscosity 
law as possible reasons for the discrepancy between his results 
and the results of Lees. However, the second reason given 
by the writer, and not mentioned by Lees or Van Driest, is 
probably more important from the standpoint of the convenient 
practical application of Lees’ analysis in cases involving large 
amounts of cooling at high Mach Numbers—that is, Lees’s as 
sumption that A [defined in Eq. (15), reference 2] is much smaller 
This device for approximately satisfying Eqs. (19) 
As the writer 


than 1 
and (20) of reference 2 bears further discussion 
indicates in work to be published by the N.A.C.A., the above 
assumption is unnecessary and in some instances is incorrect. 
As an alternative, Eqs. (19) and (20) of reference 2 may be com- 
bined to give the following expression : 


v{l (24, L)A + (@,? + ;7)d?/(1 + A)] = 


359 


FORUM 


With the values @; = 0.580 and #, = 1.48, corresponding to the 


minimum critical Reynolds Number, Eq. (1) assumes the form 


v{1 1.96 + 2.53d7/(1 + A)] = 0.580 (2) 
This equation replaces Eq. (64) of refererfce 2 with no further com- 
plication and holds for all values of \. For \ < 1, it is equiva- 
lent to Eq. (64) of reference 2. The iteration procedure sug- 
gested in connection with Eqs. (28) and (29) of reference 2 is 
thereby obviated. This iteration, moreover, is often inconven 
ient when \ < 1, since the second approximation for 4, often ex 
Trial-and-error 


Van Driest did 


ceeds 0.580 and, hence, one cannot proceed 
methods may, of course, be applied in this event. 
not indicate how he dealt with this problem. 

The values of \ which will be encountered when large amounts 
of cooling are required at high Mach Numbers can be judged 
from the following expression of Appendix C, reference 2: 

A = (0.332/7)) (n/c) l 

For significant amounts of stabilization, the factor n/c will not 
vary greatly and is of order one. However, the nondimensional 
ized surface temperature 7; may be small in these instances, so 
that \ may be large. 

Some values of \ encountered by the writer are given in Table 
1. These values are based on the assumption of a Prandtl Num 


ber of 1 and a linear viscosity law. 


TABLE | 
Ratio of Surface to 


Free-Stream Free-Stream 


Mach number Absolute Temperature (7)) r 
3.1 1.4 0.29 
4.4 0.93 Complete [.3 
5.6 0.11 stabilization 18 
5.7 0.008 240 
6 1.0 1.8 
6 0.05 43 


It may be remarked, in addition, that Figs. 7 and 8 and the 
tabulations in Appendix C and E (right-hand column) of refer- 
ence 2 must be corrected. Moreover, in view of Lin’s recent 
observations® that the stability theory developed in references 
6 and 2 is applicable, in a strict sense, only to cases involving 
moderately low Mach Numbers, results obtained by using the 
theory in its present form for higher Mach Numbers should be 


considered only for possible qualitative trends 
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